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ABSTRACT 



Composite structure repair methodology has been developed to 
specific applications (typically in small area and limited to secondary 
structure) and is being extended to Large Area Composite Structure 
Repair (with target extension to primary structures). Therefore, the 
repair becomes more critical because we get redistribution of stresses 
that can also affect the zones outside of the repair area. For this reason, an 
analytic evaluation of the repair’s reliability has to be performed to 
define a parameter which reflects on the effectiveness of the repair. 

In this work, we establish a principal guideline to evaluate the 
redundancy and compare the reliability of the repair to the reliability 
of the parent structure (i.e., the structure in the undamaged state). 
The approach adopted is to utilize structural finite element analysis to 
compute the state of stress at all the spatial elements of structure of 
the damaged state and of the candidate repaired state. The reliability 
of these two spatially non-uniform stresses is computed by a probabil- 
istic failure criterion. Thus, we can optimize the repair configuration by 
varying the strength and the stiffness of any element in the repair site 
by varying lamination angles, and selectively using hybrid materials. 
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I. INTRODUCTION 



A. BACKGROUND 

Composite materials are being utilized for commercial and mili- 
tary applications. Their combination of superior mechanical properties 
(especially the redundancy), design versatility, and low mass density 
make them attractive structural materials. Principal applications of 
composites have been on aircraft, especially on fixed-wing aircraft 
(e.g., F-18, F- 16, F- 15, F-14). 

As the structural use of composite material on aircraft increases, 
procedures to repair damage become essential. Composite repair pro- 
grams have been undertaken both by commercial and military sources 
[Ref. 1]. Repair concepts have been developed which include bolted 
patches repairs and bonded patches repairs. At the beginning, the 
repair technology was targeted for specific applications, typically lim- 
ited to small areas in secondary structures. More recently, repair 
techniques have been extended to larger areas, with target extension 
to primary structures (which inherently have high strength and dura- 
bility requirements). The particular techniques developed to repair a 
given component are dependent on the availability of specific facilities 
and the level of structural integrity. 

Under Reference 1, repair procedures were developed for two 
repair configurations. The configurations that were considered are a 
scarf inclusion and an external patch. A portable repair kit was 
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developed to provide vacuum pressure and heat in an integral rubber 
blanket. These techniques permit repairs, verified by testing, that are 
capable of restoring between 80 percent and 100 percent of the par- 
ent laminate ultimate allowable loading. 

B. STRUCTURAL REPAIR DESIGN 

Repair methodology development through testing has several 
inherent limitations. Only a specific structural configuration can be 
examined and optimized for variation of repair configuration in terms 
of the resultant effectiveness. A large number of experiments are fre- 
quently impractical due to economic and time constraints. This calls 
for the establishment of an analytical principle to guide repair 
optimization. 

In designing repair of composite structures, the external loading 
is known and the objective is to restore the structural functionality 
while recovering the overall structural reliability, which includes both 
the repair and the parent structure. The stiffness of the repair seg- 
ment and the strength of its attachment may both influence the overall 
structural reliability by altering the state of stresses in the parent 
structure. Thus, a repair segment with high stiffness and strong 
attachment may cause structural failure elsew r here, resulting in a 
strong repair but a weaker overall structure. On the other hand, a 
repair with proper stiffness but with insufficient attachment strength 
may fail in the attachment, also resulting in a low r structural reliability. 

For these reasons, an appropriate repair methodology for restor- 
ing the functionality and recovering the maximum reliability has to be 
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identified. The joint reliability of the structure— the repair segment, 
the attachment, and the parent structure— must be simultaneously 
considered by computing the reliability of each part associated with its 
respective stress states. The joint reliability of the structure can be 
used as a parameter which characterizes the overall effectiveness of 
the repair. Based on this parameter, any general configuration of 
structure and repair can now be evaluated and optimized. 

This investigation identifies and outlines the principal steps asso- 
ciated with such an analytical optimization of general structural repair. 
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II. STRUCTURAL REPAIR APPROACH 



A. GENERAL CONCEPT 

A vital attribute of composite materials is structural redundancy. 
In order to identify the appropriate repair methodology for composite 
structures, the underlying redundancy must be considered. 

The objective of repair is to restore the redundancy, even at the 
expense of increasing the stresses. In other words, our goal is to 
restore the functionality and to maximize the recoverable reliability of 
the whole structure. The effect of the repair on the entire structure 
needs to be considered (i.e., the redistribution of the stresses caused 
by the damage, and later the repair itself which may transfer the fail- 
ure to a different location). 

B. FUNCTIONALITY AND RELIABILITY 

To restore the functionality, elements such as the type of aircraft 
component and aerodynamic (smoothness) need to be examined. 

From the strength theory for composite materials, the main fac- 
tors that are peculiar to laminate strength analysis (neglecting the 
thermal effects) are: 

• laminae strength 

• laminae stiffness 

• laminae orientation 

• laminae thickness 

• stacking sequence 
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Thus, for restoring the reliability (redundancy), we need to opti- 
mize the patch stiffness through the materials selection and the lay-up 
orientation (including stacking sequence), together with different lap- 
joint configurations (scarf, step, or external patch). Moreover, the 
optimization objective is to achieve equal reliability in the parent 
structure, the joint (the adhesive), and the repair patch, thereby 
achieving the highest possible total reliability for the entire structure. 
Upon the development of analytical computation methods for the joint 
reliability of the entire structure, the optimal repair configuration may 
be analytically sought by “iterating the choice of material (including 
hybrids), the orientations of the repair patch, and the joint 
configurations.” Figure 2.1 is an example. 

A composite material structure with a repaired area is subjected 
to a complex loading, as shown in Figure 2.1. To evaluate the reliabil- 
ity, the stresses in each spatial location need to be determined 
through the stiffness and material iterations. 

This example is a boundary value problem, and for that the linear 
elastic field equations must be satisfied (see Ref. 2). 

Due to the external load (Tj-surface traction vector), the stress 
boundary' equation is considered, 

i, j = x, y, z (2.1) 

wdiere (ay) is the normal and shear stress components and (vj) is the 
unit outer normal vector to the surface. 
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«,y 



(b) 



(c) 



Figure 2.1. Typical Repaired Structure 
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Due to the kinematics, the strain-displacement relations are 
defined (for small displacements): 

e ij = i(d i u j + 5 ^) i ( j =xyz (2.2) 

where (ey) is the strain components and ui is the displacement. 
Substituting the strains into the constitutive relation, 



O =C,., £. 

lj Ijlm lm 



i, j, 1, m = x, y, z (2.3) 



where cyim is the stiffness matrix. 

Hence, the resulting stresses are determined by satisfying the 
equilibrium equation. 



3.c (i + X. = 0 

> 1J J 



i, j = x, y, z 



(2.4) 



where Xj is the body force (usually ignored in the analysis). 

Now. recall from the constitutive relation (equation 2.3) that the 
stiffness matrix (Cyim) characterizes the material property. In contrast 
to isotropic material, where the stiffness consists of two independent 
constants, and more importantly, the stresses are continuous, in an 
orthotropic composite material the stiffness matrix consists of nine 
independent constants, which reduce to five for special orthotropic 
orientation (because with respect to material properties Z and Y axes 
are equivalent). In a laminate, when the material principal axes are off 
the structural axes and there is symmetry'’ about z = 0, the stiffness 
matrix consists of 13 independent constants. Moreover, the stresses 
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are continuous across the thickness (t), and the stiffness is a function 
of the lamina orientation and the stacking sequence. 



For a thin lamina under plane stress, the constitutive relation can 
be expressed as (Ref. 2]: 

The stress-strain relation is 



where the Qy, the so-called reduced stiffnesses, are the Cjj] m stiff- 
nesses specialized for plane stress condition in an orthotropic mate- 
rial, and written in contracted notation (i, j = 1, 2, 6). 

The stress-strain relation can be related to the k th lamina, and 
later (for FEM) to the n th element, 



From the laminate plate theory [Ref. 2], the resultant external 
force acting on the laminate shown in Figure 2.2 is obtained by inte- 
gration of the stress in each lamina through the laminate thickness. 
For example. 




(2.5) 



{».} =w 

k K k 



i, j = x, y, z r, s = 1, 2, 6 



( 2 . 6 ) 
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(2.7) 
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where N x is the force per unit length of the cross-section of the lami- 
nate. Similarly, M x is a moment per unit length. 




Figure 2.2. In-Plane Forces and Moments 



Finally, by substituting equations 2.5 and 2.2 into equation 2.7, 
and summing the forces (Ny) and the moments (My) for the whole 
structure, the following relations are derived [Ref. 2], 

{N # } = [A„]{e°,} + [ B 

{M, J } = [B„]{e»} + [D„]{k 1J } ( 28) 

where 

N 

A rs = 2f(Qrs) ( Z k _Z k-]) 

k=l K 



B rs = |X(9rs) (Z 2 k "ZD 

k.l k 



D rs = |S(Qrs) (z 3 u -z k _l) 

k=l k 
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a w 0 




(2.9) 



£° and ky are middle surface values (not functions of z), £° is the 
strain component, ky is the curvature, and wo is the displacement of 
the middle surface in the z direction, t is the thickness and Zk is 
defined in Figure 2.1(c). In equations 2.8 and 2.9, A rs is called the 
extensional stiffness, B rs is called coupling stiffness (B rs = 0 if the 
laminate is symmetrical), and D r s is called the bending stiffness. 

Therefore, the stiffness in composite material is not dependent 
only on the material stiffness (Qrs) but also on the laminae thickness 
(t = Zk - Zk - 1) and the stacking sequence (zk). The stiffness matrix 
(Qrs) changes due to the orientation of the laminae as it corresponds 
to the structure coordinates. 

In short, to iterate the material stiffness for determining stresses, 
numerous independent parameters are involved. The parent structure 
configuration is given (the stiffnesses are known). But in addition, in 
composite repair, two configurations are concurrently being consid- 
ered— the parent structure and the repair (shown in Figure 2.1). For 
each of those configurations, the stresses must be evaluated and by the 
constitutive relation (equation 2.3) they must match in their bound- 
aries. Before repairing the damage, the parent structure with a "hole" 
configuration needs to be iterated, and then the repaired structure 
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configuration. Through the resulting spatial non-homogeneous 
stresses, it is impractical to a priori determine which configuration is 
optimal (i.e., as can be done for isotropic material). Thus, by calculat- 
ing the total reliability from the state of stresses, we unify all the iden- 
tified parameters into one single parameter for optimizing the repair 
configuration. 

C. CONSIDERATION 

In this investigation, we assume that the functionality of a 
repaired structure must be the same as before damage, so the repair 
effectiveness is measured by the post-repair reliability. To evaluate the 
reliability, comparisons of the structure before damage, before repair, 
and after the repair need to be done (see Figure 2.3). 



Parent 

structure 



Damaged Repaired 

structure structure 




R'>R"<<R' 



Figure 2.3. Repair Configurations 



The merits of the different elements of design of the repair can 
be consolidated into one single unified parameter and measured by the 
structural reliability. 

From the general concept of reliability, the number of probable 
failures is associated with the magnitude of the stress level (the ran- 
dom variable). For a structure of given dimensions, the higher the stress 
level the higher the probability of failure. The quantitative probability of 
failure can be computed if the cumulative distribution of failure is 
known. For a structure under non-uniform stress (spatial heteroge- 
neous stress), the structure may be visualized as being divided into 
contiguous geometric elements within which the stresses are uniform. 
The reliability of each element corresponds to the respective stress 
magnitude of that element. Equivalently, the magnitude of the stress at 
each spatial location of the structure can be thought of as being 
mapped into a respective point on the stress space (by appropriate 
stress analysis). If all the stress magnitudes are within the domain 
bounded by the strengths, then the entire structure is considered to 
be safe, as shown in Figure 2.4. 

The mapping of the physical spatial location to the stress domain 
is carried out in stress analysis. The definition of the strength in the 
stress domain is a failure criterion which can be defined either de- 
terministically or probabilistically. 
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1. A deterministic failure criterion can expressed in terms of the 
scalar stress components, which are themselves a function of the 
spatial location of the structure. 



f ( a ij ( x k)) <1 safe 

f ( a ij( x k)) >1 fail (2.10) 



where Gy is the stress component in tensor notation (i = 1, 2, 3) 
and Xk is the spatial coordinates within the structure. F(.) is a 
function representing the contour separating all the safe stress 
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states from all the non-safe states. For example, for coupled 
strength. F(.) may be expressed in terms of a polynomial crite- 
rion. as described in Reference 3. 

2. A probabilistic failure criterion expressed in terms of scalar 
stress components as random variables. The magnitudes of the 
scalar stress are themselves a function of the spatial location of 
the structure. 



where R is the reliability of a physical element located at spatial 
coordinates Xk within the structure, under stress components 
(cfjj)k. «k and (3k are the statistical parameters of the probabilistic 
failure model associated with the component stress ay. The 
entire structure is safe if the joint reliability for assemblage of all 
the spatial locations of the structure is greater than or equal to 
the reliability specification. 

D. TRADE OFF 

In the spirit of the consideration above, the joint reliability is 
consorted measure of the all the parameters wdiich need to be opti- 
mized for a given repair. To explore the effects of the repair, w T hich 
include the size of the damage area (after the clean-up) and the patch 
type (material, stiffness, and type of joint), trade-offs between stiffness 
and strength of the specific parent structure, the joining segment, and 
the repair patch need to be considered. 

If, w r hen the repair design is completed, the structural reliability 
is not met, the design parameters need to be iterated to an achieve an 
acceptable solution. This iteration involves the redesign of the struc- 
ture/ repair configuration and the selection of new material and/or a 
new^ lay-up configuration. 




( 2 . 11 ) 
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III. STRUCTURAL RELIABILITY (CONCEPT) 



A. GENERAL CONCEPT 

In Chapter II, we developed the repair design concept as it 
relates to reliability. In this chapter, we will explore the specifics of 
stresses and strength calculations and examine their relationship to 
reliability. This relation allows us to compute reliability based upon the 
failure criterion which is most appropriate to the specific composite 
under consideration. 

Historically, formulation of many failure criteria was intuitively 
taken to be mechanistically based on failure mechanisms, but in fact 
the criteria are phenomenological (because frequently, the plane of 
failure does not necessarily correspond to applied stress). Phenome- 
nological failure criteria represent (describe) failure states but do not 
explain failure causes (or mechanism). There is no uniquely “correct” 
formulation, and there is no a priori known formula. 

The failure modes observed for composite materials are quite 
unlike those observed for isotropic homogeneous materials because of 
their anisotropic and heterogeneous nature. Nevertheless, it is essen- 
tial for design to predict the capability of the structural element under 
a complex loading. Judgment must be made as to whether the theory 
(analytical results) adequately represents (for a wide class of problems) 
the physical phenomenon (i.e., experimental results). 
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In the design of repair, we are looking at the parent structure 
characteristic and the repair characteristic as a whole system’s ability 
to withstand a given load. Then, by the reliability characterization w r e 
determine the stress level (coming from the applied load), for which 
there is an acceptable number of probable failures. 

In a composite material, the lamina is the basic building block 
which consists of fibers imbedded in a matrix. In composite structure, 
especially in repair areas, the stresses are spatially heterogeneous, 
both in forms of stress components and in their magnitudes. More- 
over, because the lamina is composed of numerous fibers and 
possesses inherent redundancy, the composite failure process is 
sequential (in stress level). At the same time, the actual behavior of the 
loaded structure often exhibits considerable data scatter in measured 
strength due to processing variability' (u'hich is unavoidable in large 
composite structures). 

(Anisotropic) failure criteria are mathematical functions (analyti- 
cal. numerical, graphical) which define the stress levels below' w'hich a 
structure is safe. If the measure of safety' is absolute, the failure criteria 
are deterministic: otherwise, the failure criteria are probabilistic. The 
magnitude of stress (one component for 1-D; six components for 3-D) 
at each spatial location of the structure is mapped into a respective 
point in the stress space. If all the points are within the domain 
bounded by the strengths, then the entire structure is safe. Ideally, a 
failure criterion should capture the physics of the specific failure pro- 
cess and the sequential events of failure. 
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Therefore, in order to calculate the reliability of a structure with 
non-uniform stress, the state of stress at every spatial location needs 
to be computed. The influence of these stresses on the respective 
reliability in those specific spatial locations can then be evaluated by an 
appropriate probability model. 

B. FAILURE MODEL- WEIBULL DISTRIBUTION 

In 1958(?), Rosen recognized that the mechanism of composite 
failure originated from the statistical variability of the fiber strength 
[Ref. 6], The failure process w r as observed to be sequential, initiating 
from failure of v r eak fiber sites; the load originally carried by the w^eak 
fibers is transfered to the adjacent strong fibers, thereby postponing 
immediate failure. Upon additional load increase, the local micro- 
redundancy is insufficient to sustain the stress concentration, leading 
to catastrophic macroscopic failure of the composite. Harlow and 
Phoenix [Ref. 4] refined the probabilistic modeling of sequential failure 
events and arrived at the “chain-of-bundle” model. The composite 
structure is modeled as a chain of bundles in series. Each bundle has 
length 5, generally referred to as the “ineffective length,” The ineffec- 
tive length is a function of the matrix stiffness and the interfacial 
bound strength; it is analytically estimated to be an order of magnitude 
greater than the fiber diameter (the diameter of graphite fiber is 
around 5 pm). Each bundle is assumed to be independent and identi- 
cally distributed, and the cumulative distribution function (CDF) can be 
expressed in terms of the fiber statistics. Harlow and Phoenix have 
also derived, by analytical recursive formulation, the statistical 
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strength behavior of large composite from small bundles. Since the 
bundles are in series, the weakest link function is applied. For fila- 
ment strength, the Weibull model is appropriate for the characteriza- 
tion of the strength statistics. For the composite, the strength may be 
considered as piece-wise Weibull, where the scale and shape parame- 
ters are modified by the Harlow-Phoenix chain of bundle calculations. 

For this investigation, we consider each spatial element as 
independent and ignore failure sequence and subsequent load 
distribution. This is equivalent to stating that if one point in any 
direction fails, the whole structure fails (failure can occur in any element 
of the structure, in either the parent structure or in the repaired 
portion of the structure). The reliability under this condition can be 
calculated given the failure distribution and the failure criterion. 

We use the Weibull distribution with parameters adjusted by the 
‘‘chain-of bundle” model. The two-parameter Weibull distribution is 
expressed as: 

For any i th spatial location (considered in our calculation as an ele- 
ment), the failure distribution function (CDF) is: 



Si— stress component at location Xj 

Pi— the location parameter, representing the strength under the 

respective stress component (oj. It is dependent on the size of 
the element. 




( 3 . 1 ) 



where 
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aj— the shape parameter, representing dispersion of the failure 
function. It is also dependent on the size of the element. 

Then, the reliability at each point, o,(x^ 



R(°i) = 1 " F (o ) :(x j )p 1 ,a 1 ) 



(3.2) 



R (°0 =exp 




(3.3) 



Assuming that the failure (Pi) event is independent of all spatial 
locations, the total structure reliability (Rt) can be represented by: 



R t =R 1 R 2 • • • R k 



R=fl R v 



where ( r k) is the reliability at each spatial location. 



(3.4) 



C. SIZE— EFFECT 

In our approach, for calculating the state of stress we are using a 
finite element method (FEM). For that we assume that a configuration 
of repaired structure is partitioned into geometric meshes n of size I 
and m of size 2. The intrinsic strength of each element may be differ- 
ent. The choice of element size is arbitrary, and because the same 
structure is under consideration, the total reliability must remain the 
same. 
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( 3 . 5 ) 




where 

Rt~ the total reliability 

Ri, R 2 — The element reliability (assume equal for simplicity) of each 
meshing. 

n— number of element for size 1 
m— number of element for size 2 
By Weibull distribution, 



We make the simplifying assumption that the two different sizes 
are not drastically different so that the same shape parameters a can 
be used. It follows that: 



Therefore, we can evaluate the reliability of each element by 
relating its length or area to the unit element. In our approach, the 
reliability' is normalized by the smallest element. The smallest element 
is assumed to be the ineffective size. In short, we need to first perform 
stress analysis, which is in sections, then evaluate the reliability 
knowing the strength of the smallest element. 




( 3 . 6 ) 




( 3 . 7 ) 
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IV. REPAIR METHODOLOGY 



A. TYPICAL REPAIR REQUIREMENT 

The previous chapters discussed the rationale of adopting the 
joint reliability as a single indicator of the repair efficiency and the 
property of a probabilistic failure criterion which is required to calcu- 
late the reliability of each spatial element under a non-uniform 
complex state of stress. Furthermore, we identify the following 
requirements which are necessary to establish the appropriate repair 
methodology for composite structures with the goal of restoring the 
structural integrity: 

1. The entire structure needs to be considered. 

2. The patch stiffness through the materials selection and the lay-up 
orientation, together with different lap joint configuration, need 
to be optimized. The optimizing objective parameter is the joint 
reliability of the structure, including the parent structure, the 
joint, and the repair patch. The objective is not to maximize the 
reliability of the repair but rather to achieve highest possible total 
reliability for the entire structure. 

3. In order to assess the necessity or the effectiveness of the repair, 
a single parameter is formulated. We select reliability as the 
parameter (refer to Chapter III). 

4. Using the reliability formulation, the parent structure reliability is 
compared to the damaged structure reliability before repair, and 
the repaired structure reliability. As mentioned previously, the 
ideal repair restores the reliability of the repaired structure. 

5. It is known that conventional repair depends on the adhesive 
joint, which always has no redundancy or poor redundancy 
(smooth-joint versus step-joint). This weakness has to be taken 
into account. In our approach, the adhesive strength affects the 
stress distribution and limits the repair stiffness (see joint con- 
figuration in appendix). Thus, the joint strength is considered 
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with the stiffness optimization (lay-up design) of the repair seg- 
ment (i.e., weak adhesive limiting the transferring of stresses to 
the repair segment compliant patch) 

6. To simplify the calculation, and to reduce the number of parame- 
ters that have to be iterated in the same time, we assume the 
adhesive reliability to be equal to 1. In other words, we assume 
infinite joint reliability. 

7. The reliability of the structure is the product of the reliability of 
the individual elements (see equation 3.4). Moreover, to achieve 
the reliability of the element, it is required that the magnitude of 
the stress tensor not exceed the failure criterion (i.e., a spatial 
reliability function with stress tensor as the random variable (see 
Chapter II). 

B. CONFIGURATIONS FOR REPAIR 

Now, in order to utilize the appropriate procedure for designing a 
repair, we have to consider three separate configurations: 

1. The parent structure 

2. The damaged structure 

3. The repair 

I. Parent Structure Configuration 

For the parent structure, a finite element model must be built 
to represent the original geometry, the lamination configuration, and 
the original material. The boundary' conditions are determined by the 
actual external load and reactions (support). 

By knowing these data, the resulting stress tensor for each 
element can be calculated. Using the suitable failure criteria, the 
parent structure reliability (RO) can be evaluated. The parent structure 
reliability is the acceptable parameter reference to compare for 
optimizing the repair configuration. 
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2. Damage Structure Configuration 

We assume that all the facilities needed to evaluate the 
damage area are available (i.e., NDT). Thus, before removing the 
damage and preparing the area for the bonding repair, the clean-up 
geometry has to be optimized. Dependent on the size of damage, the 
spatial location, the repair capability, and the state of load in this area 
(considering the fibers’ orientation), several configurations can be 
evaluated. Again, going through the same procedure introduced for the 
parent structure, the state of stress is calculated and thereby the 
reliability for each configuration, Rl, R2, R3. Now the configuration 
with the largest reliability (say R3) is selected. 

3. Repair Configuration 

Starting with the R3 configuration, and assuming infinite 
joint reliability, the joint strength is considered with the optimized 
stiffness. Thus, only the reliability of the parent structure (outside of 
the repair zone) and the reliability of the repair structure (the patched 
zone) are considered for calculating the total repaired structure relia- 
bility. To optimize the repair configuration in the lay-up design, one 
must control the repair stiffness and redistribution of the stresses. 
Either one type of material can be chosen and only the lamination can 
van,' (lamina orientation) or more than one material and lamination can 
be applied. 

C. FINITE ELEMENT APPLICATION 

It has been showm that to calculate the reliability, w r e need to ana- 
lyze the state of stress in all spatial locations— analytical boundary value 
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problems or finite element methods (FEM) can be used. Because of the 
complexity of the structure geometry, it is essential to use FEM (refer 
to ADINA code in Appendix A). 

The analysis of stress in a laminate of laminae is straightforward. 
For 2-D problems, we are using the laminate plate theory (plane 
stress). For 3-D problems, we are using generalized calculations for 
orthotropic materials. Assuming the laminae exhibits linear elastic 
behavior, the stresses seen to be a linear function of the applied loads. 
That is, if all laminate stresses in each laminae (i.e., corresponding to 
an element in ADINA code) are known, then the stresses in each 
laminae can be compared with the lamina failure criterion (see Chap- 
ter III) and uniformly scaled to determine the load at which failure 
occurs and the spatial location. 

In accordance with the reliability and failure criteria, all failure 
criteria for laminates depend on the strength in the principal material 
direction, which likely does not coincide with principal stress direc- 
tions. Therefore, the strength in each lamina in a laminate must be 
assessed in a coordinate system that is likely different from those of 
its neighboring laminae. (This coordinate mismatch is one of the com- 
plications that characterizes the determination of the failure criteria.) 

Hence, for the interpretation of stresses resulting from the FEM 
calculation, all the stress tensors have to be transformed to the prin- 
cipal material direction per lamina (see Appendix C). 

Furthermore, in order to normalize the reliability and to make it 
geometrically independent, the area of each element has to be related 
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to an acceptable minimum area. The minimum area is derived from 
the minimum size effect. 
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V. METHODOLOGY FORMULATION 



The object of this chapter is to explore the effects of the damage, 
and finally the repair, on the stresses of the entire structure. Appro- 
priate simplifying assumptions have been represented in the the 
repair methodology (Chapter II). 

We wish to capture the essence of the underlying mechanics and 
physics of our problem. The field equations in solid mechanics are 
specialized for calculating the state of stress in the repaired structure. 
Then, the acceptable probabilistic failure criteria, expressed as a 
mathematical function which defines the stress level below which the 
structure is safe, is utilized. 

Therefore, in the context of preliminary sizing and design trade- 
offs of the repair, and also for confirmation of output of more accurate 
analysis carried out by the computer code, we represent the repair 
methodology by the following formulations. 

The formulations are classified by : 

1. The physical space— a line, a plane, a volume 

2. The stress space— one component, three components, six 
components 

3. The spatial distribution of stress within the structure— homoge- 
neous or heterogeneous 

4. The definition of safety— for our case, probabilistic failure criteria 

Whereas in one-dimensional theory (1-D uniaxial stress), each 

scalar stress is mapped into a respective point on a line (1-D) in the 
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stress space, under planar combined stress, the magnitude of stress 
(three components) at each spatial location of the structure is mapped 
into a respective point on the 3-D stress space. Under general com- 
bined stress, the magnitude of stresses (six components) at each spa- 
tial location of the structure is mapped into a respective point on the 
6-D stress hyperspace. 

Stress interior to domain bounded by the failure surface — > safe. 

Stress exterior to domain bounded by the failure surface — > not safe. 



A. ONE-DIMENSIONAL EXAMPLE 

In order to understand the methodology principle, a simple stress 
analysis, consisting of a scalar calculation for single component stress, 
is shown in the following example. 

A long segment (shown in Figure 5.1) of non-uniform cross-sectional 
area (A(x)) is subjected to a uniaxial tension load (P). Then, to cal- 
culate the reliability: 

1. The long segment is divided to small and unequal elements, 
because of the usual complexity of the repair geometry (geometry 
mapping) and because of different stress magnitude distributors 
(i.e., the segment is spatial stress heterogeneous). 

2. The stresses (o(x)) are uniaxial stress (1-D) and must be evalu- 
ated in each spatial location. Thus, because of uniaxial stress (one 
component, 1-D), the stress analysis is simple and represented 

by 



o(x) = 



p 

A(x) 



(5.1) 
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Figure 5.1. Long Segment, Non-Uniform Stresses 



Therefore, because of the series process assumed for composite 
materia] structure (refer to Chapter II), the total reliability of the 
structure (Rt) is presented as a function of the stress and the 
element size. 



R t = n R ((^,(x)).Ax.) 



i=l, 2 ... n 



(5.2) 



where 

n— the number of elements 
Axj— is the i th element length 

Furthermore, using the Weibull distribution function (weakest 
link of chain), the total reliability is expressed as, 



R T = riexp 



(2i 

u, 



(5.3) 



3. To compare the structure configuration (i.e., for different load, 
different type of repair), and to avoid size effects on the strength 
of the material, the reliability must be geometry independent. 
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Thus, the reliability is normalized by the smallest length, as 
follows. 

Using the weakest link rule [Ref. 4], and assuming the shape 
function is constant (in composite material a = 5~10), the size 
effect is represented by, 



1 




(5.4) 




(5.5) 



where 

5— the length of the smallest element 

(3— the mean strength value of the smallest element (evaluated 
by experiments) 

(3) — the mean strength value of the i th element. 

Then, by substituting equation 5.5 into equation 5.3, the total 
reliability is obtained as: 




(5.6) 



4. Finally, a simulation needs to be performed to demonstrate the 
probability of different loading configurations. This simulation was 
carried out via a statistical analysis using the spreadsheet Micro- 
soft Excel on a Macintosh computer (see Reference 9). There, 
the strength is evaluated by the Weibull distribution function 
(weakest link), and if any stress is larger than the suitable 
strength, then the structure fails (1— fail, 0— no fail). 

Thus, after several runs of the simulation, we can observe that the 
location of highest stress has the highest probability of failure but 
not a certainty of failure. In other words, the failure is not con- 
stant due to the spatial location— it is changing according to the 
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probabilistic failure criteria. The resulting reliability distribution 
is compared so the configuration of the structure loading can be 
optimized. 



B. TWO-DIMENSIONAL EXAMPLE 

The first approach that corresponds to the macro-mechanical 
behavior of a composite laminate is the 2- D stress calculation using 
the classical lamination theory [Ref. 2]. Thereby, under the same prin- 
ciple as 1-D calculations, 2-D calculations are performed, considering 
three stress components: 



{CTj ,a 2 ,c 6 } 

A typical analysis of a thin plate with an elliptical hole subjected to 
tensile load in a non-principal direction is described. This analysis 
models the damaged area of a structure after “clean-up” and before 
the bonding process of the “patch” repair. 

Again, going through the same procedure as the 1-D example in 
order to calculate the total reliability, 

1. The physical space considered is the space area of the entire 
plate (laminate). The space is divided into unequal areas. The 
rectangular rather than the polar division is preferable because of 
the orthogonality of the composite material, and because it is 
easier to compare the structure reliability as a spreadsheet. 

2. Assuming an unequal area for each element, the stress analysis is 
performed. 

Note that for using reliability equations, the stresses in each spa- 
tial location must transform to the principal material direction, 
as explained in Chapter II. For this, stress transformation equa- 
tions are used [Ref. 2]. 

Now, having the stresses in each spatial location, the same pro- 
cedure as in the 1-D example is utilized. We must remember that. 
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here, the weakest link rule refers to each component of stress in 
each element (three components). Therefore, the weakest link 
rule refers to each element in the entire structure (plate). 

3. To normalize the reliability, the lengths in the previous example 
(equation 5.5), are replaced by areas: 



k— relates to the element notation, 
i— relates to the stress component (i=l,2,6). 

A m — the smallest element area. 

Ak,i— the k th element area. 

4. Therefore, because of the weakest link approach, there are two 
options to calculate the total reliability (Rt)- 

Option 1: Calculate the reliability of the whole structure with 
respect to each component in each element (shear or normal 
stress), then calculate the total reliability, as shown below. 

The reliability of the k th element due to the i lh stress compo- 
nent is 



Note that the stresses are in contract notation. 

The whole n element's reliability with respect to the i th stress 
component is: 




(5.7) 




(5.8) 



n 



R t.i = IK, k = 1 ... n 



(5.9) 



k=l 



Finally, the total reliability is: 
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R. 



=riR 



T.i 



i = 1, 2, 6 



(5.10) 



i=i 



6 n 



6 n 



R T =nnR k ,=nnex P 



1=1 k =1 



1=1 k=l 



f<5 



k.l 



^k.l 



i = 1. 2, 6 

k = 1, 2 ... n (5.11) 



Option 2: calculate the reliability of each element and then cal- 
culate the total reliability as follows. 

The reliability per element (RE,k) is the product of the stress 
component reliabilities (Rkj— expressed in equation 5.8). 

Thus, 



R 



E.k 



=riR 



k.l 



i = 1. 2, 6 



(5.12) 



1=1 



Then, the total reliability is. 



R T = flR 



E.k 



k =1 



(5.13) 



And, after substituting equations 5.10 and 5.14 into equation 
5.15, the total reliability is: 



1 l V 

r t =nriexp 



k = l i=l 




i = 1, 2 6 

k = 1,2 n (5.14) 



Although the total reliability must be the same for any option, a 
difference exists between the two options and is obtained by the 
interpretation of optimizing the repair. Through Option 1, we can 
visualize the effect of a specific stress in a specific direction (i.e., 
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improve fiber, if the total reliability— [Rx.j— fiber direction], is small). 
Or, through Option 2, the weakest element which has the lowest reli- 
ability (RE.k) can be investigated. 

C. ANALYSIS IN THREE DIMENSIONS 

The next obvious approach is to extend the analysis to a laminate 
in 3-D physical space considering combined and spatial heterogeneous 
stresses (six components). Obviously, the stress analysis is more 
complex and cannot be solved “by hand,”— a numerical analysis 
method must be used (see application of ADINA code in Appendix A). 
Moreover, following the same procedure that is used in the 1-D exam- 
ple and in the 2-D example, the total reliability is calculated. 

1. The state of stresses resulting from the 3-D analysis (generalized 
orthotropic laminate— Reference 2), must relate to the mean 
strength known for laminae (from experiments). In other words, 
the stress must be calculated per element in the lamina. There- 
fore, unit volume will have an identical thickness for each ele- 
ment. Hence, the area normalization in 2-D (see equation 5.7) is 
sufficient also for the 3-D calculation. 

2. Furthermore, to calculate the reliability, instead of three compo- 
nents in the 2-D example, we consider six stress components: 



For the first option mentioned in the 2-D example, equation 5.10 
extends to: 





i = 1, 2 ... 6 (5.15) 
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The whole n element reliability with respect to the i th stress 
remains the same (see equation 5.9). Hence, the total reliability 
is: 



6 n 



6 n 



R T =nnR k ,=nnex P 



i = 1 k = 1 



i = 1 k = l 






v 



Pk.i J _ 



i = 1, 2, 6 

k = 1,2, .... n (5.16) 



For the second option mentioned in the 2-D example, equation 
5.12 extends to a product of six instead of three, as shown below. 



R 



E. k 



= I1r 



k.i 



i = 1, 2 ... 6 



i=i 



(5.17) 



Thus, similarly to the 2-D example, the total reliability is: 



n 6 



r t =nriexp - 



SO. A 

k .i 



k =1 1 = 1 



0k., 



i = 1, 2 6 

k = 1, 2 n (5.18) 



D. SUMMARY 

1. The 1-D example is only a simple presentation of the methodol- 
ogy concept used for calculating the total reliability. 

2. In the 2-D example, the whole structure analysis is simplified 
using laminate plate theory. This approach is used in practice to 
do preliminary calculation and sizing, and to check the results of 
2-D and especially 3-D numerical analysis carried out by com- 
puter codes (see Appendices A and B). 

3. For the 3-D example, the state of stress in each spatial location, 
considering the generalized orthotropic laminate, is numerically 
calculated using FEM by computer a code (i.e., ADINA). 

4. Comparing 2-D to 3-D as the generalized orthotropic laminate 
analysis, in 2-D analysis the interlaminar stresses are ignored. 
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VI. CONCLUSIONS 



This work has proposed a primary approach for the designing of 
composite material repair, using reliability as a single parameter to 
indicate the most effective repair. Finite element analysis was used to 
calculate the state of stress at any spatial location (non-uniform and 
combined), which in turn was used to calculate the resultant reliabil- 
ity. This methodology can be used to assist maintenance personnel in 
the selection of the optimal repair configuration. 

The Weibull distribution function with modified shape and scale 
parameters was used for estimating the reliability, based on the obser- 
vation that composite material failure is adequately represented by the 
weakest link chain-of-bundle model. The best selection of failure cri- 
teria has to be confirmed by experiment. 

Simulations for different patches under a specific loading or for 
specific patches under different loading need to be included in the 
future extension of the repair optimization method. Moreover, an ade- 
quate number of samplings (experiments) needs to be performed in 
order to estimate the underlying Weibull parameters and then to 
approve the appropriate failure criteria. 

Only bonded repairs with infinitely strong adhesives have been 
considered. Further studies may be extended to apply this approach to 
consider adhesives with finite reliability and bolted repairs, which are 
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considered as 3-D reliability enhanced (small transverse fibers are 
included under this category), and to scarf joint geometry approval. 
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APPENDIX A 



DESCRIPTION OF ADINA 
A. GENERAL DESCRIPTION 

ADINA is a finite element code for automatic dynamic incremental 
nonlinear analysis. The ADINA system is composed of four main pro- 
grams: ADINA-IN, ADINA, ADINA-PLOT, and ADINAT. Version 2 (Jan. 
1989) was used for this thesis. The sequence of execution of these 
programs is shown in Figure A-l [Ref. 6:pp. 1-2]. 

For our investigation, the ADINA code was used to conduct the 
stress/strain analysis (3-D/2-D) of an advance composite material 
structure. Thermal effects were ignored, so the ADINAT program (for 
heat-transfer analysis) was not used. 

The ADINA-IN program generates the input data for the ADINA 
program and has the ability to graphically check the model. All or part 
of the model can be drawn from any view desired, including element 
and nodal numbering. 

Coordinates can be input in a local or global coordinate system 
using cartesian, cylindrical, or spherical coordinates. 

ADINA is used for displacement and stress/strain execution, pro- 
duces the printout file for displacement and stress/strain, and also 
produce the porthole file for ADINA-PLOT. 
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Figure A-l. Sequence of Execution of ADINA Programs 
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ADINA-PLOT performs the post-processing of the output data 
from ADINA. It also has graphics capability, allowing a view of both the 
deformed and undeformed shapes, and can do vector plotting of 
stresses. 

B. COMPOSITE MATERIAL IMPLEMENTATION 

The program defines the parameters (modulus of elasticity (Ey), 
Poisson ratio (vy) for the orthotropic linear elastic material model 
[Ref. 6:pp. 5.5-6). The material defined is assigned an identifying 
number, which can be referenced by several or all of the element 
groups (i.e., parent structure elements, repair patch elements). The 
material directions correspond to the global coordinates of the finite 
element model [Ref. 6:pp. 5.5-7). 

In ADINA. the composite material model can be used by several 
element types: 2-D solid element, 3-D solid element, plate element, 
and shell element [Ref. 6:p. 5.5-3). In this thesis, the 2-D element and 
3-D element t}-pe were used. 

C. VERIFICATION 

Part of the establishment of the repair guideline is the verification 
of the finite element code. This is done by modeling a problem for 
which the results are known and comparing numerical with these 
known results. 

The verification includes the geometry verification, material axes 
definition, and their transformation to the global coordinates, and 
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compares the 2-D (plane stress) solution and the 3-D solution with 
theoretical solution done by hand. 

A simply supported anisotropic square plate, under uniform ten- 
sile load along one edge, is considered as a verification problem and 
discussed in Appendix B. 
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APPENDIX B 



VERIFICATION PROBLEM 

A. PHYSICAL PROBLEM 

A simply supported anisotropic square plate, under uniform ten- 
sile load along one edge, is considered. The plate is built of eight plies 
and is symmetrical, laminated once with stacking sequence [902/752ls 
and once with stacking sequence [904] s . 

B. OBJECTIVE 

To verify the performance of the 2-D solid plane stress element, 
and then the performance of the 3-D solid element, when the material 
is anisotropic. 

C. FINITE ELEMENT MODEL 

For 2-D analysis, eight-plv thickness plate (laminate) is consid- 
ered. The mesh consists of three 2-D solid elements (nine nodes). An 
orthotropic linear elastic material model is used with the principal 
material axes a and b, making angle beta with the global coordinate 
axes X and V, respectively. The laminate stiffness is calculated by the 
MIC-MAC program [Ref. 7:ch. 7] from the specific sequence and 
orientation of the laminae in the laminate. 

For 3-D analysis, each ply (lamina) is considered separately across 
the thickness. The mesh consists of 12 3-D solid elements (27 nodes). 
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The same orthotropic linear elastic material as in the 2-D model is 
used. The lamina stiffness is used directly in the finite element model. 

For any model, through ADINA execution, a printout file is cre- 
ated. The printout file generates a listing of the displacements and the 
stress/strain for each element. Only the stress/strain results of nodes 
number 9, for the 2-D problem, or node number 21, for the 3-D 
problem, are printed. Node 9 and node 21 are located at the element 
geometric centroid and represent the element mean values for calcu- 
lating the reliability. 

D. SOLUTION RESULTS 

The solution results are as expected. For the lay-up of [904] s , the 
stresses in the applied loads direction are the same as the applied 
loads. No interlaminar stresses are shown (tend to zero). The results 
of 2-D analysis are in good agreement with the results of 3-D analysis. 

For the lay-up of [902, 752)s in 3-D analysis, the layers in 90° ori- 
entation are earning more stresses than the layers in 75° orientation, 
and the interlaminar stresses are shown (they cannot be neglected). In 
2-D analysis, because of plane stress assumption, the interlaminar 
stresses are not shown, and the stresses have the same value as 
before— they are representing the mean stress. 
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AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 



PROGRAM TRANS 2 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

INTEGER NLA , NEL 

PARAMETER (NLA=6,PI = 3.1415927 ) 

PARAMETER ( ROWS=100 ,C0LS=6 ) 

REALMS A(ROWS,COLS) , B ( ROWS , COLS ) 

REAL MOl , M02 ,NU12 , NU21 ,EX , EY ,EXY 
REAL THETA ( NLA) ,M,N 

REAL Q11,Q12,022,066,G12,E1 (NLA) ,E2 (NLA) ,E6(NLA) 
CHARACTERS 40 , FNOUT 

A 

WRITE ( 5 , A ) ' OUTPUT FILE , NAME : ' 

READ (5, ' (A) ' ) FNOUT 

WRITE ( 5 , s ) 'INPUT MODULUS & POISSON MOl ,M02 ,NU21 ,G12 ?' 
READ ( 5 , a ) MOl ,M02 ,NU21 ,G12 
NU12 = M02SNU21 / MOl 
Qll = MOl / (1-NU12SNU21) 

Q12 = Q11SNU12 

Q22 = M02 / (1-NU12SNU21) 

Q66 = G12 

A 

WRITE ( 5 , s ) 'INPUT STRAIN EX,EY,EXY ?' 

READ (5, A) EX,EY,EXY 

OPEN (UNIT=13 ,FILE=FNOUT , STATUS = 'NEW' ) 

WRITE (13, s) 'INPUT : ' 

WRITE (5,10) EX , EY , EXY 
WRITE (13,10) EX , EY , EXY 

A 

WRITE (5,A) 'INPUT THE ORIENT. OF EA. ELEM. IN DEG. ?' 
WRITE (13, A) 'INPUT: THE ORIENT. OF EA. ELEM. IN DEG. ; 
DO 50 1=1, NLA 

WRITE (5,20) I 
READ ( 5 , a ) THETA ( I ) 

WRITE (13,25) I ,THETA( I ) 

THETA ( I ) = THETA(I)API / 180. 0 

50 CONTINUE 

A 

WRITE (5, A) 'OUTPUT : ' 

WRITE (13, A) 'OUTPUT : ', FNOUT 

WRITE (5,60) ' NLA' , ' SI ' , ' S2 ' ,'S3' , ' S6 ' , ' S5 ' ,'S4' 

WRITE (13,60) 'NLA' , ' SI ' , ' S2 ' , ' S3 ' , ' S6 ' , ' S5 ' , ' S4 ' 

A 

NEL = 0 

DO 90 J=1 , 3 
DO 80 1=1, NLA 
NEL=NEL+1 

a STRAINS TRANSFORMATION : 

M=C0S ( THETA ( I ) ) 

N=SIN (THETA( I ) ) 

El (I) = EXAMAA2 + E Y a N a a 2 + EXYAMAN 
E2 ( I ) = EXANAA2 + EYAMAA2 + EXYAMAN 
E6 ( I ) = -EXA2AMAN + EYA2AMAN + EXYa(Maa2 - NAA2) 
a STRESSES CALCULATION : 



A ( 


NEL , 


.1) 


= OH A E1 ( 


I) 


+ Q12AE2(I) 


A ( 


NEL , 


2 ) 


= Q12AEK 


I) 


-f Q22AE2(I) 


A( 


NEL , 


, 6 ) 


= 066AE6( 


I) 




A( 


NEL, 


r 4 ) 


= 0. 






A ( 


NEL , 


.5) 


= 0. 






A ( 


NEL , 


r 3 ) 


= 0. 
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80 CONTINUE 

90 CONTINUE 

DO 100 1=1 , NEL 

WHITE (13,70) NEL ,A(I,1) , A ( 1 , 2 ) ,A(I,3) ,A(I,6) , A ( I , 5 ) , A ( I , 4 ) 
100 CONTINUE 

A 

CALL BETA ( A, B, NEL) 

WRITE (5, A) ' INPUT ALFA' 

READ ( 5 , a ) ALF 
WRITE (13,40) ALF 
CALL ALFA ( A , B , NEL , ALF ) 

CALL RELI( A,B,NEL) 

CALL ELRE( B ,NEL ) 

A 

10 FORMAT ( // ,1X, # STRAINX= ' ,E11 . 3 , SX , ' STRAINY= ' ,E11.3, 

& 9X , ' STRAINXY= ',E11.3) 

20 FORMAT ( IX , ' THETA (',12,')= ') 

25 FORMAT ( IX, ' THETA ( ' , 12 , ' ) = ',F6.1) 

40 FORMAT (//, IX SHAP FUNCTION (ALFA) = ' ,F5 . 1 / IX, 40 ( ' - ' ) ) 

60 FORMAT ( / / , IX , A10 , 6 ( A7 , 4X) / 1X,80('-')) 

70 FORMAT ( / , 4X , 16 , 6E1 1 . 3 ) 

A 

STOP 

END 

AAAAAAAAAAAAAAAAAAAAA+ A A A A A A A A 

SUBROUTINE BETA( A,B ,NEL) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( ROWS=100 ,C0LS=6 ) 

REALMS A( ROWS, COLS ) ,B( ROWS, COLS) 

REAL X(COLS) ,X1C,X2C,X3C 

A 

* PLY STRESS - STRENGTH, RAT 10 CALCULATION : 

A 

WRITE (13,10) 'NEL' , ' SRI ' , ' SR2 ' , ' SR3 ' , ' SR6 ' , ' SR5 ' , ' SR4 ' 

DO 90 1=1, NEL 

X(l) = 15.0E4 
X( 2 ) = 0.4E4 

X( 3 ) = X ( 2 ) 

X ( 6 ) = 0.68E4 

X( 4 ) = X ( 6 ) 

X( 5 ) = X( 6 ) 

X1C = -15.0E4 
X2C = -2.46E4 
X3C = X2C 

IF ( A( I , 1 ) . LT . 0 ) X(l) = X1C 
IF ( A ( I , 2 ) . LT . 0 ) X( 2 ) = X2C 
IF (A( 1,3) .LT. 0) X ( 3 ) = X3C 
DO 80 J= 1 , COLS 

B ( I , J ) = ABS ( A( I , J ) / X(J) ) 

60 CONTINUE 

WRITE (13,60) I ,B ( I , 1 ) ,B(I,2) ,B(I,3) ,B(I,6) ,B(I,5) ,B(I,4) 

90 CONTINUE 

10 FORMAT ( // ,1X,A1Q,6(A7,4X) / 1X,80('-')) 

60 FORMAT ( / , 4X, 16 , 6E11 . 3 ) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE ALFA( A , B , NEL , ALF ) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
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PARAMETER (ROWS=100,COLS=6) 

REAL* 6 A (ROWS, COLS) , E (ROWS , COLS ) 

A 

a PLY STRESS - STRENGTH, RATIO RISED BY ALFA CALCULATION : 

A 

WRITE (13,10) 'NEL' , 'VI' , 'V2' , 'V3' , 'V6' , 'V5' , 'V4' 

DO 90 1=1, NEL 
DO 80 J=l,COLS 

A ( I , J ) = B ( I , J ) **ALF 
BO CONTINUE 

WRITE (13,60)I,A(I,1) ,A(I,2) , A ( I , 3 ) ,A(I,6) , A ( I , 5 ) , A ( I , 4 ) 

90 CONTINUE 

10 FORMAT (//,1X,A10,6(A7,4X) / 1X,80('-')) 

60 FORMAT ( / ,4X, I6,6E11. 3) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE RELI ( A,B,NEL) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER (P.OWS=10O,COLS=6) 

REAL* 8 A ( ROWS , COLS ) , B ( ROWS , COLS ) 

A 

a PLY STERSS TENSOR RELIABILITY CALCULATION (WEIBUL) : 

A 

WRITE (13,10) 'NEL','R1\'R2','R3','R6','R5','R4' 

DO 90 1=1, NEL 
DO BO J= 1 , COLS 

E ( I , J ) = EXP(-A(I,J) ) 

BO CONTINUE 

WRITE (13,60) I , B ( I , 1 ) ,B(I,2) ,B(I,3),E(I,6) ,B(I,5) ,E(I,4) 

90 CONTINUE 

10 FORMAT (//,1X,A10,6(A7,4X) / 1X,80('-')) 

60 FORMAT ( / , 4X, 16 ,6E11 . 3 ) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAj-AAAAAAAAAA 

SUBROUTINE ELRE( B ,NEL) 

AAAAAAAAAAAAAAA^AAAA^-^-AAA^AA*' 

PARAMETER ( ROWS =100 , COLS =6 ) 

REALMS B (ROWS , COLS ) , ELR ( ROWS ) ,STR 

A 

* ELEMENTS RELIABILITY CALCUL. 6 STRUCTURE RELIABILITY CALCUL. 

A 

WRITE (13,10) ' NEL ' , ' ELM. RE' 

SIR = 1.0 
DO 90 1=1, NEL 
ELR ( I ) =1.0 
DO 80 J=1 ,COLS 

ELR ( I ) = ELR ( I ) * B( I , J ) 

80 CONTINUE 

WRITE (5,60) I ,ELP ( I ) 

WRITE (13,60) I , ELR ( I ) 

STR = STR a ELR ( I ) 

90 CONTINUE 

WRITE (5,70) STR 
WRITE (13,70) STR 

10 FORMAT ( / / , IX , 2A10 / 1X,25('-')) 

60 FORMAT ( / , 4X, 16 , Ell . 3 ) 

70 FORMAT (//, IX ,' STRUCTURE RELIABILITY = ',E11.3 / IX , 40 ( ' = ' ) ) 



4 5 



RETURN 

end 



APPENDIX C 



POST PROCESSOR PROGRAMS 

We have the state of stress in the spatial locations defined from 
the FEM (ADINA code). The next steps for evaluating the reliability- 
are: 

1. To normalize the size of the element. 

2. To perform, per element, the transformation of stresses, read 
from the FEM output file, to the lamina (composite material) 
principal directions. 

3. To compute the reliabilty based on a specific reliability model. 
This process is done by two programs— CLAREA and TRANS3 or 

TRANS 2. 

CLAREA is the program which reads in the coordinates as well as 
the element area defining the nodes. It then computes the element 
area and the smallest element area. 

TRANS3 is the program used for 3-D analysis. It reads in the 
stress components of each element as well as the element area and 
the angle theta, coresponding to the angle between the global coordi- 
nate system and the layer pricipal coordinate system. It then performs 
the transformation of stress to the principal direction. In addition, 
this program calls the subroutine BETA, which in turn calls subroutine 
ALFA, which in turn calls subroutine RELI, and which finally in turn 
calls subroutine ELRE. 
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TRANS2 is the program used for 2-D analysis. It reads in the 
strain components as well as the modulus of elasticity utilized in the 
repair or the parent structure, the angle theta, and the element area. 
It then calculates the stress tensor using the constitutive relation. 
From this point onward, the program is the same as the TRANS3 
program. 

These programs are illustrated by solving the verification exam- 
ple, as shown below. 
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2-D PROBLEM 

(Plane Stress 

[ 90 4 ] s ) 
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FILE 



a ADINA-IN INPUT 

A 

A THV2D. IN ANALYSIS OF UNISOTROPIC PLATE 

A 



FILEUNITS LIST=8 LOG=7 ECHO=7 
CONTROL PLOTUNIT=CM 
COLORS BCODE=BLUE 

A 

DATABASE CREATE 
WORKSTATION DEVICE=0 

A 

HEAD 'THV2D.IN ANALYSIS OF UNISOTROPIC PLATE ( 90,90 DEG )' 
A 

MASTER IDOF =000111 

PRINTOUT VOLUME=MAXIMUM P=STRAINS 

PORTHOLE FORMA TTED= YES FILE=60 

A 

COORDINATES 
ENTRIES NODE Y 



1 

2 

3 

4 

5 

6 

7 

8 
9 

A 

A 

FIXE 123 / 
FIXB 13 / 

A 

MATERIAL 1 



20 . 
- 20 . 
20 . 
- 20 . 
20 . 
- 20 . 
20 . 
- 20 . 
0. 



9 

7 8 



Z 

20 . 

20 . 

5. 

5. 

- 5 . 

- 5 . 

- 20 . 

- 20 . 

- 20 . 



ORTHOTROPIC EA=181.E05 EB=10.3E05 EC=10.3E05, 

NUAB = 0 .0159 NUAC = 0 .0159 NUBC = 0.28 GAB = 7.17E05 



EGROUP 1 TWODSOLID STRESS2 TABLES 



3 EL1=1 EL2=1 NODES =9 
7 EL1= 1 EL2= 1 NODES =9 



BET TABLE 
90 1 

90 1 



GSURFACE 124 
GSURFACE 568 
STRESSTAELE 1 9 
EDATA 

ENTRIES" EL THICK 
1 0.1 
2 0.1 
A 

LOADS ELEMENT 
1 2 -3.0E04 

A 

EGROUP 2 TWODSOLID STRESS2 TABLES 
GSURFACE 3465 EL1=1 EL2=1 NODES=9 
STRESSTABLE 1 9 
ED7^TA 

ENTRIES EL THICK 
1 0.1 
A 
A 

VIEW ID= 1 XV I EW = 1 YV I EW = 0 ZVIEW = 0 
EGZONE NAME=TW0D30LID 
1/2 



BET TABLE 
90 1 
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FRAME HEADING=UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH ZONE=TWODSOLID N0DES=11 ELEMENT =2 BCODE=ALL VIEW=1 

A 

ADINA 

A 

END 
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PROGRAM ADINA VERSION ADINA 5.0/NL2 THV2D. IN ANALYSIS OF UNISOTROPIC PLATE ( 90,90 DEG 
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2-D PROBLEM 

(Plane Stress 
[90 2 /75 2 ] s ) 



*Also includes output of MIC-MAC program 
showing the laminate stiffness Ej 0 (page 2. 
column w, rows 22-28) 



53 





A 


B 1 


c 


D 


E 1 


~r~T 


G 


H 1 


1 


1 


MIC-MAC/IH-PlArC: {[theta/'*], . . .}S j 


IPly mat: 1 


T3/N52[SI] | 




2 


(theta) j 


-82.5; 


82.5 = 


82.5! 


-82.5! 


[repeat] i 


* 1 


M-3 J 


[Rotate] 


3 


(•/group) ! 


10! 


1 .0! 


i .of 


1.0! 


1.0 1 


8.0! 


10! 


0.00 


4 


1 • 1 • 


ill! 




5 


R /intact i 


6.08! 


6.08: 


6.08! 


6.08! 


R/FPF 1 


6.08: 


safety 1 


1.50 


6 


R /degraded! 


627| 


627j 


627: 


627j 


R/LPF | 


6.27j 


R/utt | 


627 


7 


[ : i l 


t 


R/lim | 


4.18! 


R/lim* 1 


4.18 


8 


{N},MN/mor k/x. | 


<sigrna> j 


<sig>1im | 


<sig>lim*| <sig>ult 1 


{E°)lim j 


E‘u/E‘1 




9 


1 


oooi 


o.i 


O.i 


O.I 


O.I 


1 0.5‘[ 


0.436! 




10 


2 | 


0.25; 


250. : 


1045.! 


1045.1 


1567.! 


173.4! 


0.997: 




11 


6 j 


o.ooj 


0.! 


0.1 


o.i 


o.i 


9.8! 


0.579! 




12 




| S | j j 


13 




<N} lim ! 


<eps>E-3 j 


<eps>1im!<eps>lim*; <eps>u1t !a)pho,E-6 j 


beta 0 




14 


1 


0.001 


-0.75- 


-3.14! 


! -3.14! 


-4.71 1 


! _21794j _ 


0.581! 




15 


2 j 


1 .04 j 


1.44 


6.03! 


603! 


904! 


r -0278!* 


-0.008 




16 


6 


0.00: 


o.oo’ 


0.00 


0.00| 


0.00; 


0 . 000 ! 


0.000; 




17 






| | | j! 


18 






c,moist 


YOl/f 


1 Em j 


Efx 


| Xrn 1 


Xfx 


iEm/Em 0 


19 


Baseline j 


22.0! 


0.005! 


! 070 


j 3.40! 


259 


40.0! 


2143 


I 0.30 


20 


[Modified] j 


22. of 


" 0.005; 


070 


3.40! 


259 


40.0! 


2143 


1 0.30 


21 


Mod/Bsin j 


\ .000! 


1.000! 


1 1.000 


L 1.0Q0; 


1 OOP 


! 1000| 


7000 


; 1.000 


22 


Mod-Bsln | 


” 6. of 


0.000! 


! Hot/Vet 


" 340; 


” 259 ’ 


r " 400 !* 


2143 




23 


/Q-iso | 


at limit j 


at ult ! 






j j 


1 


24 


stiffness j 


2.30j 


2.30 






: : 


1 


25 


strength i 


•REF if 


’ •RE?! 








I 


26 


j 1 | 




j j j 




27 


Ply data linked to Ply Data File: i 






j 


28 


T3/N52[SI| 


18lj 


10.3 


j 0.28 


7.17 


3.4 


I 122! 


160 


i 


29 


1500: 


15661 


40 


246 


1 68 


1.2 


I 22; 


0.005 




30 


-6.5! 


" 1 25 j 


07 


1.6 


| 0.516! 


0.5 


! 0.2! 


0.9 




31 


0.02! 


22.5! 


0 


0.6 


1 0.316 


j 0.004 


j 0.004; 


2000 


j 


32 


69.7! 


! 325; 


0.916 


1.56 


0 


| 0 


i o| 


0 


: 


33 














i i i 


34 


i • : • i 






35 




; : ; : ! 




! 1 i 


36 




! • J • 


[ 


: j | 


37 




| ! j 1 




t : 


38 










39 




i : j : 


i 


i j j 


40 




i ; : 


i • : t 
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K 




s 1 


T 


u 


V 1 


V X 


Y 


1 


INTACT LAMINATE MODULUS MODULE - elastic and hggrothermal constants 1 




2 




\ \ \ \ \ \ 




3 


{Angle] i 


theta/1 1 


theta 12 ! 


theta /3 j 


theta /4 j 






4 


[theta] I 


-82 .5 1 


82.5j 


82.5; 


-82.51 


[REPTJ+ ! 




5 


l»/grp) j 


i.oi 


1.0; 


1.0! 


1.0) 


1 1 I 




6 


2X/ad | 


-3E+00| 


3E+Oo| 


3E+00; 


-3E+00) 


h+/r+,* j | 




7 


4X,rad 


-6£+00j 


6E+00! 
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-6E+00! 


4 j i 




8 




i [ : : 


: : 




9 


Top z* 1 
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: ) 




10 
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o5o: 


! 




11 
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t 


12 






IE-03) ] 




13 






j i i 




14 
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(01/t | 


(01/2 j 


(0]/3 ] 


101/4 ! 


! [A] ) [A*] j 




15 


11 : 
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10.631 


10.63! 


10.63: 


0.6i j 8E-03 




16 


22 
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176251 


1 76.25 j 


17625] 
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17 
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5.54 1 
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0.01 1 4E~03 
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66 __ ' 
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19 
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20 


62=26 | 
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21 








lAl 
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22 


Compl 


i la]/n/GN j 


U*E 




j Ei 0 < i 




23 


ii 
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IE-01 j 
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24 


22 
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25 
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26 
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27 
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Baseline j 


22.0: 


0.005: 


0.70 


3.40 


259 


40.0 


2143 


| 0.30 


20 


[Modified] I 


22.0 


! 0.005! 


0.70 


3.40 


259 


40.0 i 


2143 


j 0.30 


21 


Mod/Bsln | 


1.000 


1.000! 


1.000 


; 1.000 


I 1.000 


1 .0001 


1.000 


ij 1.000 


22 


Mod-Bsln j 


O 

d 


0.000! Hot/Vet : 


3.40 


259 


40.0! 


2143! 


23 


/Q-iso 


at limit ; 


| at ult 1 




| J 


24 


stiffness j 


1.36; 


0.96! 


||| 


j : 


25 


strength ! 


*REF! 


•REF!! 




: : 


26 




i ♦ 


\ \ \ 


i 1 


27 


Ply data linked to Pig Data File : j 


j : : 1 • 


28 


T3/N52[SlJ 


181 


! 10.3: 


0.28 


1 7,17 


! 3.4 


122 


160! 


29 


1500! 


’ 1500 


40[ 


246' 


68 


12 


22 


0.005! 


30 


-0.5! 


125 


I 0.7! 


1.6 


! 0.516 


0.5 


02 


0.9; 


31 


0.02! 


22.5 


0! 


0.6 


| 0.316 


| 0.004 


0.004 


2000! 


32 


69.7; 


325 


! 0.916! 


1.56 


0 


! 0 


0 


0 


ij 


33 






34 




: ; r 1 j [ j 


35 




: : j I : j | 


36 
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| j : j : j | 


38 


j 


\ \ \ \ \ 1 : 


39 


j 


: : 




j 


| 


; 


\ 


j 


40 


[:!;:[ || 
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R 


S 


T 


U 


TT 


V 


x 1 


Y 


1 


INTACT LAMINATE MODULUS MODULE - elastic and hygrothermal constants 




2 


i | 


j : j • s 




3 


l Angle! 


1 theta /I 1 


theta/2 1 


theta/3 j 


theta/4 ! 


j 


i 




4 


[theta] 


| 90.0! 


90.0! 


75.0! 


75.0! 


[REPT ]+ 1 


: 




5 


[*/grp] 


1.01 


1 .01 


1.0! 


1 .of 


1 [ 






6 


2X/ad 


1 3E+00; 


3E+00! 


3E+00: 


3E+Oo| 


llt/r +*A 


j 




7 


4X,rad 


1 6E+00j 


6E+00! 


5E+00! 


5E+00! 


4 j 




8 




} j : j 


r 


9 


Top z* 


! 1 .00! 


0.75! 


0.50! 


0.25! 




j 


10 


Bott z* 


0.75; 


0.50! 


025! 


o.ooj 




j 




11 


de1(z*) 


! 0.25j 


! i 
1 
1 

| i 
io i 
Ko | 


025| 


0.25! 


h+ | 






12 




: 




1 


IE-03! 






13 




1 


j | : 






14 


Stiff 


! loi/i i 


[01/2 | 


[01/3 j 


IQ1/4 | 


[A] I 


[A*] | 




15 


11 


! i 0.35 : 


10.35! 


1 1 .98 1 


1 1 .98! 


0.01! 


9E-03! 




16 


22 


i 181.81* 


181.81! 


160.47! 


160.47!" 


0.171 


IE-01! 


\ 


17 


21 = 12 
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2.90! 


12.75! 


1 2.75 •: 
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7E-03 


1 


18 


66 


7.17 
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0.01 1 
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19 


61 = 16 


J OJOOj 


0.00! 
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o.ool 
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20 


62=26 


j" 0.00 i 
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38.50j 


38.50! 
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2E-02 




21 


i : : 


: 


IAI 


2E-05! 






22 


Compl 


j [al/n/GN 




(ill. L 


: 


Eio | 






23 


ii 


~1 94.04j 




9E-02! 


i 


10.63 




• ' 

| 


24 


22 


T" 7.21 
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! 


138.75 






25 


21 = 12 


| -2.91 


■ i 
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j 102.30; 
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| j E,GPa 
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Nonmechanical stress(Pa) and stra 


in" 1 1 V*/iA 
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! sig‘n fl j 
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INPUT 



FILE 



a ADINA-IN 

A 

a THV2D. IN ANALYSIS OF UNISOTROPIC PLATE 

A 



A 

FILEUNITS LIST=8 LOG=7 ECHO =7 
CONTROL PLOTUNIT=CM 
COLORS BCODE=BLUE 

A 

DATABASE CREATE 
WORKSTATION DEVICE=0 
A 

HEAD 'THV2D.IN ANALYSIS OF UNISOTROPIC PLATE ( 90,75 DEC )' 

A 



MASTER IDOF=OOOH1 

PRINTOUT VOLUME = MAXI MUM P=STRAINS 

PORTHOLE FORMATTED=YES FILE=60 



A 

COORDINATES 
ENTRIES NODE Y 
1 20 . 
2 - 20 . 

3 20. 

4 -20. 

5 20. 

6 - 20 . 

7 20. 

8 - 20 . 

9 0. 

A 

A 

FIXB 123 / 9 
FIXB 13 / 7 8 

A 



20 . 

20 . 

5 . 

5. 

- 5 . 

-5. 

- 20 . 

- 20 .' 

- 20 . 



MATERIAL 1 ORTHOTROPIC EA=173 . 36E05 EB=10.45E05 EC=10.3E05, 

NUAB=0 . 03 NUAC=0 . 03 NUBC=0.28 GAB=9.81E05 

A 



EGROUP 1 TWODSOLID STRESS2 TABLES 
GSURFACE 1243 EL1=1 EL2=1 NODES=9 
GSURFACE 5687 EL1=1 EL2=1 NODES=9 
STRESSTABLE 1 9 
EDATA 



ENTRIES EL THICK 
1 0.1 
2 0.1 



BET TABLE 
82.5 1 

82.5 1 



A 



LOADS ELEMENT 
1 2 -3.0E04 

A 



EGROUP 2 TWODSOLID STRESS2 TABLES 
GSURFACE 3465 EL1=1 EL2=1 NODES=9 
STRESSTABLE 1 9 
EDATA 

ENTRIES EL THICK BET TABLE 

1 0.1 82.5 1 

A 



A 

VIEW ID=1’XVIEW=1 YVIEW=0 ZVIEW=0 
EGZONE NAME=TWODSOLID 
1/2 
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FRAME HEADING=UPPER XSF=1.0 Y5F=0.5 XFMAX=22.5 YFMAX=17.0 
MESH 20NE=TW0DS0LID N0DES=11 ELEMENT=2 BCODE=ALL VIEW=1 

A 

ADINA 

A 

END 
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PROGRAM ADINA - VERSION ADINA 5.0/NL2 THV2D.IN ANALYSIS OF UNISOTROPIC PLATE ( 90,75 DEG 
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STRESSES ARE MEASURED IN GLOBAL COORDINATE SYSTEM 
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3-D PROBLEM 

[90 4 ] s 
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INPUT 



FILE 



A ADINA-IN 

A 

a THV3D.IN EIGHT SYMMETRIC LAYER USING 3-D. 

A 

FILEUNITS LIST= 8 L0G=7 ECH0=7 
CONTROL PLOTUNIT=CM 
COLORS BCODE=BLUE 

A 

DATABASE CREATE 
WORKSTATION DEVICE=0 

A 

HEADING ' THV3D.IN; aaaeIGHT LAYERS SOLUTION 90,90 DEG aaa 1 

A 

MASTER IDOF=OOOH1 

PRINTOUT VOLUME = MAX I MUM IPD = 1 IPRIC=0 P=STRAINS 
PORTHOLE FORMATTED=YES FILE=60 

A 

COORDINATES / ENTRIES NODE X Y Z 



1 


-20. 


-20. 


0 . 1 


2 


20. 


-20. 


0.1 


3 


-20. 


-5. 


0.1 


4 


20. 


-5. 


0.1 


5 


-20. 


-20. 


0 . 


6 


20. 


-20. 


0 . 


7 


-20. 


-5. 


0 . 


8 


20. 


-5. 


0 . 


9 


-20. * 


5. 


0.1 


10 


20. 


5. 


0.1 


11 


-20. 


20. 


0.1 


12 


20. 


20. 


0.1 


13 


-20. 


5. 


0 . 


14 


20. 


5. 


0 . 


15 


-20. 


20. 


0 . 


16 


20. 


20. 


0 . 


17 


0 . 


20. 


0.1 


18 


0 . 


20. 


0 . 


101 


5. 


0 . 


0 . 


102 


0 . 


10. 


0 . 


103 


5. 


10. 


0 . 


104 


15. 


2.68 


0 . 


105 


7.68 


10. 


0 . 


106 


15. 


10 . 


0 . 


107 


15. 


-10. 


0 . 


; NOTE THE MAXES MEANS : 1=90DEG 


2=15 3= 


75 4=45 5=- 


45 



MATERIAL 1 ORTHOTROPIC 1B1.E05 I0.3E05 I0.3E05 0.0159 

GAB = 7 . 17E0 5 7.17E05 7.17E05 

A 

EGROUP 1 THREEDSOLID MATERIALS TABLES 

GVOLUME 12435687 EL1=1 EL2 = 1 EL3 = 4 NODES=27 
GVOLUME 9 10 12 11 13 14 16 15 EL1=1 EL2=1 EL3 = 4 N0DES=27 
AGVOLUME 12435687 EL1=1 EL2 = 1 EL3=4 N0DES=8 
A GVOLUME 9 10 12 11 13 14 16 15 EL1 = 1 EL2 = 1 EL3=4 NODES=B 
GSURFACE 11 12 16 15 EL1 = 1 EL2 = 4 N0 = 9 
LINE STRAIGHT 17 18 EL=4 M=1 NC=ALL 

1 NODI = 101 NODJ = 103 NODK=102 

2 NODI= 10 1 N0DJ=1 04 NODK=102 

3 N0DI = 101 NODJ = 105 NODK=102 

4 NODI = l 01 NODJ = 106 N0DK=102 



0.0159 0.28, 



AXES-ORTHOTROPIC 
AXES -ORTHOTROPIC 
AXES-ORTHOTROPIC 
AXES-ORTHOTROPIC 
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AXES-ORTHOTROPIC 5 NODI=101 NODJ=107 NODK=102 
*STRESSTABLE 1 12 3 4 5 6 7 B 21 26 27 
STRESSTABLE 1 21 



EDATA 


/ ENTRIES 


EL 


TABLE MAXES 








1 


I I 








2 


1 1 








3 


1 1 








4 


1 1 








5 


1 1 








6 


1 1 








7 


1 I 








8 


1 1 


A 










LOADS 


ELEMENT 


1 


2 


- 3 . 0E04 




TO 


4 


2 


-3.0E04 



A 

FIXB 123 LINES / 17 18 

FIXB 23 SURFACES / 11 12 16 15 

AFIXB 23 LINES / 12 16 / 16 11 / 11 15 / 15 16 . 

A 

A 

EGROUP 2 THREEDSOLID MATERIALS TABLES 

GVOLUME 3 4 10 9 7 8 14 13 EL1=1 EL2=1 EL3=4 NODES=27 

*GV0LUME 3 4 10 9 7 8 14 13 EL1=1 EL2=1 EL3=4 N0DES=8 

AXES-ORTHOTROPIC 1 NODI=101 N0DJ=103 NODK=102 
AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 
AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 
AXES-ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 
AXES-ORTHOTROPIC 5 NODI=101 NODJ=107 NODK=102 
^STRESSTABLE 112345678 21 26 27 
STRESSTABLE 1 21 

EDATA / ENTRIES EL TABLE MAXES 
1 1 1 

2 11 

3 1 1 

4 1 1 

* 

A 

EGZONE NAME=THREEDSOLID / 1 / 2 

FRAME HEADING =UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH ZONE=THREEDSOLID NODES = 01 ELEMENT =1 H I DDEN = DASHED BCODE=ALL 

A 

ADINA 

A 

END 
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r N P U T 



FILE 



a A D I N A - IN 

A 

A THV3D.IN EIGHT SYMMETRIC LAYER USING 3-D. 

A 

FILEUNITS LIST=8 L0G = 7 ECH0=7 
CONTROL PLOTUNIT=CM 
COLORS BCODE=BLUE 

A 

DATABASE CREATE 
WORKSTATION DEVICE=0 

A 

HEADING ' THV3D.IN; aaaeIGHT LAYERS SOLUTION 

A 

MASTER IDOF =000111 

PRINTOUT VOLUME=MAXIMUM IPD=1 IPRIC=0 P=STRAINS 
PORTHOLE FORMATTED=YES FILE=60 

A 

COORDINATES / ENTRIES NODE X Y 

1 - 20 . - 20 . 

2 20 . - 20 . 

3 -20. -5. 

4 20. -5. 

5 -20. -20. 

6 20 . - 20 . 

7 -20. -5. 

8 20. -5. 

9 - 20 . 5 . 

10 20. 5. 

11 - 20 . 20 . 

12 20 . 20 . 

13 -20. 5. 

14 20. 5. 

15 -20. 20. 

16 20. 20. 

17 0. 20. 

18 0. 20. 

101 5. 0. 

102 0 . 10 . 

103 5. 10. 

104 15. 2.68 

105 7.68 10. 

106 15. 10. 

107 15. -10. 



90,75 DEG aaa' 



Z 

0 

0 

0 

0 

0 . 

0. 

0 . 

0. 

0. 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 

0 . 

0 . 

0 . 



1 



ANOTE THE MAXES MEANS : 1=90DEG 2=15 3=75 4=45 5=-45 

A=====================================================: 

MATERIAL 1 ORTHOTROPIC 181. E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 

GAB =7 . 1 7E05 7.17E05 7.17E05 

A » 

EGROUP 1 THREEDSOLID MATERIALS TABLES 

GVOLUME 12435687 EL1=1 EL2=1 EL3=4 N0DES=27 
GVOLUME 9 10 12 11 13 14 16 15 EL1=1 EL2=1 EL3=4 N0DES=27 
’‘'GVOLUME 12435687 EL1 = 1 EL2 = 1 EL3 = 4 N0DES=8 
’“'GVOLUME 9 10 12 11 13 14 16 15 EL1 = 1 EL2 = 1 EL3 = 4 N0DES = 8 
GSURFACE 11 12 16 15 EL1=1 EL2=4 NO=9 
LINE STRAIGHT 17 18 EL=4 M=1 NC=ALL 
AXES-ORTHOTROPIC 1 N0DI=101 N0DJ=103 N0DK=102 

AXES-ORTHOTROPIC 2 NODI=101 N0DJ=104 NODK=102 

AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 N0DK=102 

AXES-ORTHOTROPIC 4 N0DI=101 NODJ=106 N0DK=102 
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AXES -ORTHOTROPIC 5 NODI =101 NODJ=107 NODK=102 
^STRESSTABLE 1 1 2 3 4 5 6 7 8 21 26 27 
STRESSTABLE 1 21 

EDATA / ENTRIES EL TABLE MAXES 







1 


1 






2 


1 






3 


1 






4 


1 






5 


1 






6 


1 






7 


1 






8 


1 


A 








LOADS ELEMENT 


1 


2 


-3.0E04 


TO 


4 


2 


-3.0E04 



A 

F1XB 123 LINES / 17 18 

FIXB 23 SURFACES / 11 12 16 15 

*FIXB 23 LINES / 12 16 / 16 11 / 11 15 / 15 16 

A 

A 

EGROUP 2 THREEDSOLID MATERIALS TABLES 

GVOLUME 3 4 10 9 7 8 14 13 EL1=1 EL2=1 EL3=4 NODES=27 

AGVOLUME 3 4 10 9 7 8 14 13 EL1=1 EL2=1 EL3=4 NODES=8 

AXES-ORTHOTROPIC 1 NODI=101 N0DJ=103 NODK=102 
AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 
AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 
AXES-ORTHOTROPIC 4 NODI=101 N0DJ=106 NODK=102 
AXES-ORTHOTROPIC 5 N0DI=101 N0DJ=107 NODK=102 
^STRESSTABLE 112345678 21 26 27 
STRESSTABLE 1 21 

EDATA / ENTRIES EL TABLE MAXES 
1 1 1 

2 13 

3 13 

4 11 
A 

A 

EGZONE NAME=THREEDSOLID / 1 / 2 

FRAME HEADING = UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH ZONE=THREEDSOLID NODES=01 ELEMENT=1 HIDDEN =DASHED BCODE=ALL 

A 

ADINA 

A 

END 
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FILE 



A ADINA-IN 

A 

a THV3D.IN EIGHT SYMMETRIC LAYER USING 3-D. 
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FILEUNITS LIST= 8 L0G = 7 ECH0=7 
CONTROL PLOTUNIT=CM 
COLORS BCODE=BLUE 
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WORKSTATION DEVICE=0 
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END 
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n n n n n o n 



AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
PROGRAM CLAREA 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 



A 



A 



THIS PROGRAM WHICH READS THE NODAL POINT COORDINATE, AS WELL READS 
THE NODES DEFINING THE ELEMENT AREA, AND READS THE ANGLE OF ELEMENT 
ORIENTATION FROM FILES CREATED FROM THE MAIN ADINA OUTPUT FILE. 

IN ADDITION THIS PROGRAM CALCULATES THE MINIMUM AREA. 



BEFORE RUNING THE PROGRAM VERIFY THE # OF ELEMENT ( NEL ) 
( NNOD) . 

INTEGER IEL, NEL, NNOD, INOD 

PARAMETER ( NEL=12 , NNOD=196 ) 

INTEGER INODE( NEL , 4 ) , NODE( 4 ) 

REAL COOR ( NNOD, 20 ) , AREA1 , AREA2 , AREA ( NEL) 

REAL B1,B2,C1,C2,BB1,BB2,CC1,CC2 

REAL X ( 4 ) , Y ( 4 ) 

CHARACTER’ 1 ' 10 FNAME1 , FNAME2 , FNOUT 



THE #OF NODES 



WRITE (5, A) 'ENTER FILES NAMES FROM ADINA OUTPUT FILE?' 
WRITE (5, A) 'NODAL POINT DATA FN ?' 

READ ( 5, ' (A) ' ) FNAME1 

OPEN ( UNIT=16 ,FILE=FNAME1 , STATUS = ' OLD' ) 

WRITE (5, A) 'NODAL POINT DATA FN : ' ,FNAME1 

WRITE (3, A) 'NODAL POINT DATA FN : ' ,FNAME1 

WRITE (5, A) 'ELEMENTS NODES FN ?' 

READ (5, ' (A) ' ) FNAME2 

OPEN ( UNIT=15 ,FILE=FNAME2 , STATUS= ' OLD' ) 

WRITE (5, A) 'ELEMENTS NODES FN : ' ,FNAME2 

WRITE (3, A) 'ELEMENTS NODES FN : ' , FNAME2 

WRITE (5, A) 'ENTER THE NAME OF ELEMENTS AREAS FILE ?' 

READ (5, ' (A) ' ) FNOUT 

OPEN ( UN IT= 1 3 , FILE=FNOUT , STATUS = ' NEW ' ) 

WRITE (5, A) 'AREAS OF ELEMENTS FILE : ', FNOUT 

WRITE (3, A) 'AREAS OF ELEMENTS FILE : ', FNOUT 



A DIFINING X,Y COORDINATES FROM NODAL POINT DATA FILE 

A 

a WRITE (5,30) 'NNOD' , 'XCOOR' , 'YCOOR' 

a WRITE (3,30) 'NNOD' , 'XCOOR' , 'YCOOR' 

DO 50 1=1, NNOD 

READ (16, A) INOD, 112,113, 114,115, 116,117, 118,119, 1110,1111, 

& 1112,1113,1114 ,C00R( I ,1) , COOR (1,2) ,A1 7, 1118,1119 

a WRITE (5, 40)1, CO OR (1,1) , COOR (I, 2) 

A WRITE (3,40)I,COOR(I,1),COOR(I,2) 

50 CONTINUE 

A 

A DIFINING ELEMENT'S NODES FROM ELEM. NODES FN. 

A 

WRITE (5,10) ' NEL ' ,. ' NODE1 ' , ' N0DE2 ' , ' N0DE3 ' , ' N0DE4 ' 

WRITE (3,10) 'NEL', ' NODE1 ' , ' N0DE2 ' , ' NODE3 ' , ' NODE4 ' 

AMIN=2000 . 0 
DO 100 IEL=1 , NEL 

READ (15, A) I IEL ,I2,I3,I4,I5,I6,I7,I8,A9,A10,I11, INODE ( IEL , 1 ) , 

& INODE (IEL, 2) , INODE( IEL , 3 ) , INODE( IEL, 4 ) ,116,117,118,119 

WRITEt 5,20) IIEL , INODE( IEL , 1 ) , INODE( IEL , 2 ) , INODE( IEL , 3 ) , 

& INODE( IEL , 4 ) 

WRITE ( 3,20) IIEL , INODE( IEL,1) , INODE ( IEL , 2 ) ,INODE( IEL, 3) , 

£ INODE( IEL , 4 ) 

NODE(l) = INODE ( IEL, 1) 

NODE( 2 ) = INODEt IEL , 2 ) 
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NODE ( 3 ) = INODE ( IEL , 3 ) 

NODE ( 4 ) = INODE ( IEL , 4 ) 

k 

a DIFINING THE X,Y COORD. TO NODE(l) - NODE(4) 

A 



WRITE (5,30) ' NNOD' , ' XCOOR ' 
WRITE (3,30) 'NNOD' , 'XCOOR' 

DO 80 J = 1 , 4 
DO 70 1=1, NNOD 

IF ( I . EQ . NODE( J ) ) THEN 

X(J) = COOR (1,1) 

Y ( J ) = COOR (1,2) 
WRITE ( 5 , 40 ) I ,X( J ) , Y ( J ) 
WRITE ( 3 , 40 ) I ,X( J ) ,Y( J) 
GO TO 80 
END IF 

70 CONTINUE 

80 CONTINUE 

A 

A ELEMENTS AREA CALCULATION. 

A 

B1 = Y ( 3 ) - Y { 4 ) 

B2 = Y ( 4 ) - Y ( 2 ) 

Cl = X( 4 ) - X( 3 ) 

C2 = X( 3 ) - X ( 2 ) 

BB1 = Y ( 4 ) - Y ( 1 ) 

BB2 = Y ( 1 ) - Y ( 2 ) 

CC1 = X ( 1 ) - X ( 4 ) 

CC2 = X( 2 ) - X ( 1 ) 

AREA1 = ABS ( C2*B1 - Cl AB2 ) / 2 
AREA 2 = ABS(CC2ABB1 - CC1ABB2) 



' YCOOR' 
' YCOOR ' 



100 



2.0 



AREA (IEL) = AREA1 +AREA2 
AEL=AREA( IEL) 

AMIN = MIN ( AMIN , AEL ) 

WRITE (5,90) I IEL, AREA ( IEL) , AMIN 
WRITE (3,90) I IEL, AREA ( IEL) , AMIN 
CONTINUE 





DO 110 
WRITE 
WRITE 


1=1 ,NEL 

(13.95) 1, AREA ( I ) , AMIN 

(5.95) 1, AREA ( I ) ,AMIN 


110 

A 


CONTINUE 


10 


FORMAT 


( / 5 ( A 7 , 4 X ) / 7 0 ( ) ) 


20 


FORMAT 


( / / 5 ( 1 7 , 4X ) / ) 


30 


FORMAT 


( /3(A7,4X) / 4 0 ( ' - ' ) ) 


40 


FORMAT 


( /I7,4X, 2E11.3) 


90 


FORMAT 


( //IX, ' AREA( ',12,') = 


95 

A 


FORMAT 

STOP 

END 


(*/ /IX, 12 , 2F10 . 2) 



, 2F10 . 2 ) 
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AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA*AAAAAAAAAAAAAAAAAAAAAAAAAA 

PROGRAM TRANS 3 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

C THIS PROGRAM TRANFORMS THE STRESSES OUTPUT FROM ADINA , 

C EXPRESSED IN THE GLOBAL STRUCTURE COORDINATE SYSTEM, TO 
C THE PRINCIPAL MATERIAL COORDINATE. AND THEN, CALCULATES THE 
C RELIABILTY PER STRESS COMPONENT AND PER ELEMENT AND FINALLY 

C EVALUATES THE TOTAL RELIABILITY. ALL THAT BY CALLING THE 

C SUBROUTINE HEREIN. 

C NEL : # OF TOTAL ELEMENTS. 

C SX , SY , . . . , SYZ : THE STRESS IN GLOBAL COORDINATE SYSTEM. 

C THETA : THE LAYUP ANGLE PER ELEMENT. 

C S1,S2 S6 s THE STRESS IN THE MATERIAL PRINCIPAL COORDINATE. 

C R1,R2,...,R6 : THE RELIABILITY PER STRESS COMPONENT. 

C RT1 ,RT2 , . . .RT6 : THE TOTAL RELIABILTY OF STRUCTURE PER COPMPONENT. 

C RE : THE RELIABILITY PER ELEMENT 

C RT : THE TOTAL RELIABILITY. 

C ALF : THE SHAPE FUNCTION PARAMER FOR WEI BULL 

C BEFORE RUNNING THE PROGRAM VERIFY IN THE PROGRAM THE ANGLE THETA 

C INPUT AND THE # OF ELEMENT. AND ALSO VERIFY THE MATERIAL PROEPRTIES 

C GIVEN IN SUBROUTINE BETA. 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

INTEGER NEL 

PARAMETER (NEL* 12, PI = 3.1415927 ) 

PARAMETER ( ROWS* 100 , COLS = 6 , COLS 1 = 7 ) 

INTEGER I THETA ( NEL) 

REAL*B A (RONS, COLS) , B (ROWS , COLS ) , C ( ROWS , COLS1 ) 

REAL SX(NEL) , SY ( NEL) ,SZ(NEL) ,SXY(NEL) ,SXZ(NEL) ,SYZ(NEL) 

REAL L , M , THETA ( NEL ) , THETAF ( NEL ) , ALF 
CHARACTER* 2 0 FNAME , FNOUT , FNAME1 , FNAME2 

A 

WRITE (5,*) 'ENTER INPUT STRESS FN ?' 

READ (5, '(A)') FNAME 

WRITE (5,*) 'INPUT STRESS FN = FNAME 
OPEN ( UNIT* 15 , FILE=FNAME, STATUS* 'OLD' ) 

WRITE (5,*) 'ENTER INPUT ELEM. ANGLE FN?' 

READ ( 5 , ' ( A) ' ) FNAME 2 

WRITE (5,*) 'INPUT ELEM. ANGLE FN = ' ,FNAME2 
OPEN ( UNIT=16 ,FILE=FNAME2 , STATUS* ' OLD' ) 

WRITE (5,*) 'ENTER AREAS OF ELEMENTS FN ?' 

READ ( 5 , ' (A) ' ) FNAME1 

OPEN (UNIT* 17 ,FILE=FNAME1 , STATUS* ' OLD' ) 

WRITE (5,*) ' ENTER OUTPUT FN ?' 

READ (5, ' (A) ' ) FNOUT 

OPEN (UNIT* 13, FILE* FNOUT, STATUS* 'NEW' ) 

WRITE (13,*) 'INPUT STRESS FN = ', FNAME 
WRITE (13,*) 'INPUT ELEM. ANGLE FN = ' ,FNAME2 
WRITE (13,*) 'INPUT AREAS OF ELEMENTS FN = ' ,FNAME1 
WRITE (13,*) ' OUTPUT FN =', FNOUT 

A 

* INPUT OF ANGLE (THETA) FOR EA. ELEMENT ; 

A 

DO 20 1=1, NEL 

READ (16,*) N,I2,I3,I4,I5,I6,ITHETA(I) , IB ,A9 ,A10 , 111 , 

& 112,113,114,115,116,117,118,119 

IF (ITHETA(I) .EQ.l) THETA ( I ) = 90.0 
IF (ITHETAt I) .EQ.3) THETA ( I ) =75.0 
WRITE (5,25) N , THETA ( I ) 

THETAF ( I ) = THETA ( I ) * PI/180.0 
20 CONTINUE 
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WRITE (13, A) 'OUTPUT :',FNOUT 

A 

a TRANSFORMATION OF STRESSES TO PRINCIPAL AXES STRESS ; 



WRITE (13,11) ' NEL ' , ' SX ' , ' SY ' , ' SZ ' f ' SXY ',' SXZ ',' SYZ ',' ANGLE ' 

DO 55 1=1, NEL 

READ (15, A) IEL , NOD , SX ( I ) ,SY(I) ,SZ(I) , SXY ( I ) , SXZ ( I ) , SYZ ( I ) 
THETA ( I ) = THETAF(I)A180.0 / PI 

WRITE (13,61) IEL , SX ( I ) ,SY(I) ,SZ(I) , SXY ( I ) , SXZ ( I ) , SYZ ( I ) , 

& THETA(I) 



A 

L=COS ( THETAF ( I ) ) 

M=SIN (THETAF( I ) ) 

A ( 1 , 1 ) = SX(I)A(LAA2) + SY(I)A(MAA2) + 2ASXY(I)ALAM 
A ( 1 , 2 ) = SX(I)A(MAA2) + SY(I)A(LAA2) - 2ASXY(I)alam 
A ( I , 3 ) = SZ(I) 

A ( I , 6 ) = -SX(I)ALAM + SY(I)amal + SXY(I)A(LAA2 - MAA2) 

A( 1 , 5 ) = SYZ(I)AM 4 SXZ ( I ) AL 
A( I , 4 ) = SYZ ( I ) AL - SXZ ( I ) AM 
55 CONTINUE 

WRITE (13,10) 'NEL' , 'SI' , 'S2' , 'S3' , 'S6' , 'S5' , 'S4' 

DO 80 1=1, NEL 

WRITE (13,60) I,A(I,1) ,A(I,2) ,A(I,3) ,A(I,6) ,A(I,5) ,A(I,4) 
80 CONTINUE 

A 



ASUBROUTINES FOR CALCULATING THE RELIABILITY ; 

A — — 



CALL BETA ( A, B, NEL) 

WRITE (5, A) 'INPUT ALFA' 
READ ( 5 , A ) ALF 
WRITE ( 13, A) 'ALFA : ' , ALF 
CALL ALFA ( A, B, NEL, ALF) 
CALL RELI ( A,B,NEL) 

CALL ELRE(B,NEL) 



10 


FORMAT 


( // , IX , A1 0 , 6 ( A7 ,4X) 


/ 


IX , 80 ( 


' ) ) 


11 


FORMAT 


( // ,1X,A10,7(A7,4X) 


/ 


IX , 80 ( ' - 


' ) ) 


15 


FORMAT 


(IX, ' THETA ( ' ,12, ' ) = 


' 


) 




25 


FORMAT 


(IX, ' THETA ( ' ,12, ' ) = 


' 


,F6. 1) 




60 


FORMAT 


( / , 4X , 16 , 6E1 1.3) 








61 


FORMAT 


( / ,4X,I6,6E11.3,3X, 


F6 


.1) 





A 



STOP 

END 



F+ + 4 - 4 - + + + + 4 - 



AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 



SUBROUTINE BETA ( A, B, NEL) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 



PARAMETER (ROWS= 100 , C0LS=6 ) 
REALA8 A(ROWS, COLS) ,B(ROWS, COLS) 
REAL X(COLS) ,X1C,X2C,X3C 



a PLY STRESS - STRENGTH, RATIO CALCULATION : 

A 



WRITE (13,10) 'NEL' , ' SRI ' , ' SR2 ' , ' SR3 ' , ' SR6 ' , ' SR5 ' , ' SR4 ' 
DO 90 1=1, NEL 

C STRENGTH PROPERTIES FOR GRAPHITE EPOXY, 

C 
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X(l) = 15.0E4 
X ( 2 ) = 0.4E4 

X(3) = X( 2 ) 

X ( 6 ) = 0.68E4 

X ( 4 ) = X(6) 

X( 5) = X ( 6 ) 

X1C = -15.0E4 
X2C = -2.46E4 
= Y?r 

IF ( A ( I , 1 ) • LT . 0 ) X ( 1 ) = X1C 
IF ( A ( 1 , 2 ) . LT . 0 ) X ( 2 ) = X2C 
IF ( A ( 1 , 3 ) . LT . 0 ) X( 3 ) = X3C 
DO 80 J=1 , COLS 

B ( I , J ) = ABS ( A( I f J ) / X(J) ) 

80 CONTINUE 

WRITE (13,60) I,B(I,1) ,B(I,2) ,B(I,3) ,B(I,6),B(I,5) ,B(I,4) 
90 CONTINUE 

10 FORMAT (// ,lX,A10,6(A7 f 4X) / 1X,80('-')) 

60 FORMAT ( / , 4X , 16 , 6E1 1 . 3 ) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE ALFA ( A , B , NEL , ALF ) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( ROWS= 1 00 , C0LS=6 ) 

REALA8 A (ROWS, COLS) ,B(ROWS, COLS) 

REAL AREA (100) , AMIN 

A 

* PLY STRESS - STRENGTH, RATIO RISED BY ALFA CALCULATION : 



WRITE (13,10) 'NEL' ,'V1' ,'V2' , 'V3‘ ,'V6' , 'V5' ,'V4* 

DO 90 1=1, NEL 

READ (17,*) I IEL , AREA ( I ) , AMIN 
AMIN=400 . 

WRITE (5,95) I IEL , AREA ( I ) , AMIN 
WRITE (3,95) I IEL, AREA ( I ) , AMIN 
DO 80 J = 1 , COLS 

* A ( I , J ) = B ( I , J ) * *ALF 

A ( I , J ) = ( B ( I , J ) **ALF) *AREA( I ) / AMIN 
80 CONTINUE 

WRITE (13,60)I,A(I,1),A(I,2),A(I,3),A(I,6),A(I,5),A(I,4) 
90 CONTINUE 

10 FORMAT (//,1X,A10,6(A7,4X) / lX,80('-')> 

60 FORMAT ( / , 4X, 16 , 6E11 . 3) 

95 FORMAT ( / / IX, ' AREA( ' , 12 , ' ) = ',2F10.2) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE RELI ( A,B,NEL) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( ROWS= 1 00 , C0LS=6 ) 

REAL* 8 A ( ROWS , COLS ) , B ( ROWS , COLS ) 

A 

* PLY STERSS TENSOR RELIABILITY CALCULATION (WEIBULL) : 



WRITE (13,10) ' NEL ' ,'R1','R2','R3','R6','R5','R4' 

DO 90 1=1, NEL 
DO 80 J = 1 , COLS 
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B(I,J) = EXP(-A(I,J ) )' 

80 CONTINUE 

WRITE (13,60) I,B(I,1) ,B(I,2> ,B(I,3) ,B(I,6) ,B(I,5) ,B(I,4) 

90 CONTINUE 

10 FORMAT (//,1X,A10,6(A7,4X) / 1X,80('-')) 

60 FORMAT ( / , 4X , 16 , 6E1 1 . 3 ) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE ELRE ( B , NEL ) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( ROWS=100 ,COLS=6 , C0LS1=7 ) 

REALAQ B(ROWS;COLS) ,C(ROWS,COLSl) ,RT 
A 

* ELEMENTS RELIABILITY CALCUL. , AS WELL AS 'TIBER" (Rl) AND "MATRIX" 
^RELIABILITY CALCUL . & TOTAL STRUCTURE RELIABILITY CALCUL. 

A 

WRITE (13,10) 'NEL' , ' Rl ' , ' R2 ' , ' R3 ' , ' R6 ' , ' R5 ' , ' R4 ' , ' RE' 

RT =1.0 

C ( NEL+1 , 1 ) =1 . 0 
DO 90 1=1, NEL 
C ( I , COLS1 ) = 1.0 
DO 80 J = 1 ,COLS 

C(I,J)=B(I,J) 

C ( I , COLS1 ) = C ( I , COLS1 ) A C(I,J) 

C(NEL+1 , J ) = C ( NEL+1 , J ) A C(I,J) 

80 CONTINUE 

RT = RT A C(I,C0LS1) 

WRITE (5,60) I,C(I,1),C(I,2),C(I,3) ,C(I,4) ,C(I,5) ,C(I,6),C(I,7) 
WRITE (13,60) I,C(I,1) ,C(I,2) ,C(I ,3) ,C(I,4 ) ,C(I,5) ,C(I,6) ,C( 1,7) 
90 CONTINUE 

WRITE (13,10) ' TOTAL REL RT1 ' , ' RT2 ' , ' RT3 ' , ' RT6 ' , ' RT5 ' , 

£ ' RT4 ' , ' RT' 

WRITE (13,65) C ( NEL+1 , 1 ) , C( NEL+1 , 2 ) ,C ( NEL+1 , 3 ) , C (NEL+1 , 4 ) , 

£ C ( NEL+1 , 5 ) ,C ( NEL+1 , 6 ) ,RT 

WRITE (5,70) RT 
WRITE (13,70) RT 

10 FORMAT ( // ,1X,A10,7(A7,4X) / 1X,95('-')) 

60 FORMAT ( / , 4X , 16 , 7E11 . 3 ) 

65 FORMAT ( / , 11X, 7E11 . 3 ) 

7 0 FORMAT (//, IX, 'TOTAL STRUCTURE RELIABILITY = ',E11.3 / IX , 4 0 ( ' = ' ) 
RETURN 
END 
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INPUT STRESS FN = THV3D.STR 
INPUT ELEM. ANGLE FN = THV3D.EL 
INPUT AREAS OF ELEMENTS FN =THV3D. AREA 
OUTPUT FN =THV3D.REL 
OUTPUT :THV3D. REL 



NEL 


SX 


SY 


SZ 


SXY 


SXZ 


SYZ 


ANGLE 


1 


-0. 279E+03 


0 . 391E+05 


-0. 286E+01 


-0 . 235E+04 


0 . 254E+01 


-0. 983E+01 


90.0 


2 


0 . 27 7E+03 


0 . 209E+05 


-0.268E+01 


0 . 235E+04 


0 . 254E+01 


-0 . 976E+01 


75.0 


3 


0.277E+03 


0. 209E+05 


- 0 . 268E+01 


0. 235E+04 


-0. 254E+01 


0 . 976E+01 


75.0 


4 


-0. 279E+03 


0. 391E+05 


-0 . 286E+01 


-0 . 235E+04 


-0 . 254E+01 


0.983E+01 


90.0 


5 


-0.316E+03 


0 . 390E+05 


-0 . 881E+02 


-0 . 235E+04 


0 . 742E-01 


-0.566E+00 


90.0 


6 


0 . 243E+03 


0. 209E+05 


-0. 754E+02 


0 . 235E+04 


0. 742E-01 


-0 . 332E+00 


75.0 


7 


0 . 243E+03 


0 . 209E+05 


-0. 754E+02 


0 . 235E+04 


-0 . 742E-01 


0. 332E+00 


75.0 


8 ■ 


-0. 316E+03 


0 . 390F.+0 5 


-0. 881E+02 


-0 . 235E+04 


-0. 742E-01 


0 . 566E+00 


90 . 0 


1 


-0.287E+03 


0 . 391E+05 


-0 . 169E+02 


-0 . 235E+04 


0. 494E+00 


-0. 216E+01 


90.0 


2 


0 . 270E+03 


0 . 209E+05 


-0 . 146E+02 


0. 235E+04 


0 . 494E+00 


-0. 196E+01 


75.0 


3 


0. 270E+03 


0. 209E+05 


-0. 146E+02 


0. 235E+04 


-0. 494E+00 


0 . 196E+01 


75.0 


4 ■ 


-0.287E+03 


0 . 391E+05 


-0.169E+02 


-0.235E+04 


-0 . 494E+00 


0 . 216E+01 


90.0 


NEL 


SI 


S2 


S3 


S6 


S5 


S4 




1 


0. 391E+05 


- 0 . 27 9E+03 


-0. 286E+01 


0. 235E+04 


-0.983E+01 


-0. 254E+01 




2 


0 . 207E+05 


0. 484E+03 


-0. 268E+01 


0 . 313E+04 


-0. 877E+01 


- 0 . 498E+01 




3 


0 . 207E+05 


0 . 484E+03 


-0 . 26BE+01 


0.313E+04 


0 . 877E+01 


0.498E+01 




4 


0.39 1E+0 5 


-0. 279E+03 


-0 . 286E+01 


0 . 235E+04 


0. 983E+01 


0.254E+01 




5 


0. 390E+05 


-0.31 6E+03 


-0 . 881E+02 


0 . 235E+04 


-0. 566E+00 


-0 . 742E-01 




6 


0. 207E+05 


0 . 451E+03 


-0 . 754E+02 


0 . 313E+04 


-0.301E+00 


-0 . 158E+00 




7 


0. 207E+05 


0 . 451E+03 


- 0 . 7 54E+02 


0. 313E+04 


0. 301E+00 


0. 158E+00 




8 


0 . 390E+05 


-0.316E+03 


-0 . 881E+02 


0 . 235E+04 


0 . 566E+00 


0.742E-01 




9 


0. 391E+05 


- 0 . 287E+03 


-0. 169E+02 


0 . 235E+04 


-0. 216E+01 


- 0 . 494E+00 




10 


0 . 207E+05 


0 . 477E+03 


-0. 146E+02 


0.31 3E+04 


-0.177E+01 


-0 . 985E+0 0 




11 


0 . 207E+05 


0. 477E+03 


-0 . 146E+02 


0 . 313E+04 


0.177E+01 


0. 985E+00 
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TOTAL STRUCTURE RELIABILITY = 0.807E+00 



APPENDIX D 



APPLICATIONS 

We are again considering the simply supported anisotropic square 
plate with uniform tensile load, described in the former appendices. 

Here, the meshing used is modeling a conventional repair of a 
damaged structure. We are thereby calculating the reliability for differ- 
ent configurations. 

Using the THP3D4 program, the verification problem (THV3D) is 
solved with different mesh sizes (i.e., different element sizes). Com- 
paring the results for the total reliability, w T e may conclude that the 
results are in good agreement with each other (Rx(THV3D) = 0.807, 
R t (THP3D4) = 0.824). 

The THM3D4 program analyzes the parent structure with a hole. 
The resulting total reliability is 0 (Rx = 0.0). In other w r ords, the 
structure will probably fail due to the given applied load. 

The THM3D4 program simulates a repair configuration using a 
w r eaker material (E-glass/epoxy instead of graphite/epoxy used for the 
original parent structure). The reliability decreases in comparison 
with the original undamaged structure, Rx = 0.1. 

To improve the total reliability, in the THA3D4 program w^e sim- 
ulate a repair configuration by using a different lay-up orientation 
([ 9 O 4 ]) for the patch and, as expected, the repair causes the total 
structure to be more reliable, Rx = 0.87. 
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In summary, from all the individual results for the computed reli- 
abilities listed in the programs below, we can interpret the critical 
element (see RE in the results) or the critical stress component (see 
RT1, RT2, etc.), which may represent the fiber or the matrix of the 
composite material. 
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INPUT 



FILE 



a ADINA-IN 

A 

a THP3D4 . IN EIGHT SYMMETRIC LAYER USING 3-D. 

A 

FILEUNITS LIST=B L0G=7 ECHO =7 
CONTROL PLOTUNIT=CM 
COLORS BCODE=BLUE 

A 

DATABASE CREATE 
WORKSTATION DEVICE=0 
A 

HEADING ' THP3D4 . IN ; aaaeiGHT LAYERS SOLUTION 90,75 DEG AAA' 

A 

MASTER IDOF = 000111 

PRINTOUT VOLUME = MAX I MUM IPD=1 IPRIC=0 P=STRAINS 
PORTHOLE FORMATTED=YES FILE=60 

A 

COORDINATES / ENTRIES NODE X Y Z 

1 - 20 . - 20 . 0.1 

2 0 . - 20 . 0.1 

3 20. -20. 0.1 

4 -20. 0. 0.1 

5 0. 0. 0.1 

6 20 . 0 . 0.1 

7 -20. 20. 0.1 

B 0. 20. 0.1 

9 20. 20. 0.1 

21 - 20 . - 20 . 0 . 

22 0 . - 20 . 0 . 

23 20. -20. 0. 

24 -20. 0. 0. 

25 0. 0. 0. 

26 20. 0. 0. 

27 -20. 20. 0. 

28 0. 20. 0. 

29 20. 20. 0. 

101 5. 0. 0. 

102 0 . 10 . 0 . 

103 5. 10. 0. 

104 15. 2.68 0. 

105 7.68 10. 0. 

106 15. 10. 0. 

107 15. -10. 0. 



ANOTE THE MAXES MEAN : 1=90DEG 2=15 3=75 4=45 5=-45 
A===================================================== 

MATERIAL 1 ORTHOTROPIC 181. E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 

GAB=7 . 17E05 7.17E05 7.17E05 

A 

SYSTEM 1 CYLINDRIC X=0 Y=0 Z=0 PHI=0 THETA=-90 . XSI=90. 



COORDINATES / ENTRIES NODES R THETA 



10 

11 

12 

13 

14 

15 

16 
17 



10 . 

10 . 

10 . 

10 . 

10 . 

10 . 

10 . 

10 . 



0. 

45. 

90. 

135. 

180. 

225. 

270. 

315. 



XL 
1 
1 
1 
1 
1 
1 
1 
1 



0 . 

0 . 

0. 

0. 

0. 

0. 

0. 

0 . 
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18 


10. 


360. 


0.1 


30 


10. 


0 . 


0 . 


31 


10. 


45. 


0 . 


32 


10. 


90. 


0 . 


33 


10. 


135. 


0 . 


34 


10. 


180. 


0 . 


35 


10. 


225. 


0 . 


36 


10. 


270. 


0 . 


37 


10. 


315. 


0 . 


30 


10. 


360. 


0 . 



A 



LINE 


CYLINDRIC 


ii 


10 


EL= 1 


MIDNODES-1 


NC-A 


LINE 


CYLINDRIC 


12 


11 


EL= 1 


MIDNODES= 1 


NC-A 


LINE 


CYLINDRIC 


13 


12 


EL=1 


HIDNODES= 1 


NC-A 


LINE 


CYLINDRIC 


14 


13 


EL= 1 


MIDNODES-1 


NC-A 


LINE 


CYLINDRIC 


15 


14 


EL = 1 


MIDNODES-1 


NC-A 


LINE 


CYLINDRIC 


16 


15 


EL=1 


MIDNODES-1 


NC-A 


LINE 


CYLINDRIC 


17 


16 


EL=1 


MIDNODES-1 


NC = A 


LINE 


CYLINDRIC 


18 


17 


EL= 1 


MIDNODES-1 


NC-A 


LINE 


CYLINDRIC 


31 


30 


EL-1 


MIDNODES-1 


NC-A 


LINE 


CYLINDRIC 


32 


31 


EL=1 


MIDNODES= 1 


NC-A 


LINE 


CYLINDRIC 


33 


32 


EL= 1 


MIDNODES-1 


NC-A 


LINE 


CYLINDRIC 


34 


33 


EL-1 


MIDNODES= 1 


NC-A 


LINE 


CYLINDRIC 


35 


34 


EL = 1 


MIDNODES-1 


NC-A 


LINE 


CYLINDRIC 


36 


35 


EL= 1 


MIDNODES- 1 


NC-A 


LINE 


CYLINDRIC 


37 


36 


EL= 1 


MIDNODES- 1 


NC-A 


LINE 


CYLINDRIC 


38 


37 


EL = 1 


MIDNODES-1 


NC-A 



A 



EGROUP 1 THREEDSOLID MATERIALS TABLES 

GVOLUME 7 4 18 17 27 24 38 37 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 4 1 11 10 24 21 31 30 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 1 2 12 11 21 22 32 31 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 2 3 13 12 22 23 33 32 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 3 6 14 13 23 26 34 33 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 6 9 15 14 26 29 35 34 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 9 8 16 15 29 28 36 35 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 8 7 17 16 28 27 37 36 EL1=1 EL2=1 EL3=4 NODES=27 

GSURFACE 7 8 28 27 EL1=1 EL2=4 NO=9 

GSURFACE 8 9 29 28 EL1=1 EL2=4 NO=9 

AXES- ORTHOTROPIC 1 N0DI=101 NODJ=103 N0DK=102 

AXES- ORTHOTROPIC 2 N0DI=101 N0DJ=104 N0DK=102 

AXES -ORTHOTROPIC 3 NODI=101 N0DJ=105 N0DK=102 

AXES- ORTHOTROPIC 4 N0DI=101 N0DJ=106 N0DK=102 

AXES-ORTHOTROPIC 5 N0DI=101 N0DJ=107 N0DK=102 

A STRES STABLE 112345678 21 26 27 

STRESSTABLE 1 21 

EDATA / ENTRIES EL TABLE MAXES 









1 


1 


1 


STEP 


4 


TO ' 


29 


1 


1 








2 


1 


3 


STEP 


4 


TO 


30 


1 


3 








3 


1 


3 


STEP 


4 


TO 


31 


1 


3 








4 


1 


1 


STEP 


4' 


‘to 


32 


1 


1 



A 



•102 



LOADS ELEMENT 



TO 



9 2 -3.0E04 

16 2 -3.0E04 



A 

FIXB 123 LINES / 8 28 

FIXB 23 SURFACES / 7 8 28 27 / 8 9 29 28 

A 

A 

EGROUP 2 THREEDSOLID MATERIALS TABLES 

GVOLUME 11 12 5 5 31 32 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 12 13 5 5 32 33 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 13 14 5 5 33 34 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 14 15 5 5 34 35 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 15 16 5 5 35 36 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 16 17 5 5 36 37 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 17 18 5 5 37 38 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 10 11 5 5 30 31 25 25 EL1=1 EL2=1 EL3=4 N0DES=27 

AXES-ORTHOTROPIC 1 NODI=101 N0DJ=103 N0DK=102 
AXES “ORTHOTROPIC 2 NODI=101 N0DJ=104 N0DK=102 
AXES “ORTHOTROPIC 3 N0DI=101 NODJ=105 NODK=102 
AXES-ORTHOTROPIC 4 NODI = 101 NODJ = 106' NODK = 102 
AXES-ORTHOTROPIC 5 NODI=101 N0DJ=107 N0DK=102 
^STRESSTABLE 112345678 21 26 27 
STRESSTABLE 1 21 

EDATA / ENTRIES EL TABLE MAXES 









1 


1 


1 


STEP 


4 


TO 


29 


1 


1 








2 


1 


3 


STEP 


4 


TO 


30 


1 


3 








3 


1 


3 


STEP 


4 


TO 


31 


1 


3 








4 


1 


1 


STEP 


4 


TO 


32 


1 


1 



A 

A 

EGZONE NAME=PARENT / 1 / 2 

FRAME HEAD I NG= UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH ZONE=PARENT N0DES=0 ELEMENT =1 HIDDEN=DASHED BCODE=ALL 
A FRAME HEADING=UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
AMESH ZONE=PATCH N0DES=0 ELEMENT =1 HIDDEN =DASHED BCODE=ALL 

A 

ADINA 

A 

END 
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>nooononnonnnnnnnnn > 



AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 



PROGRAM TRANS3 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

THIS PROGRAM TRANFORMS THE STRESSES OUTPUT FROM ADINA, 

EXPRESSED IN THE GLOBAL STRUCTURE COORDINATE SYSTEM, TO 
THE PRINCIPAL MATERIAL COORDINATE. AND THEN, CALCULATES THE 
RELIABILTY PER STRESS COMPONENT AND PER ELEMENT AND FINALLY 
EVALUATES THE TOTAL RELIABILITY. ALL THAT BY CALLING THE 
SUBROUTINE HEREIN. 

NEL : # OF TOTAL ELEMENTS. 

SX , SY SYZ s THE STRESS IN GLOBAL COORDINATE SYSTEM. 

THETA : THE LAYUP ANGLE PER ELEMENT. 

SI , S2 , . . . S6 : THE STRESS IN THE MATERIAL PRINCIPAL COORDINATE. 

R1,R2 ,R6 : THE RELIABILITY PER STRESS COMPONENT. 

RT1 , RT2 , . . . RT6 : THE TOTAL RELIABILTY OF STRUCTURE PER COPMPONENT. 
RE : THE RELIABILITY PER ELEMENT 

RT : THE TOTAL RELIABILITY. 

ALF : THE SHAPE FUNCTION PARAMER FOR WEIBULL 

BEFORE RUNNING THE PROGRAM VERIFY IN THE PROGRAM THE ANGLE THETA 
INPUT AND THE # OF ELEMENT. AND ALSO VERIFY THE MATERIAL PROEPRTIES 
GIVEN IN SUBROUTINE BETA. 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 



INTEGER NEL 

PARAMETER <NEL=64,PI = 3.1415927 ) 

PARAMETER ( R0WS=100 , C0LS=6 , COLSl=7 ) 

INTEGER ITHETA ( NEL) 

REAL* 8 A(ROWS,COLS) , B ( ROWS , COLS ) , C ( ROWS , COLS1 ) 

REAL SX(NEL) , SY ( NEL) ,SZ(NEL) ,SXY(NEL) , SXZ ( NEL ) ,SYZ(NEL) 
REAL L , M , THETA ( NEL ) , THETAF ( NEL ) , ALF 
CHARACTERS 20 FNAME , FNOUT , FNAME1 ,FNAME2 



WRITE ( 5 , * ) 'ENTER INPUT STRESS FN ?' 

READ ( 5, ' (A) ' ) FNAME 

WRITE (5, A) 'INPUT STRESS FN = ', FNAME 
OPEN ( UNIT=1 5 , FILE=FNAME , STATU S= ' OLD ' ) 

WRITE (5, A) 'ENTER INPUT ELEM. ANGLE FN? ' 

READ ( 5, ' ( A) ' ) FNAME2 

WRITE (5, A) 'INPUT ELEM. ANGLE FN = ' ,FNAME2 
OPEN ( UNIT= 16 ,FILE=FNAME2 , STATUS= ' OLD ' ) 

WRITE <5,A) 'ENTER AREAS OF ELEMENTS FN ?' 

READ (5, ' (A) ' ) FNAME1 

OPEN (UNIT= 17 ,FILE=FNAME1 , STATUS= ' OLD ' ) 

WRITE (5,A) ' ENTER OUTPUT FN ?' 

READ (5, ' (A) ' ) FNOUT 

OPEN ( UNIT= 13 ,FILE=FNOUT, STATUS = ' NEW ' ) 

WRITE (13, A) 'INPUT STRESS FN = ', FNAME 
WRITE (13, A) 'INPUT ELEM. ANGLE FN = ' ,FNAME2 
WRITE (13, a) 'INPUT AREAS OF ELEMENTS FN = ' ,FNAME1 
WRITE (-13, A) ' OUTPUT FN =', FNOUT 

A 

A INPUT OF ANGLE (THETA) FOR EA. ELEMENT ; 



DO 20 1=1, NEL 

READ (16, A) N, 12, 13, 14, 15, 16, ITHETAf I) ,I8,A9,A10,I11, 
& I 12, I 13, 114, 115, 116, 117, I IB, 119 

IF ( ITHETA ( I ) . EQ • 1 ) THETA ( I ) = 90.0 
IF ( ITHETA ( I ) . EQ . 3 ) THETA ( I ) = 75.0 
WRITE (5,25) N , THETA ( I ) 

THETAF ( I ) = THETA ( I ) A PI/180.0 
20 CONTINUE 
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A 

WRITE (13,*) ' OUTPUT :',FNOUT 

A 

* TRANSFORMATION OF STRESSES TO PRINCIPAL AXES STRESS ; 

A — — — — 

WRITE (13,11) ' NEL ' , ' SX ' , ' SY ' , ' SZ ',' SXY ' SXZ SYZ ',' ANGLE' 
DO 55 1=1, NEL 

READ (15,*) IEL , NOD, SX(I) ,SY(I) ,SZ(I) , SXY ( I ) , SXZ ( I ) , SYZ ( I ) 
THETA ( I ) = THETAF ( I ) * 1 8 0 . 0 / PI 

WRITE (13,61) IEL , SX( I ) , SY ( I ) , SZ ( I ) , SXY ( I ) , SXZ ( I ) , SYZ ( I ) , 

& THETA ( I ) 

A 

L=COS ( THETAF ( I ) ) 

M=SIN ( THETAF ( I ) ) 

A ( I , 1 ) = SX( I ) * ( L**2 ) + SY ( I ) * ( M**2 ) + 2*SXY(I)*L*M 
A( 1 , 2 ) = SX( I ) * ( M**2 ) -f SY ( I ) * ( L**2 ) - 2*SXY(I)*L*M 
A( 1 , 3 ) = SZ(I) 

A ( 1 , 6 ) = -SX ( I ) *L*M -f SY ( I ) *M*L + SXY(I)*(L**2 - M**2) 

A ( I , 5 ) = SYZ ( I ) *M + SXZ ( I ) *L 
A ( I , 4 ) = SYZ ( I ) *L - SXZ ( I ) *M 
55 CONTINUE 

WRITE (13,10) ' NEL ' ,'S1' ,'52' , ' S3 ' , ' S6 ' , ' S5 ' , ' S4 ' 

DO 80 1=1, NEL 

WRITE (13,60) I , A( I , 1 ) ,A(I,2) ,A(I,3) ,A(I,6) ,A(I,5) ,A(I,4) 

80 CONTINUE 

A 

* SUBROUTINES FOR CALCULATING THE RELIABILITY ; 

A 

CALL BETA(A,B,NEL) 

WRITE (5,*) 'INPUT ALFA' 

READ ( 5 , * ) ALF 

WRITE (13,*) 'ALFA : ' , ALF 

CALL ALFA ( A, B, NEL, ALF) 

CALL RELI ( A,B,NEL) 

CALL ELRE ( B , NEL ) 

A 



10 


FORMAT 


( // ,1X,A10,6(A7, 


4X ) 


/ 


IX , 80 ( 


11 


FORMAT 


( // , IX , A10 , 7 ( A7 , 


4X) 


/ 


IX , 80 ( 


15 


FORMAT 


( IX , ' THETA (',12, 


' ) = 


/ 


) 


25 


FORMAT 


(IX, ' THETA ( ' ,12, 


' ) = 


' 


, F6 . 1 ) 


60 


FORMAT 


( / , 4X , 16 , 6E1 1.3) 








61 


FORMAT 


( / ,4X,I6,6E11.3, 


3X,F6 


. 1) 



A 



STOP 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE BETA ( A , B , NEL ) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

• PARAMETER ( ROWS =100 , COLS =6 ) 

REALMS A(ROWS,COLS) ,B (ROWS, COLS) 

REAL X(COLS) ,X1C,X2C,X3C 

A 

A PLY STRESS - STRENGTH, RATIO CALCULATION : 

A 

WRITE (13,10) ' NEL ' , ' SRI ' , ' SR2 ' , ' SR3 ' , ' SR6 ' , ' SR5 ' , ' SR4 ' 
DO 90 1=1, NEL 

C STRENGTH PROPERTIES FOR GRAPHITE EPOXY, 

C 
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X ( 1 ) = 15.0E4 

X( 2) = 0.4E4 

X( 3 ) = X ( 2 } 

X ( 6 ) = 0.68E4 

X ( 4 ) = X ( 6 ) 

X( 5 ) = X(6) 

X1C = -15.0E4 
X2C = -2.46E4 
X3C = X2C 

IF (A(I,1) . LT . 0 ) X(l) = X1C 
IF ( A( I , 2 ) . LT . 0 ) X ( 2 ) = X2C 
IF ( A ( 1 , 3 ) . LT . 0 ) X ( 3 ) = X3C 
DO 80 J = 1 , COLS 

B ( I , J ) = ABS ( A( I , J ) / X( J) ) 

80 CONTINUE 

WRITE (13,60) I , B ( I , 1 ) ,B(I,2) , B ( I , 3 ) , B ( I , 6 ) , B ( I , 5 ) , B ( I , 4 ) 
90 CONTINUE 

10 FORMAT (// ,1X,A10,6(A7,4X) / 1X,80('-')) 

60 FORMAT ( / , 4X, 16 , 6E11 . 3 ) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE ALFA ( A , B , NEL , ALF ) 

AAAAA.AAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER (ROWS =100, COLS =6) 

REAL* B A ( ROWS , COLS ) , B ( ROWS , COLS ) 

REAL AREA ( 100) ,AMIN 

A 

* PLY STRESS - STRENGTH, RAT 10 RISED BY ALFA CALCULATION : 



WRITE (13,10) 'NEL' , 'VI' , 'V2' , 'V3' , 'V6' , 'V5' ,'V4' 

DO 90 1=1, NEL 

READ (17, *) IIEL,AREA(I) ,AMIN 
AMIN=400 . 

WRITE (5,95) I I EL , AREA ( I ) , AM I N 
WRITE (3,95) I IEL , AREA ( I ) , AMIN 
DO 80 J=1 ,CQLS 

* A ( I , J ) = B ( I , J ) a aaLF 

A ( I , J ) = ( B ( I ,J)**ALF)*AREA(I)/AMIN 
80 CONTINUE 

WRITE (13,60)I,A(I,1) ,A(I,2),A(I,3) , A ( I , 6 ) , A ( I , 5 ) , A ( I , 4 ) 
90 CONTINUE 

10 FORMAT ( // ,1X,A10,6(A7 ,4X) / 1X,80('-')) 

60 FORMAT ( / , 4X, 16 , 6E11 . 3 ) 

95 FORMAT ( //IX, 'AREA( ' ,12, ' ) = ',2F10.2) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

- SUBROUTINE RELI (A, B, NEL) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( ROWS=100 , C0LS=6 ) 

REAL* 8 A (ROWS, COLS) , B ( ROWS , COLS ) 

A 

* PLY STERSS TENSOR RELIABILITY CALCULATION ( WEIBULL) : 



WRITE (13,10) 'NEL ' , ' R1 ' , ' R2 ' , ' R3 ' , ' R6 ' , ' R5 ' , ' R4 ' 
DO 90 1=1, NEL 
DO 80 J= 1 , COLS 
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EXP ( -A( I , J ) ) ' * 



B(1,0) = 

80 CONTINUE 

WRITE (13,60) I , B ( 1 , 1 ) , B ( I , 2 ) , B ( 1 , 3 ) , B ( I , 6 ) , B ( 1 , 5 ) , B ( 1 , 4 ) 

90 CONTINUE 

10 FORMAT ( // ,1X,A10,6(A7,4X> / 1X,80('-')) 

60 FORMAT ( / , 4X, 16 , 6E11 . 3 ) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE ELRE(B,NEL) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( R0WS=100 , COLS =6 , C0LS1=7 ) 

REAL*8 B( ROWS, COLS) , C ( ROWS , COLS1 ) , RT 

A 

^ELEMENTS RELIABILITY CALCUL., AS WELL AS "FIBER " (Rl) AND "MATRIX" 

^ RELIABILITY CALCUL. 8 TOT AL STRUCTURE RELIABILITY CALCUL. 

A 

WRITE (13,10) 'NEL' , ' Rl ' , ' R2 ' , ' R3 ' , ' R6 ' , ' R5 # , ' R4 ' , ' RE ' 

RT =1.0 

C(NEL+1 ,1) =1 . 0 
DO 90 1=1, NEL 
C ( I , C0LS1 ) = 1.0 
DO 80 J=1 , COLS 

C(I,J)=B(I,J) 

C ( I , COLS1 ) = C ( I , COLS1 ) ^ C(I,J) 

C ( NEL+1 , J ) = C ( NEL+1 , J ) * C(I,J) 

80 CONTINUE 

RT = RT * C ( I , C0LS1 ) 

WRITE (5,60) I,C(I,1> ,C(I,2) ,C(I,3) ,C(I,4) ,C(I,5) ,C(I,6) ,C(I,7) 
WRITE (13,60) I , C ( 1 , 1 ) , C ( 1 , 2 ) ,C(I,3) , C ( 1 , 4 ) , C ( 1 , 5 ) ,C(I,6) , C ( 1 , 7 ) 
90 CONTINUE 

WRITE (13,10) 'TOTAL REL . : ' , ' RT1 ' , ' RT2 ' , ' RT3 ' , ' RT6 ' , ' RT5 ' , 

8 ' RT4 ' , ' RT ' 

WRITE (13,65) C ( NEL X 1 , 1 ) ,C (NEL+1, 2) ,C(NEL+1,3) ,C(NEL+1,4) , 

8 C( NEL+1, 5) ,C (NEL+1 ,6) , RT 

WRITE (5,70) RT 
WRITE (13,70) RT 

10 FORMAT ( // , IX, AID, 7 (A7,4X) / 1X,95('-')) 

60 FORMAT ( / , 4X , 16 , 7E1 1 . 3 ) 

65 FORMAT ( / , 11X , 7E1 1 . 3 ) 

70 FORMAT (//, IX ,' TOTAL STRUCTURE RELIABILITY = ',E11.3 /1X,40('=')) 
RETURN 
END 
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INPUT STRESS FN = THP3D4 . STP. 

INPUT ELEM. ANGLE FN = THP3D4.EL 

INPUT AREAS OF ELEMENTS FN =THP3D4 . AREA 

OUTPUT :THP3D4 . REL 



NEL 


SX 


SY 


sz 


SXY 


SXZ 


SYZ 


ANGLE 


1 


-0 . 277E+03 


0 . 391E+05 


0. 187E+02 


-0 . 235E+04 


0 . 909E+00 


0 . 622E+01 


90.0 


2 


0. 279E+03 


0. 209E+05 


0. 160E+02 


0 . 236E+04 


0 . 849E+00 


0 . 604E+01 


75.0 


3 


0 . 279E+03 


0. 209E+05 


0 . 160E+02 


0. 236E+04 


-0 . B49E+00 


-0.604E+01 


75.0 


4 


-0 . 277E+03 


0. 391E+05 


0. 187E+02 


-0 . 235E+04 


-0 . 909E+00 


- 0 . 622E+01 


90.0 


5 


-0. 278E+03 


0 . 391E+05 


0 . 248E+01 


-0 . 235E+04 


0 . 399E+00 


0 . 675E+01 


90.0 


6 


0 . 279E+03 


0 . 209E+05 


0 . 217E+01 


0 . 235E+04 


0 . 361E+00 


0 . 683E+01 


75.0 


7 


0 . 279E+03 


0. 209E+05 


0 . 217E+01 


0. 235E+04 


-0.361E+00 


- 0 . 683E+01 


75.0 


8 


-0 . 278E+03 


0 . 391E+05 


0 . 248E+0 1 


-0 . 235E+04 


-0 . 399E+00 


-0 . 675E+01 


90.0 


9 


-0 . 278E+03 


0 . 391E+05 


0 . 653E-01 


-0 . 235E+04 


0 . 552E+0 1 


-0 . 216E+02 


90.0 


10 


0 . 278E+03 


0 . 209E+05 


-0 . 746E+00 


0. 235E+04 


0 . 553E+01 


-0.216E+02 


75.0 


11 


0 . 278E+03 


0. 209E+05 


-0. 746E+00 


0. 235E+04 


-0 . 553E+0 1 


0 . 216E+02 


75.0 


12 


-0 . 278E+03 


0. 391E+05 


0.653E-01 


- 0 . 235E+04 


-0 . 552E+0 1 


0 . 216E+02 


90.0 


13 


-0. 278E+03 


0 . 391E+05 


0. 113E+00 


-0.235E+04 


0. 559E+01 


-0 . 210E+02 


90.0 


14 


0 . 278E+03 


0 . 209E+05 


- 0 . 692E+00 


0. 235E+04 


0. 559E+01 


-0. 210E+02 


75.0 


15 


0. 278E+03 


0. 209E+05 


- 0 . 692E+00 


0 . 235E+04 


-0 . 559E+01 


0 . 210E+02 


75.0 


16 


- 0 . 278E+03 


0 . 391E+05 


0. 113E+00 


- 0 . 235E+04 


-0 . 559E+0 1 


0 . 210E+02 


90.0 


17 


-0. 278E+03 


0. 391E+05 


0. 250E+01 


-0 . 235E+04 


-0 . 992E+0 0 


-0 . 4 34E+01 


90.0 


18 


0. 279E+03 


0 . 209E+05 


0. 218E+01 


0 . 235E+04 


-0.967E+00 


- 0 . 4 38E+01 


. 75.0 


19 


0 . 279E+03 


0 . 209E+ 05 


0 . 218E+ 01 


0 . 235E+04 


0. 967E+00 


0.438E+01 


75.0 


20 


-0. 278E+03 


0 . 391E+05 


0. 250E+01 


-0 . 235E+04 


0.992E+00 


0 . 434E+01 


90.0 


21 


-0. 278E+03 


0. 391E+05 


0. 190E+02 


-0 . 236E+04 


-0 . 608E+00 


-0 . 566E+01 


90.0 


22 


0. 277E+03 


0 . 209E+05 


0 . 163E+02 


0 . 235E+04 


-0 . 550E+00 


-0 . 573E+01 


75.0 


23 


0. 277E+03 


0 . 209E+05 


0. 163E+02 


0 . 235E+04 


0. 550E+00 


0 . 573E+01 


75.0 


24 


-0 . 278E+03 


0 . 391E+05 


0. 190E+02 


- 0 . 236E+04 


0 . 608E+00 


0 . 566E+01 


90.0 


25 


-0.31 8E+ 0 3 


0. 390E+05 


-0. 115E+03 


-0 . 235E+04 


-0. 598E+01 


-0 . 109E+01 


90.0 


26 


0 . 243E+03 


0 . 209E+05 


-0. 982E+02 


0 . 235E+04 


-0 . 607E+01 


-0 . 414E+00 


75.0 



108 



27 


0 . 243E+03 


0 . 209E+05 


- 0 . 982E+02 


28 


-0.31 8E+03 


0. 390E+05 


-0. 115E+03 


29 


-0.320E+03 


0 . 390E+05 


-0 . 115E+03 


30 


0.241E+03 


0 . 209E+05 


-0 . 980E+02 


31 


0 . 241E+03 


0 . 209E+05 


-0 . 980E+02 


32 


- 0 . 320E+03 


0 . 390E+05 


-0 . 1 15E+03 


1 


-0. 281E+03 


0 . 391E+05 


-0 . 594E+00 


2 


0 . 27 5E+03 


0 . 209E+05 


-0.107E+01 


3 


0 . 275E+03 


0 . 209E+05 


-0. 107E+01 


4 


-0. 281E+03 


0 . 391E+05 


-0. 594E+00 


5 


-0 . 281E+03 


0. 391E+05 


-0 . 646E+00 


6 


0 . 27 5E+03 


0 . 209E+05 


-0. 117E+01 


7 


0 . 275E+03 


0 . 209E+05 


-0.117E+01 


8 


- 0 . 281E+03 


0 . 391E+05 


-0.646E+00 


9 


-0 . 281E+03 


0.391E+05 


0 . 320E+00 


10 


0.27 5E+03 


0.209E+05 


0. 274E+00 


11 


0 . 275E+03 


0 . 209E+05 


0. 274E+00 


12 


-0. 281E+03 


0. 391E+05 


0. 320E+00 


13 


-0.282E+03 


0 . 39 1E+05 


0 . 64 1E-+01 


14 


0 . 274E+03 


0 . 209E+05 


0 . 552E+01 


15 


0. 274E+03 


0 . 209E+05 


0 . 552E+01 


16 


-0. 282E+03 


0. 391E+05 


0 . 64 1E+01 


17 


-0.294E+03 


0 . 391E+05 


-0. 422E+02 


18 


0. 263E+03 


0 . 209E+05 


- 0 . 361E+02 


19 


0. 263E+03 


0 . 209E+05 


-0. 361E+02 


20 


-0.294E+03 


0.391E+05 


-0 . 422E+02 


21 


-0.295E+03 


0. 390E+05 


-0. 416E+02 


22 


0 . 264E+03 


0. 209E+05 


-0. 356E+02 


23 


0. 264E+03 


0 . 209E+05 


-0.356E+02 


24 


-0 . 295E+03 


0 . 390E+05 


- 0 . 4 16E+02 



. 235E+04 


0 . 607E+01 


0 . 4 14E+00 


75.0 


. 235E+04 


0. 598E+01 


0. 109E+01 


90.0 


. 235E+04 


-0.634E+01 


-0.996E+00 


90.0 


. 235E+04 


-0 . 621E+01 


-0. 305E+00 


75.0 


. 235E+04 


0 . 621E+01 


0 . 305E+00 


75.0 


. 235E+04 


0.634E+01 


0 . 996E+00 


90.0 


. 235E+04 


0 . 205E+01 


-0. 802E+01 


90.0 


. 235E+04 


0 . 204E+0 1 


-0 . 800E+01 


75.0 


. 235E+04 


- 0 . 204E+01 


0 . 800E+01 


75.0 


. 235E+04 


-0 . 205E+01 


0 . 802E+01 


90.0 


. 235E+04 


0.206E+01 


-0 . 778E+01 


90.0 


. 235E+04 


0 . 206E+01 


-0.781E+01 


75.0 


. 235E+04 


-0.206E+01 


0 . 781E+01 


75.0 


. 235E+04 


-0 . 206E+01 


0 . 778E+01 


90.0 


. 235E+04 


-0 . 370E+00 


”0. 167E+01 


90.0 


. 235E+04 


-0 . 359E+00 


-0 . 170E+01 


75.0 


. 235E+04 


0 . 359E+00 


0. 170E+01 


75.0 


. 235E+04 


0.370E+00 


0. 167E+01 


90.0 


. 235E+04 


-0. 167E+00 


-0.211E+01 


90.0 


. 236E+04 


-0 . 170E+00 


-0 . 218E+01 


75.0 


. 236E+04 


0. 170E+00 


0 . 218E+01 


75.0 


. 235E+04 


0. 167E+00 


0.211E+01 


90.0 


. 236E+04 


-0 . 238E+0 1 


-0 . 128E+01 


90.0 


. 235E+04 


-0 . 229E+01 


-0. 481E+00 

\ 


75.0 


. 235E+04 


0 . 229E+01 


0 . 481E+00 


75.0 


. 236E+04 


0 . 238E+01 


0 . 128E+01 


90.0 


. 235E+04 


-0 . 202E+01 


-0.119E+01 


90.0 


. 236E+04 


-0 . 220E+01 


-0 . 446E+00 


75.0 


. 236E+04 


0 . 220E+01 


0 . 446E+00 


75.0 


. 235E+04 


0 . 202E+01 


0.119E+01 


90.0 



0 , 

- 0 , 

- 0 . 

0 

0 

-0 

- 0 , 

0 , 

0 

-0 

-0 

0 

0 

>0 

-0 
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0 

-0 

-0 

0 

0 

-0 

-0 

0 

0 

-0 

-0 
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0 

-0 
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25 


-0 . 281E+03 


0 . 39 1E+05 


0.629E+01 


26 


0 . 275E+03 


0.209E+05 


0.535E+01 


27 


0. 275E+03 


0 . 209E+05 


0 . 535E+01 


28 


-0. 281E+03 


0 . 391E+05 


0 . 629E+01 


29 


-0 . 280E+03 


0. 391E+05 


0 . 293E+00 


30 


0. 276E+03 


0 . 209E+05 


0.287E+00 


31 


0 . 276E+03 


0 . 209E+05 


0. 287E+00 


32 


-0.280E+03 


0 . 391E+05 


0 . 293E+0 0 


NEL 


SI 


S2 


S3 


1 


0. 391E+05 


-0. 277E+03 


0 . 187E+02 


2 


0. 207E+05 


0 . 485E+03 


0 . 160E+02 


3 


0.20 7E+05 


0 . 485E+03 


0 . 160E+02 


4 


0 . 391E+05 


-0 . 277E+03 


0 . 187E+02 


5 


0 . 391E+05 


-0 . 278E+03 


0. 248E+01 


6 


0 . 207E+05 


0 . 486E+03 


0.217E+01 


7 


0 . 20 7E+05 


0 . 486E+03 


0 . 217E+01 


8 


0 . 391E+05 


-0. 278E+03 


0 . 248E+0 1 


Q 


0 . 391E+05 


-0. 278E+03 


0 . 653E-01 


10 


0. 207E+05 


0 . 485E+03 


-0. 746E+00 


11 


0 . 207E+05 


0 . 485E+03 


-0.746E+00 


12 


0.391E+05 


-0.278E+03 


0 . 653E-01 


13 


C.391E+05 


-0.278E+03 


0 . 1 13E+00 


14 


0 . 207E+05 


0.485E-J-03 


-0 .692E+00 


15 


0. 207E+05 


0 . 485E+03 


-0 . 692E+00 


16 


0 . 391E+05 


- 0 . 278E+03 


0 . 1 13E+00 


17 


0 . 391E+05 


-0 . 278E+03 


0 . 250E+0 1 


18 


0 . 207E+05 


0 . 485E+03 


0 . 218E+0 1 


19 


0 . 207E+05 


0 . 485E+03 


0 . 218E+01 


20 


0.391E+05 


-0.278E+03 


0 . 250E+01 



235E+04 


0 . 323E+0 0 


0 . 1 79E+0 1 


90.0 


235E+04 


0 . 2 9 3 E+ 0 0 


0.211E+01 


75.0 


235E+04 


-0.293E+00 


-0 . 211E+01 


75.0 


235E+04 


-0 . 323E+00 


-0 . 179E+01 


90.0 


235E+04 


0. 130E+00 


0 . 284E+01 


90 . 0 


235E+04 


0.140E+00 


0 . 243E+01 


75.0 


235E+04 


-0.140E+00 


-0 . 243E+01 


75.0 


235E+04 


-0. 130E+00 


-0.284E+01 


90.0 


S6 


S5 


S4 




235E+04 


0 . 622E+0 1 


-0 . 909E+00 




313E+04 


0.605E+01 


0 . 743E+00 




313E+04 


-0 . 605E+01 


-0 . 743E+00 




, 235E+04 


-0.622E+01 


0 . 909E+00 




235E+04 


0 . 675E+0 1 


-0 . 399E+00 




, 313E+04 


0 . 669E+01 


0 . 142E+01 




, 313E+04 


-0 . 669E+0 1 


-0 . 142E+01 




. 235E+04 


-0.675E+01 


0 . 3 9 9 E+ 0 0 




, 235E+04 


-0 . 216E+02 


-0 . 552E+01 




. 313E+04 


-0. 194E+02 


-0 . 109E+02 




. 313E+04 


0 . 194E+02 


0. 109E+02 




. 235E+04 


0 . 216E+02 


0 . 552E+0 1 




. 235E+04 


-0 . 210E+02 


- 0 . 559E+0 1 




. 313E+04 


-0 . 189E+02 


-0. 10EE+02 




. 313E+04 


0 . 189E+02 


0. 108E+02 




, 235E+04 


0. 210E+02 


0 . 559E+01 




, 235E+04 


-0 . 434E+01 


0 . 992E+00 




. 3 13E+04 


-0 . 448E+01 


- 0 . 198E+00 




. 313E+04 


0 . 448E+01 


0 . 198E+00 




, 235E+04 


0 . 434E+01 


-0 . 992E+00 





0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

3B 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 



0 . 391E+05 - 
0. 207E+05 
0. 207E+05 
0 . 391E+05 - 
0 . 390E+05 - 
0.207E+05 
0 . 207E+05 
0 . 390E+05 - 
0 . 390E+05 - 
0 . 207E+05 
0 . 207E+05 
0. 390E+05 - 
0 . 391E+05 - 
0 . 207E+05 
0.207E+05 
0.391E+05 - 
0 . 391E+05 - 
0 . 2 07E+05 
0 . 207E+05 
0 . 391E+05 - 
0.391E+05 - 
0. 207E+05 
0 . 207E+05 
0 . 391E+05 - 
0 . 391E+05 - 
0.207E+05 
0 . 2 07E+05 
0.391E+05 - 
0 . 3 91E+05 - 
0 . 207E+05 



0.278E+03 
0 . 486E+03 
0 . 486E+03 
0 . 278E+03 
0.318E+03 
0 . 450E+03 
0 . 450E+03 
0. 318E+03 
0. 320E+03 
0 . 449E+03 
0. 449E+03 
0. 320E+03 
0 . 281E+03 
0. 482E+03 
0. 482E+03 
0. 281E+03 
0 . 281E+03 
0. 482E+03 
0 . 482E+03 
0. 281E+03 
0 . 281E+03 
0 . 482E+03 
0 . 482E+03 
0. 281E+03 
0 . 282E+03 
0 . 4 81E+03 
0 . 4 81E+03 
0. 282E+03 
0. 294E+03 
0 . 471E+03 



0. 190E+02 
0 . 163E+02 
0.163E+02 
0 . 190E+02 
-0. 115E+03 
-0 . 982E+02 
-0 . 982E+02 
-0.115E+03 
-0. 115E+03 
-0 . 980E+02 
-0 . 980E+02 
-0. 115E+03 
-0 . 594E+00 
-0. 107E+01 
-0.107E+01 
-0. 594E+00 
-0 . 646E+00 
-0.117E+01 
-0. 117E+01 
-0.646E+00 
0. 320E+00 
0. 274E+00 
0.274E+00 
0 . 320E+00 
0 . 641E+01 
0. 552E+01 
0 . 552E+01 
0 . 641E+01 
-0. 422E+02 
-0. 361E+02 



0 . 236E+04 
0 . 313E+04 
0 . 313E+04 
0 . 236E+04 
0 . 235E+04 
0 . 313E+04 
0 . 313E+04 
0 . 235E+04 
0 . 235E+04 
0. 313E+04 
0 . 313E+04 
0. 235E+04 
0 . 235E+04 
0. 313E+04 
0.313E+04 
0. 235E+04 
0 . 235E+04 
0.313E+04 
0 . 313E+04 
0 . 235E+04 
0 . 235E+04 
0 . 313E+04 
0 . 313E+04 
0.235E+04 
0 . 235E+04 
0 . 313E+04 
0 . 313E+04 
0 . 235E+04 
0 . 235E+04 
0.313E+04 



- 0 . 566E+01 
-0. 568E+01 

0 . 568E+01 
0 . 566E+01 
-0 . 109E+01 
-0.197E+01 
0 . 197E+01 
0. 109E+01 
-0 . 996E+00 
-0. 190E+01 
0. 190E+01 
0 . 996E+00 
-0 . 802E+01 
-0. 720E+01 
0. 720E+01 
0 . 802E+0 1 
-0 . 778E+01 
-0 . 701E+01 
0. 701E+01 
0. 778E+01 
-0.167E+01 
-0. 173E+01 
0 . 173E+01 
0 . 167E+01 
-0.211E+01 
-0 . 215E+01 
0 . 215E+01 
0.211E+01 

- 0 . 128E+01 
-0 . 106E+01 
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0 . 608E+00 
0 . 952E+00 
0. 952E+00 
0 . 608E+00 
0 . 598E+01 
0 . 575E+01 
0 . 575E+01 
0 . 598E+01 
0 . 634E+01 
0 . 592E+01 
0 . 592E+01 
0 . 634E+01 
0 . 205E+01 
0 . 404E+01 
0 . 404E+01 
0 . 205E+01 
0 . 206E+01 
0 . 401E+01 
0 . 401E+01 
0 . 206E+01 
0.370E+00 
0 . 927E-01 
0.927E-01 
0 . 370E+00 
0 . 167E+00 
0 . 401E+00 
0.401E+00 
0 . 167E+00 
0 . 238E+01 
0 . 208E+01 



51 


0 . 207E+05 


0 . 4 7 1E+0 3 


“0 . 361E+02 


52 


0.391E+05 


-0.294E+03 


- 0 . 422E+02 


53 


0 . 390E+05 


-0.295E+03 


-0 . 416E+02 


54 


0 . 207E+05 


0 . 470E+03 


-0 . 356E+02 


55 


0 . 207E+05 


0 . 47 0E+03 


-0 . 356E+02 


56 


0 . 390E+05 


-0.295E+03 


-0 . 416E+02 


57 


0.391 E*.- 0 5 


-0 . 281E+03 


0 . 629E+01 


58 


0 . 207E+05 


0 . 482E+03 


0 . 535E+01 


59 


0 . 207E+05 


0 . 482E+03 


0 . 535E+01 


60 


0 . 391E+05 


-0 . 281E+03 


0 . 629E+01 


61 


0 . 3 91E+05 


-0 . 280E+03 


0 . 293E+00 


62 


0 . 207E+05 


0 . 483E+03 


0 . 287E+00 


63 


0 . 207E+05 


0 . 483E+03 


0 . 287E+00 


64 


0 . 391E+05 


-0 . 280E+03 


0. 293E+00 


ALFA : 5 . 


.000000 






NEL 


R1 


R2 


R3 


1 


0. 999E+00 


0 . 1 00E+01 


0. 100E+01 


2 


0. 100E+01 


0 . 100E+01 


0.100E+01 


3 


0. 100E+01 


0. 100E+01 


0. 100E+01 


4 


0 . 999E+00 


0. 100E+01 


0. 100E+01 


5 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


6 


0. 100E+01 


0. 100E+01 


0 . 1 00E+01 


7 


0. 100E+01 


0.100E+01 


0 . 100E+01 


8 


0. 100E+01 


0. 100E+01 


0. 100E+01 


9 


0 . 999E+00 


0. 100E+01 


0. 100E+01 


10 


0. 100E+01 


0.100E+01 


0 . 100E+01 


11 


0.100E+01 


0. 100E+01 


0. 100E+01 


12 


0 . 999E+00 


0. 100E+01 


0. 100E+01 


13 


0. 100E+01 


0. 100E+01 


0 . 100E+01 



0.313E+04 


0. 106E+01 


“0 . 208E+01 




0 . 235E+04 


0 . 128E+01 


-0 . 238E+01 




0 . 235E+04 


-0. 119E+01 


0 . 202E+01 




0 . 313E+04 


-0 . 100E+01 


0 . 201E+01 




0 . 313E+04 


0 . 100E+01 


-0 . 201E+01 




0 . 235E+04 


0.119E+01 


-0.202E+01 




0. 235E+04 


0 . 179E+01 


-0 . 323E+00 




0 . 313E+04 


0 . 2 12E+01 


0 . 264E+00 




0 . 313E+04 


-0.212E+01 


“0 . 264E+00 




0 . 235E+04 


-0.179E+01 


0 . 323E+00 




0 . 235E+04 


0 . 284E+01 


-0.130E+00 




0 . 31 3E+04 


0.238E+01 


0 . 493E+00 




0 . 313E+04 


-0 . 238E+01 


-0 . 493E+00 




0 . 235E+04 


-0.284E+01 


0. 130E+00 




R6 


R5 


R4 


RE 


0. 100E+01 


0. 100E+01 


0 . 9 98E+00 


0.997E+00 


0 . 100E+01 


0 . 100E+01 


0 . 991E+00 


0 . 991E+00 


0. 100E+01 


0. 100E+01 


0 . 991E+00 


0 . 99 1E+00 


0. 100E+01 


0 . 100E+01 


0 . 998E+00 


0 . 997E+00 


0 . 100E+01 


0 . 100E+01 


0 . 9 99E+00 


0 . 9 98E+00 


0 . 100E+01 


0. 100E+01 


0 . 994E+00 


0*. 994E+00 


0 . 1 00E+01 


0 . 100E+01 


0 . 994E+00 


0 . 994E+00 


0 . 100E+01 


0 . 100E+01 


0 . 999E+00 


0 . 998E+00 


0. 100E+01 


0 . 100E+01 


0 .998E+00 


0 . 997E+00 


0 . 100E+01 


0 . 100E+01 


0 . 991E+00 


0 . 991E+00 


0. 100E+01 


0 . 100E+01 


0.991E+00 


0.991E+00 


0. 100E+01 


0. 100E+01 


0 . 998E+00 


0 . 997E+00 


0 . 100E+01 


0 . 100E+01 


0 . 999E+00 


0 . 998E+0 0 



14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 



0.100E+01 
0.100E+01 
0 . 100E+01 
0 . 999E+00 
0. 100E+01 
0. 100E+01 
0 . 999E+00 
0.100E+01 
0. 100E+01 
0.100E+01 
0. 100E+01 
0.999E400 
0 . 100E+01 
0 . 100E+01 
0 . 999E+00 
0. 100E+01 
0 . 1 OOE+Ol 
0 . 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0 . 100E+01 
0.100E+01 
0 . 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0.10 0E+01 
0. 100E+01 
0.100E+01 



0.100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0.100E+01 
0. 100E+01 
0. 100E+01 
0 . 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0 . 1 00E+01 
0 . 100E+01 
0 . 100E+01 
0.100E+01 
0 . 100E+01 
0 . 1 00E+01 
0 . 1 00E+01 



0. 100E+01 
0 . 100E+01 
0.100E+01 
0. 100E+01 
0 . 100E+01 
0.100E+01 
0. 100E+01 
0 . 100E+01 
0 . 1 00E+01 
0.10 OE+0 1 
0. 100E+01 
0 . 100E+01 
0. 100E+01 
0. 100E-^01 
0.100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0.100E+01 
0. 100E+01 
0. 100E+01 
0.100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0 . 100E+01 
0 . iOOE+O 1 
0.100E+01 
0.100E+01 
0. 100E+01 



0 . 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0 . 100E+01 
0. 100E+01 
0.100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0 . 1 00E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0 . 1 00E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0.100E+01 
0 . 100E+01 
0. 100E+01 
0. 100E+01 
0 . 100E+01 
0. 100E+01 
0. 100E+01 
0 . 100E+01 
0. 100E+C1 



0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0 . 100E+01 
0.100E+01 
0 . 100E+01 
0. 100E+01 
0 . 1 OOE+Ol 
0.100E+01 
0.100E+01 
0 . 100E+01 
0. 100E+01 
0 . 100E+01 
0.100E+01 
0. 100E+01 
0. 100E+01 
0. 100E+01 
0. 1 OOE+O 1 
0. 100E+01 
0. 1 OOE+O 1 
0. 100E+01 
0. 100E+01 
0 . 1 OOE+O 1 
0. 100E+01 
0. 100E+01 
0 . 1 OOE+O 1 
0. 100E+01 
0. 100E+01 
0. 100E+01 



0.994E+00 
0 . 994E+00 
0 . 999E+00 
0.998E+00 
0 . 991E+00 
0 . 991E+00 
0. 998E+00 
0 . 999E+00 
0.994E+00 
0. 994E+00 
0 . 999E+00 
0 . 998E+00 
0.991E+00 
0. 991E+00 
0 . 998E+00 
0. 999E+00 
0 . 994E+00 
0.994E+00 
0 . 999E+00 
0. 100E+01 
0 . 998E+00 
0 . 998E+00 
0. 100E+01 
0 . 100E+01 
0.998E+00 
0 . 998E+00 
0 . 1 OOE+Ol 
0 . 100E+01 
0. 998E+00 
0. 99BE+00 



0 . 994E+00 
0 . 994E+00 
0. 998E+00 
0 . 997E+00 
0 . 991E+00 
0 . 991E+00 
0 . 9 97E+00 
0 . 998E+00 
0 . 994E+00 
0 . 994E+0 0 
0 . 998E+00 
0 . 997E+00 
0 . 991E+00 
0 . 991E+00 
0 . 997E+00 
0 . 998E+00 
0 . 994E+00 
0 . 994E+00 
0 . 9 98E+00 
0 . 999E+00 
0 . 998E+00 
0. 998E+00 
O'. 999E+00 
0 . 999E+00 
0. 998E+00 
0 . 998E+00 
0.999E+00 
0. 999E+00 
0. 99BE+00 
0. 99BE+00 
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44 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


45 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


46 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


47 


0 . 100E+01 


0 . 100E+01 


0. 100E+01 


48 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


49 


0 . 1 OOE+Ol 


0 . 100E+01 


0. 100E+01 


50 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


51 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


52 


0. 100E+01 


0 . 100E+01 


0.100E+01 


53 


0. 100E+01 


0. 100E+01 


0. 100E+01 


54 


0. 100E+01 


0. 100E+01 


0. 100E+01 


55 


0. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


56 


0. 100E+01 


0 . 1 OOE+Ol 


0 . 1 OOE+Ol 


57 


0. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


58 


0. 100E+01 


0. 100E+01 


0. 100E+01 


59 


0. 100E+01 


0. 100E+01 


0 . 1 OOE+O 1 


60 


0. 100E+01 


0. 100E+01 


0. 100E+01 


61 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


62 


0. 100E+01 


0 . 1 OOE+Ol 


0. 100E+01 


63 


0. 100E+01 


0 . 1 OOE+Ol 


0. 100E+01 


64 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


TOTAL REL. 


RT1 


RT2 


RT3 




0 . 991E+00 


0. OOOE+OO 


0. OOOE+OO 



TOTAL STRUCTURE RELIABILITY = 0.824E+00 



. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 999E+00 


. 100E+01 


0. 100E+01 


0. 100E+01 


0. 999E+00 


. 1 OOE+O 1 


0. 100E+01 


0 . 998E+00 


0. 998E+00 


. 100E+01 


0. 100E+01 


0 . 998E+00 


0. 998E+00 


. 100E+01 


0 . 100E+01 


0 . 100E+01 


0 . 999E+00 


. 100E+01 


0 . 100E+01 


0. 100E+01 


0. 999E+00 


.iooe:+oi 


0 . 1 OOE+Ol 


0 . 998E+00 


0 . 998E+00 


. 100E+01 


0. 100E+01 


0 . 998E+00 


0 . 998E+00 


. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


0. 999E+00 


. 100E+01 


0 . 1 OOE+O 1 


0 . 100E+01 


0 . 999E+00 


. 100E+01 


0. 100E+01 


0 . 998E+00 


0 . 998E+00 


. 1 OOE+O 1 


0. 100E+01 


0. 998E+00 


0.998E+00 


. 100E+01 


0 . 100E+01 


0. 100E+01 


0. 999E+00 


. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


0. 999E+00 


. 100E+01 


0. 100E+01 


0. 998E+00 


0. 998E+00 


. 100E+01 


0 . 100E+01 


0 . 998E+00 


0. 998E+00 


. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


0. 999E+00 


. 100E+01 


0 . 1 OOE+O 1 


0. 100E+01 


0. 999E+00 


. 100E+01 


0. 100E+01 


0 . 998E+00 


0 . 998E+00 


. 100E+01 


0. 100E+01 


0 . 998E+00 


0. 998E+00 


. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 999E+00 


RT6 


RT5 


RT4 


RT 


0. OOOE+OO 


0. OOOE+OO 


0. OOOE+OO 


0 . 824E+00 



0 , 

0 , 

0 , 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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INPUT 



FILE 



A ADI N A - IN 

A 

A THD3D4 . IN EIGHT SYMMETRIC LAYER USING 3-D. 

A 

FILEUNITS LIST=8 L0G=7 ECH0=7 
CONTROL PLOTUNIT=CM 
COLORS BCODE=BLUE 

A 

DATABASE CREATE 
WORKSTATION DEVICE=0 
A 

HEADING ' THD3D4 . IN ; aaa£IGTH LAYERS C90,753 , DAMAGED STRUCTURE aaa' 

A 

MASTER IDOF=0001 1 1 

PRINTOUT VOLUME = MAX I MUM IPD=1 IPRIC=0 P=STRAINS 
PORTHOLE FORMATTED=YES FILE=60 

A 



OORDINATES / ENTRIES NODE 


X 


Y 


Z 


1 


-20. 


-20. 


0.1 


2 


0 . 


-20. 


0.1 


3 


20. 


-20. 


0.1 


4 


-20. 


0 . 


0.1 


5 


0 . 


0 . 


0.1 


6 


20. 


0 . 


0.1 


7 


-20. 


20. 


0.1 


8 


0 . 


• 20. 


0.1 


9 


20. 


20. 


0.1 


21 


-20. 


-20. 


0 . 


22 


0 . 


-20. 


0 . 


23 


20. 


-20. 


0 . 


24 


-20. 


0 . 


0 . 


25 


0 . 


0 . 


0 . 


26 


20. 


0 . 


0 . 


27 


-20. 


20. 


0 . 


28 


0 . 


20. 


0 . 


29 


20. 


20. 


0 . 


101 


5. 


0 . 


0 . 


102 


0 . 


10. 


0 . 


103 


5. 


10. 


0 . 


104 


15. 


2. 


68 0. 


105 


7.68 


10. 


0 . 


106 


15. 


10. 


0 . 


107 


15. 


-10. 


0 . 


'NOTE THE MAXES MEAN : 1=90DEG 2=15 


3 = 75 


4 = 45 5 = - 4 5 



MATERIAL 1 ORTHOTROPIC 1EJ1.E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 

GAB=7 . 17E05 7.17E05 7.17E05 

A 

SYSTEM 1 CYLINDRIC X=0 Y = 0 Z = 0 PHI = 0 THETA=-90. XSI=90. 



COORDINATES / ENTRIES NODES 


R THETA 


XL 


10 


10. 


0 . 


0.1 


11 


10. 


45. 


0.1 


12 


10. 


90. 


0.1 


13 


10. 


135. 


0.1 


14 


10. 


180. 


0.1 


15 


10. 


225. 


0.1 


16 


10. 


270. 


0.1 


17 


10. 


315. 


0.1 
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IB 


10. 


360. 


0.1 


30 


10. 


0 . 


0 . 


31 


10. 


45. 


0 . 


32 


10. 


90. 


0 . 


33 


10. 


135. 


0 . 


34 


10. 


1B0 . 


0 . 


35 


10. 


225. 


0 . 


36 


10. 


270. 


0 . 


37 


10. 


315. 


0 . 


3B 


10. 


360. 


0 . 



A 



LINE 


CYLINDRIC 


11 


10 


EL = 1 


MIDN0DES=1 


NC=A 


LINE 


CYLINDRIC 


12 


11 


EL=1 


MIDN0DES=1 


NC=A 


LINE 


CYLINDRIC 


13 


12 


EL=1 


MIDNODES =1 


NC=A 


LINE 


CYLINDRIC 


14 


13 


EL = 1 


MIDNODES =1 


NC = A 


LINE 


CYLINDRIC 


15 


14 


EL = 1 


MIDNODES =1 


NC=A 


LINE 


CYLINDRIC 


16 


15 


EL = 1 


MIDNODES= 1 


NC = A 


LINE 


CYLINDRIC 


17 


16 


EL= 1 


MIDN0DES=1 


NC=A 


LINE 


CYLINDRIC 


IB 


17 


EL = 1 


MIDNODES= 1 


NC=A 


LINE 


CYLINDRIC 


31 


30 


EL=1 


MIDNODES = 1 


NC = A 


LINE 


CYLINDRIC 


32 


31 


EL= 1 


MIDNODES =1 


NC = A 


LINE 


CYLINDRIC 


33 


32 


EL= 1 


MIDN0DES=1 


NC=A 


LINE 


CYLINDRIC 


34 


33 


EL = 1 


MIDNODES =1 


NC = A 


LINE 


CYLINDRIC 


35 


34 


EL=1 


MIDNODES= 1 


NC=A 


LINE 


CYLINDRIC 


36 


35 


EL = 1 


MIDNODES =1 


NC = A 


LINE 


CYLINDRIC 


37 


36 


EL= 1 


MIDNODES = 1 


NC = A 


LINE 


CYLINDRIC 


38 


37 


EL = 1 


MIDNODES =1 


NC=A 



A 



EGROUP 1 THREEDSOLID MATERIALS TABLES 

GVOLUME 7 4 18 17 27 24 38 37 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 4 1 11 10 24 21 31 30 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 1 2 12 11 21 22 32 31 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 2 3 13 12 22 23 33 32 EL1 = 1 EL2 = 1 EL3 = 4 N0DES = 27 

GVOLUME 3 6 14 13 23 26 34 33 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 6 9 15 14 26 29 35 34 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 9 8 16 15 29 28 36 35 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 0 7 17 16 28 27 37 36 EL1=1 EL2=1 EL3=4 N0DES=27 

GSURFACE 7 8 28 27 EL1=1 EL2=4 N0=9 

GSURFACE 8 9 29 28 EL1=1 EL2=4 N0=9 

AXES -ORTHOTROPIC 1 NODI=101 NODJ=103 NODK=102 

AXES -ORTHOTROPIC 2 N0DI=101 N0DJ=104 NODK=102 

AXES -ORTHOTROPIC 3 NODI=101 N0DJ=105 NODK=102 

AXES -ORTHOTROPIC 4 N0DI=101 NODJ=106 N0DK=102 

AXES-ORTHOTROPIC 5 N0DI=101 N0DJ=107 N0DK=102 



ASTRESSTABLE 1 


1 2 


3 4 5 


6 7 8 


STRES STABLE 1 


21 






EDATA / 


ENTRIES 


EL 


TABLE 


MAXES 








1 


1 


1 


STEP 


4 


TO * 


29 


1 


1 


, 






2 


1 


3 


STEP 


4 


TO 


30 


1 


3 








3 


1 


3 


STEP 


4 


TO 


31 


1 


3 








4 


1 


1 


STEP 


4 


"TO 


32 


1 


1 



A 
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LOADS ELEMENT 

9 2 -3.0E04 

TO 

16 2 -3.0E04 

A 

FIXB 123 LINES / 8 28 

FIXB 23 SURFACES / 7 8 28 27 / 8 9 29 28 

A 

A 

EGZONE NAME=PARENT / 1 

FRAME HEADING =UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH ZONE=PARENT N0DES = 0 ELEMENT =1 H I DDEN = D A S HED BCODE=ALL 

A 

ADINA 

A 

END 
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INPUT STRESS FN = THD3D4.STR 
INPUT ELEM. ANGLE FN = THD3D4.EL 
INPUT AREAS OF ELEMENTS FN =THD3D4 . ARE 
OUTPUT FN =THD3D4 . REL 
OUTPUT :THD3D4.REL 



NEL 


SX 


5Y 


SZ 


SXY 


SXZ 


SYZ 


ANGLE 


1 


-0. 136E+04 


0 . 660E+05 


0. 477E+02 


-0 . 669E+04 


0 . 161E+01 


-0.203E+02 


90.0 


2 


-0.220E+04 


0. 199E+05 


0 . 4 96E+02 


-0. 279E+04 


-0.422E+01 


-0 . 381E+02 


75.0 


3 


-0.220E+04 


0. 199E+05 


0 . 496E+02 


-0 . 279E+04 


0 . 422E+0 1 


0 . 381E+02 


75.0 


4 


-0 . 136E+04 


0.660E+05 


0. 477E+02 


-0.669E+04 


-0. 161E+01 


0.203E+02 


90.0 


5 


0 . 322E+04 


0. 725E+05 


0. 120E+03 


-0.68BE+04 


0 . 267E+02 


0.119E+02 


90.0 


6 


0 . 3 26E+04 


0 . 291E+05 


0. 119E+03 


-0. 245E+03 


0.105E+01 


0. 110E+02 


75.0 


7 


0 . 326E+04 


0. 291E+05 


0. 119E+C3 


-0.24 5E+03 


-0. 105E+01 


-0.110E+02 


75.0 


6 


0 . 322E+04 


0. 725E+05 


0. 120E+03 


-0 .688S+04 


-0 . 267E+02 


-0. 119E+02 


90 . 0 


c 


0 .409E+04 


0. 343E+05 


-0. 165E+03 


- 0 . 54 2E+04 


0. 774E+02 


-0.334E+02 


90. 0 


10 


0 . 3 19E+04 


0 . 478E+04 


-0 . 113E+03 


-0. 427E+04 


0. 115E+02 


0.387E+01 


75.0 


11 


0.319E+04 


0. 478E+04 


-0 . 113E+03 


-0 . 427E+04 


-0 . 115E+ 0 2 


-0. 387E+01 


75.0 


12 


0 . 4 09E+04 


0. 343E+05 


-0 . 165E+03 


-0.54 2E+04 


-0. 774E+02 


0.334E+02 


90 . 0 


13 


0. 168E+04 


0. 145E+05 


-0. 985E+02 


0. 161E+04 


-0.421E+02 


-0.4 0 9E-*- 02 


90.0 


14 


0. 376E+04 


0. 240E X 05 


-0. 156E+03 


0. 792E+04 


0.97 7E+0 1 


-0. 744E+01 


75.0 


15 


C. 376E+04 


0.240E+05 


-0 . 156E+03 


C. 792E+04 


-0 . 977E+01 


0. 744E+01 


75.0 


ie 


C.168E+04 


0.145E-05 


-0 . 985E+02 


0. 181E+04 


0 . 421E+C2 


0 . 4G9E+02 


90. 0 


17 


0. 132E+04 


C.444E+05 


0. 14SE+03 


-0. 180E+04 


-0 . 172E+02 


0. 362E+02 


90.0 


ie 


0 . 271E+04 


0 . 30 1E+05 


0.152E+03 


0. 533E+G4 


0. 516E+01 


0 . 432E+02 


• 75.0 


is 


0. 271E+04 


0 . 3C1E+C5 


0 . 152E+03 


0. 533E+04 


- 0 . 518E+01 


-0 . 4 32E+02 


75.0 


20 


0. 132E^04 


0.44 4E X 0 5 


0. 149E+03 


-0. 180E+04 


0. 172E+02 


-0 . 362E+02 


90.0 


21 


-0. 342E+02 


0 . 607E+Q 5 


0 . 257E+02 


-0. 126E+04 


-0.96SE+01 


0 . 169E+0 2 


90.0 


22 


0. 234E+04 


0 . 458E+05 


0. 148E+02 


0 . 948E+04 


0. 214E+01 


0 . 727E+01 


75.0 


23 


0.234E+04 


0 . 4 56E+0 5 


C. 148E^02 


0. 94 8E+0 4 


-0.214E+01 


-0 . 727E+01 


75.0 


24 


-0. 342E+02 


0. 607E+05 


0 . 257E+02 


-0. 126E+04 


0. 969E+01 


- 0 . 169E+02 


90.0 


25 


-0. 440E+04 


0. 238E+05 


0. 16SE+03 


-0. 125E+04 


-0.B36E+O1 


0.377E+02 


90.0 



26 


-0. 456E+04 


0. 950E+04 


0. 196E+03 


27 


-0. 456E+04 


0. 950E+04 


0 . 196E+03 


21 


0 . 238E+05 


0. 169E+03 


-0 . 125E+04 


29 


-0 . 410E+04 


0 . 265E+05 


-0.145E+03 - 


30 


-0.447E+04 


0 . 858E+04 


-0 . 126E+03 - 


31 


-0.447E+04 


0 . 858E+04 


-0.126E+03 - 


32 


-0 . 410E+04 


0 . 265E+05 


-0.145E+03 - 


NEL 


SI 


S2 


S3 


1 


0 . 660E+05 


-0. 136E+04 


0 . 4 77E+02 


2 


0. 171E+05 


0.677E+03 


0.496E+02 


3 


0. 171E+05 


0 . 677E+03 


0. 496E+02 


4 


0.660E+05 


-0. 136E+04 


0 . 4 77E+02 


5 


0. 725E+05 


0 . 322E+04 


0. 120E+03 


6 


0. 272E+05 


0. 511E+04 


0 . 119E+03 


7 


0 . 272E+C5 


0.511E+04 


0 . 119E+03 


8 


0. 725E+05 


0. 322E+04 


0 . 120E+03 


q 


0. 343E+05 


0.40 9E+04 


-0 . 165E+0 3 


10 


0. 254E+04 


0. 543E+04 


-0. 113E+03 


11 


0. 254E+04 


0. 543E+04 


-0 . 113E+03 


12 


0. 343E+05 


0 . 4 0SE+04 


-0. 165E+03 


13 


0. 145E+05 


0. 168E+04 


- 0 . 985E+02 - 


14 


0 . 266E+0 5 


0. 116E+04 


-0 . 156E+03 - 


15 


0. 266E+05 


0. 116E+04 


- 0 . 156E+03 - 


16 


0. 145E+05 


0. 168E+04 


- 0 . 985E+02 - 


17 


0 . 444E+05 


0.132E+04 


0. 149E+03 


18 


0. 309E+05 


0. 188E+04 


0 . 152E+0 3 


19 


0. 309E+05 


0 . 188E+04 


0. 152E+03 


20 


0 . 444E+05 


0 . 132E+04 


0. 149E+03 


21 


0. 607E+05 


-0.342E+02 


0 . 257E+02 



298E+03 


-0. 589E+01 


-0. 127E+01 


75.0 


298E+03 


0 . 589E+01 


0 . 12 7E+0 1 


75.0 


836E+01 


-0. 377E+02 


-0.469E-02 


90.0 


193E+04 


-0 . 202E+02 


0.241E+02 


90.0 


600E+03 


-0 . 94 3E+01 


-0 . 503E+0 1 


75.0 


600E+03 


0.943E+01 


0 . 503E+01 


75.0 


1 93E+04 


0 . 202E+02 


-0.24 1E+02 


90.0 


S6 


S5 


S4 




. 669E+04 


- 0 . 2 03E+02 


-0.161E+01 




, 795E+04 


-0. 379E+02 


-0 . 578E+01 




, 79 5E+04 


0. 379E+02 


0.578E+01 




, 669E+04 


0 . 203E+0 2 


0. 161E+01 




, 688E+04 


0. 119E+02 


-0 . 267E+02 




, 666E+04 


0 . 1 09E+0 2 


0. 184E+01 




, 666E+04 


-0.109E+02 


-0. 184E+01 




, 688E+04 


-0. 119E+02 


0 . 267E+02 




542E+04 


-0. 334E+02 


-0 . 774E+02 




410E+04 


0 . 671E+01 


-0. 101E+02 




. 410E+04 


-0 . 671E+01 


0. 101E+02 




, 542E+04 


0. 334E+02 


0 . 774E+02 




, 181E+04 


-0. 409E+02 


0 . 421E+02 




, 1 80E+04 


- 0 . 465E+01 


-0 . 114E+02 




, 180E+04 


0. 465E+01 


0. 114E+02 




, 181E+04 


0. 409E+02 


-0.421E+02 




, 180E+04 


0. 362E+02 


0. 172E+02 




, 223E+04 


0. 430E+02 


0.617E+01 




, 223E+04 


-0.430E+02 


-0 . 617E+01 




180E+ 04 


-0. 362E+02 


-0 . 172E+02 




126E+04 


0. 169E+02 


0. 969E+01 





0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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22 


0 . 476E+05 


0. 507E+03 


0 . 148E+02 


23 


0. 476E+05 


0. 507E+03 


0. 148E+02 


24 


0 . 607E+05 


-0. 342E+02 


0 . 257E+02 


25 


0. 238E+05 


-0. 440E+04 


0. 169E+03 


26 


0. 870E+04 


-0. 377E+04 


0. 196E+03 


27 


0.870E+04 


-0. 377E+04 


0. 196E+03 


28 


0 . 169E+03 


0. 238E+05 


-0. 125E+04 


29 


0 . 265E+05 


-0.41 0E+O4 


-0.145E+03 


30 


0. 741E+04 


-0. 330E+04 


-0.126E+03 


31 


0. 741E+04 


-0 . 330E+04 


-0.126E+03 


32 


0 . 265E+05 


-0. 410E+04 


-0 . 145E+03 



ALFA : 5.000000 

NEL R1 R2 R3 



1 


0. 993E+00 


0 . 100E+01 


0 . 100E+01 


2 


0. 100E+01 


0. 100E+01 


0. 100E+01 


3 


0. 100E+01 


0. 100E+01 


0. 100E+01 


4 


0. 993E+00 


0. 100E+01 


0. 100E+01 


5 


0. 992E+00 


0.906E+00 


0 . 100E+01 


6 


0. 100E+01 


0. 372E+00 


0. 100E+01 


7 


0. lOOE^Ol 


0. 372E+00 


0. 100E+01 


e 


0. 9S2S+00 


0. 906E+00 


0. 100E+01 


9 


0.10 OE+0 1 


0 . 614E+00 


0. 100E+01 


10 


0. 100E+01 


0. 133E+00 


0 . 100E+01 


n 


0. 100E+01 


0.133E+00 


0. 100E+01 


12 


0. 100E+01 


0 . 614E+0 0 


0. 100E+01 


13 


0. 100E+01 


0 . 996E+0 0 


0. 100E+01 


14 


0. 100E+01 


0. 999E+00 


0. 100E+01 


15 


0. 100E+01 


0. 999E+00 


0. 100E+01 


16 


0. 100E+01 


0. 996E+00 


0. 100E+01 



265E+04 


0. 757E+01 


-0. 187E+00 


265E+04 


-0. 757E+01 


0. 187E+00 


126E+04 


-0. 169E+02 


-0 . 969E+01 


125E+04 


0 . 377E+02 


0 . 836E+01 


326E+04 


-0 . 276E+01 


0 . 536E+01 


326E+04 


0. 276E+01 


-0 . 536E+01 


836E+01 


-0. 468E-02 


0 . 37 7E+02 


193E+04 


0. 241E+02 


0 . 202E+02 


378E+04 


-0.7 30E+01 


0. 781E+01 


378E+04 


0 . 730E+01 


-0. 781E+01 


1 93E+04 


-0 . 241E+02 


-0. 202E+02 



R6 


R5 


R4 


RE 


. 1 00E+01 


0. 100E+01 


0 . 668E+00 


0 . 663E+00 


. 100E+01 


0. 100E+01 


0. 384E+00 


0. 384E+00 


. 100E+01 


0 . 100E+01 


0.384E+00 


0. 384E+00 


. 100E+01 


0 . 100E+01 


0 . 66 8E+0 0 


0 . 663E+00 


. 100E+01 


0 . 100E+01 


0 . 735E+00 


0 . 661E+0 0 


. 100E+01 


0 . 100E+01 


0. 769E+00 


0. 286E+00 


. 100E+01 


0. 100E+01 


0. 769E+00 


0. 286E+00 


. 100E+01 


0. 100E+01 


0 . 735E+00 


0 . 661E+00. 


. 100E+01 


0. 100E+01 


0. 869E+00 


0. 533E+00 


. 100E+01 


0. 100E+01 


0 . 966E+00 


0 . 129E+00 


. 100E+01 


0. 100E+01 


0. 966E+00 


0 . 129E+00 


. 100E+01 


0. 100E+01 


0. 869E+00 


0. 533E+00 


. 100E+01 


0 . 100E+01 


0 . 1 00E+01 


0. 996E+00 


. 100E+01 


0. 100E+01 


0. 100E+01 


0. 999E+00 


. 100E+01 


0. 100E+01 


0. 100E+01 


0. 999E+00 


. 100E+01 


0. 100E+01 


0. 100E+01 


0. 996E+00 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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17 


0. 999E+00 


0 . 998E+00 


0. 100E+01 


18 


0. 100E+01 


0. 990E+00 


0. 100E+01 


19 


0 . 100E+01 


0. 990E+00 


0. 100E+01 


20 


0 . 999E+00 


0. 998E+00 


0. 100E+01 


21 


0 . 997E+00 


0. 100E+01 


0 . 100E+01 


22 


0 . 999E+00 


0. 100E+01 


0 . 100E+0 1 


23 


0. 999E+00 


0 . 100E+01 


0 . 100E+01 


24 


0 . 997E+00 


0. 100E+01 


0 . 100E+01 


25 


0 . 100E+01 


0 . 100E+01 


0.100E+01 


26 


0 . 100E+01 


0. 100E+01 


0 . 1 00E+01 


27 


0. 100E+01 


0 - 100E+01 


0 . 1 00E+01 


28 


0.100E+01 


O.OOOE+OO 


0. 1 00E+01 


29 


0. 100E+01 


0 . 100E+0 1 


0 . 100E+01 


30 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


31 


0. 100E+01 


0. 100E+01 


0. 100E+01 


32 


0. 100E+C1 


0. 100E+01 


0 . 100E+01 


TOTAL REL. 


RT1 


RT2 


RT3 




0. 960E+00 


O.OOOE+OO 


O.OOOE+OO 



TOTAL STRUCTURE RELIABILITY = O.OOOE+OO 



. 100E+01 


0 . 100E+01 


0 . 999E+00 


0.99 7E+00 


. 100E+01 


0. 100E+01 


0. 998E+00 


0 . 9 88E+00 


. 100E+0 1 


0. 100E+0 1 


0. 998E+00 


0.988E+00 


. 100E+01 


0 . 1 00E+01 


0 . 9 99E+00 


0 . 997E+00 


. 100E+0 1 


0 . 1 OOE+O 1 


0. 100E+01 


0 . 997E+00 


. 100E+01 


0. 100E+01 


0 . 99 7E+00 


0.996E+00 


. 100E+01 


0. 100E+01 


0 . 997E+00 


0 . 996E+00 


. 100E+01 


0 . 100E+01 


0.100E+01 


0. 997E+00 


. 100E+01 


0 . 1 00E+01 


0 . 100E+01 


0 . 100E+01 


. 1 OOE+O 1 


0 . 1 OOE+O 1 


0. 989E+00 


0.989E+00 


. 100E+01 


0 . 100E+01 


0 . 989E+00 


0 . 9 89E+00 


. 100E+01 


0. 100E+01 


0. 100E+01 


O.OOOE+OO 


. 1 OOE+O 1 


0 . 100E+01 


0 . 999E+00 


0 . 999E+00 


. 1 00E+01 


0 . 100E+01 


0. 985E+00 


0 . 985E+00 


. 100E+01 


0. 100E+01 


0. 985E+00 


0.985E+00 


. 100E+01 


0 . 100E+01 


0 . 999E+0 0 


0 . 999E+00 


RT6 


RT5 


RT4 


RT 


O.OOOE+OO 


O.OOOE+OO 


O.OOOE+OO 


0 . OOOE + OO 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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INPUT 



FILE 



a ADINA-IN 

A 

A THM3D4.IN EIGHT SYMMETRIC LAYER USING 3-D. 

A 

FILEUNITS LI ST=B L0G = 7 ECH0=7 
CONTROL PLOTUNIT=CM 
COLORS BCODE=BLUE 

A 

DATABASE CREATE 
WORKSTATION DEVICE=0 

A 

HEADING ' THM3D4 . IN ; aaaeIGHT LAYERS SOLUTION 90,75 DEG aaa' 

A 

MASTER IDOF=0001 1 1 

PRINTOUT VOLUME=MAXIMUM IPD=1 IPRIC=0 P=STRAINS 
PORTHOLE FORMATTED=YES FILE=60 
A 

COORDINATES / ENTRIES NODE X Y Z 

1 - 20 . - 20 . 0.1 
2 0 . - 20 . 0.1 

3 20. -20. 0.1 

4 -20. 0. 0.1 

5 0. 0.0.1 

6 20 . 0 . 0.1 

7 -20. 20. 0.1 

8 0 . 20 . 0.1 

9 20. 20. 0.1 

21 - 20 . - 20 . 0 . 

22 0 . - 20 . 0 . 

23 20. -20. 0. 

24 -20. 0. 0. 

25 0. 0. 0. 

26 20. 0. 0. 

27 -20. 20. 0. 

28 0. 20. 0. 

29 20. 20. 0. 

101 5. 0. 0. 

102 0 . 10 . 0 . 

103 5. 10. 0. 

104 15. 2.68 0. 

105 7.68 10. 0. 

106 15. 10. 0. 

107 15. -10. 0. 



ANOTE THE MAXES MEAN : 1=90DEG 2=15 3=75 4=45 5=-45 

A================:===================================== 

MATERIAL 1 ORTHOTROPIC 181. E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 

GAB= 7 . 1 7E0 5 7.17E05 7.17E05 

MATERIAL 2 ORTHOTROPIC 53.8E05 17.9E05 17.9E05 0.083 0.083 0.25, 

GAB=8 . 6E05 8.6E05 8.6E05 

A 

SYSTEM 1 CYLINDRIC X=0 Y = 0 Z = 0 PHI = 0 THETA= - 90 . XSI = 90. 



COORDINATES / ENTRIES NODES R THETA 



10 

11 

12 

13 

14 

15 



10 . 

10 . 

10 . 

10 . 

10 . 

10 . 



0 . 
45. 
90. 
135 . 
180 . 
225, 



XL 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 
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16 


10 . 


270. 


0.1 


17 


10. 


315. 


0.1 


18 


10. 


360. 


0.1 


30 


10. 


0 . 


0 . 


31 


10. 


45. 


0 . 


32 


10. 


90. 


0 . 


33 


10. 


135. 


0 . 


34 


10. 


180. 


0 . 


35 


10. 


225. 


0 . 


36 


10. 


270. 


0 . 


37 


10. 


315. 


0 . 


38 


10. 


360. 


0 . 



A 



LINE 


CYLINDRIC 


11 


10 


EL = 1 


MIDNODES = 1 


NC=A 


LINE 


CYLINDRIC 


12 


11 


EL=1 


MIDN0DES=1 


NC = A 


LINE 


CYLINDRIC 


13 


12 


EL=1 


MIDNODES= 1 


NC=A 


LINE 


CYLINDRIC 


14 


13 


EL = 1 


MIDN0DES=1 


NC = A 


LINE 


CYLINDRIC 


15 


14 


EL = 1 


MIDNODES = 1 


NC = A 


LINE 


CYLINDRIC 


16 


15 


EL= 1 


MIDNODES= 1 


NC=A 


LINE 


CYLINDRIC 


17 


16 


EL = 1 


MIDNODES =1 


NC=A 


LINE 


CYLINDRIC 


18 


17 


EL = 1 


MIDNODES =1 


NC = A 


LINE 


CYLINDRIC 


31 


30 


EL= 1 


MIDNODES= 1 


NC=A 


LINE 


CYLINDRIC 


32 


31 


EL= 1 


MIDN0DES=1 


NC = A 


LINE 


CYLINDRIC 


33 


32 


EL= 1 


MIDNODES= 1 


NC=A 


LINE 


CYLINDRIC 


34 


33 


EL= 1 


MIDN0DES=1 


NC =A 


LINE 


CYLINDRIC 


35 


34 


EL= 1 


MIDNODES = 1 


NC = A 


LINE 


CYLINDRIC 


36 


35 


EL = 1 


MIDNODES =1 


NC=A 


LINE 


CYLINDRIC 


37 


36 


EL= 1 


MIDNODES = 1 


NC = A 


LINE 


CYLINDRIC 


38 


37 


EL = 1 


MIDNODES = 1 


NC=A 



A 



EGROUP 1 THREEDSOLID MATERIALS TABLES 

GVOLUME 7 4 18 17 27 24 38 37 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 4 1 11 10 24 21 31 30 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 1 2 12 11 21 22 32 31 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 2 3 13 12 22 23 33 32 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 3 6 14 13 23 26 34 33 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 6 9 15 14 26 29 35 34 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 9 8 16 15 29 28 36 35 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 8 7 17 16 28 27 37 36 EL1=1 EL2=1 EL3=4 N0DES=27 

GSURFACE 7 8 28 27 EL1=1 EL2=4 NO=9 

GSURFACE 8 9 29 28 EL1=1 EL2=4 NO=9 

AXES-ORTHOTROPIC 1 N0DI=101 N0DJ=103 NODK=102 

AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 

AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 

. AXES-ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 

AXES-ORTHOTROPIC 5 N0DI=101 N0DJ=107 N0DK=102 

*STRES STABLE 112345678 21 26 27 
STRES STABLE 1 21 

EDATA / ENTRIES EL TABLE MAXES 







1 


1 


1 


STEP 


o 

H 


29 


1 


1 






2 


1 


3 


STEP 


4 TO 


30 


1 


3 






3 


1 


3 


STEP 


4 ,TO 


31 


1 


3 






4 


1 


1 


STEP 


o 

H 
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32 1 1 

A 



LOADS ELEMENT 

9 2 -3.0E04 

TO 

16 2 -3.0E04 



A 

FIXB 123 LINES / 8 28 

FIXB 23 SURFACES / 7 8 28 27 / 8 9 29 28 

A 

A 

EGROUP 2 THREEDSOLID MATERIAL =2 TABLES 



GVOLUME 


11 


12 


5 


5 


31 


32 


25 


25 


EL1 = 1 


EL2 = 1 


EL3 = 4 


N0DES=27 


GVOLUME 


12 


13 


5 


5 


32 


33 


25 


25 


EL1 = 1 


EL2 = 1 


EL3 = 4 


N0DES=27 


GVOLUME 


13 


14 


5 


5 


33 


34 


25 


25 


EL1 = 1 


EL2 = 1 


EL3 =4 


N0DES=27 


GVOLUME 


14 


15 


5 


5 


34 


35 


25 


25 


EL1 = 1 


EL2 = 1 


EL3 = 4 


NODES =27 


GVOLUME 


15 


16 


5 


5 


35 


36 


25 


25 


EL1 = 1 


EL2 = 1 


EL3 = 4 


NODES=27 


GVOLUME 


16 


17 


5 


5 


36 


37 


25 


25 


ELI = 1 


EL2 = 1 


EL3 = 4 


NODES=27 


GVOLUME 


17 


18 


5 


5 


37 


38 


25 


25 


EL1 = 1 


EL2 = 1 


EL3 = 4 


NODES =27 


GVOLUME 


10 


11 


5 


5 


30 


31 


25 


25 


EL1 = 1 


EL2 = 1 


EL3 = 4 


NODES =27 



AXES-ORTHOTROPIC 1 N0DI=101 N0DJ=103 N0DK=102 
AXES -ORTHOTROPIC 2 N0DI=101 NODJ= 1 04 ' NODK= 102 
AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 
AXES-ORTHOTROPIC 4 N0DI=101 NODJ=106 NODK=102 
AXES-ORTHOTROPIC 5 N0DI=101 NODJ=107 NODK=102 



ASTRESSTABLE 1 


1 2 


3 4 5 


6 7 8 


STRESSTABLE 1 


21 






EDATA / ENTRIES 


EL 


TABLE 


MAXES 




1 


1 


1 


STEP 4 TO 


29 


1 


1 




2 


1 


3 


STEP 4 TO 


30 


1 


3 




3 


1 


3 


STEP 4 TO 


31 


1 


3 




4 


1 


1 


STEP 4 TO 


32 


1 


1 



A 

A 



. EGZONE NAME= PARENT / 1 
EGZONE NAME= PATCH / 2 

FRAME HEAD I NG= UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH ZONE= PARENT NODES=0 ELEMENT = 1 HIDDEN= DASHED BCODE=ALL 
FRAME HEADING =UPPER XSF=1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH Z ONE = PATCH NODES=0 ELEMENT =1 HIDDEN =DASHED BCODE=ALL 
A 

ADINA 

A 

END 
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>nnnnnnnnnoonnnnnnno 



AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 



PROGRAM TRANS3 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

THIS PROGRAM TRANFORMS THE STRESSES OUTPUT FROM ADINA, 

EXPRESSED IN THE GLOBAL STRUCTURE COORDINATE SYSTEM, TO 
THE PRINCIPAL MATERIAL COORDINATE. AND THEN, CALCULATES THE 
RELIABILTY PER STRESS COMPONENT AND PER ELEMENT AND FINALLY 
EVALUATES THE TOTAL RELIABILITY. ALL THAT BY CALLING THE 
SUBROUTINE HEREIN. 

NEL : # OF WHOLE SRUCTURE ELEMENTS. 

NELP : ft OF PARENT STRUCTURE ELEMENTS 

SX,SY, . . . ,SYZ : THE STRESS IN GLOBAL COORDINATE SYSTEM. 

THETA : THE LAYUP ANGLE PER ELEMENT. 

S1,S2 S6 : THE STRESS IN THE MATERIAL PRINCIPAL COORDINATE. 

R1,R2 R6 : THE RELIABILITY PER STRESS COMPONENT. 

RT1 ,RT2 , . . . RT6 : THE TOTAL RELIABILTY OF STRUCTURE PER COPMPONENT. 
RE : THE RELIABILITY PER ELEMENT 

RT : THE TOTAL RELIABILITY. . 

ALF : THE SHAPE FUNCTION PARAMER FOR WEI BULL 

BEFORE RUNNING THE PROGRAM VERIFY IN THE PROGRAM THE ANGLE THETA 
INPUT AND THE # OF ELEMENT. AND ALSO VERIFY THE MATERIAL PROEPRTIES 
GIVEN IN SUBROUTINE BETA. 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 



INTEGER NEL 

PARAMETER ( NEL=64 ,NELP=32 ,PI = 3.1415927 ) 

PARAMETER ( ROWS = l 00 f COLS = 6 , C0LS1 = 7 ) 

INTEGER ITHETA(NEL) 

REAL A 8 A (ROWS, COLS) r B ( ROWS , COLS ) , C ( ROWS , COLS 1 ) 

REAL SX(NEL) , SY ( NEL) , SZ ( NEL) r SXY ( NEL ) ,SXZ(NEL) ,SYZ(NEL) 
REAL L,M, THETA ( NEL ) , THETAF ( NEL ) , ALF 
CHARACTER* 20 FNAME , FNOUT , FNAME1 , FNAME2 



WRITE (5,*) 'ENTER INPUT STRESS FN ?' 

READ (5, ' (A) ' ) FNAME 

WRITE (5,*) 'INPUT STRESS FN = FNAME 
OPEN (UNIT=15 ,FILE=FNAME, STATU S= ' OLD' ) 

WRITE (5,*) 'ENTER INPUT ELEM. ANGLE FN? ' 

READ (5, ' (A) ' ) FNAME2 

WRITE (5,*) 'INPUT ELEM. ANGLE FN = ' , FNAME 2 
OPEN ( UNIT=16 ,FILE=FNAME2 , STATUS = ' OLD' ) 

WRITE (5, *) 'ENTER AREAS OF ELEMENTS FN ?' 

READ ( 5 , ' ( A ) ' ) FNAME1 

OPEN ( UNIT=1 7, FILE=FNAME1, STATUS =' OLD' ) 

WRITE (5, A) ' ENTER OUTPUT FN ?' 

READ (5, ' (A) ' ) FNOUT 

OPEN ( UNIT= 13 , FILE=FNOUT , STATUS = ' NEW' ) 

WRITE (13, A) 'INPUT STRESS FN = ', FNAME 
WRITE (13, A) 'INPUT ELEM. ANGLE FN = ' ,FNAME2 
WRITE (13, A) 'INPUT AREAS OF ELEMENTS FN = ' ,FNAME1 
WRITE (13, A) ' OUTPUT FN =', FNOUT 

A , 

a INPUT OF ANGLE ( THETA ) FOR EA. ELEMENT ; 



DO 20 1=1, NEL 

READ (16, A) N,I2,I3,I4,I5,I6,ITHETA(I) , IB ,A9 ,A10 , 111 , 
& 112,113,114,115,116,117,118,119 

IF. (ITHETA(I) .EQ.l) THETA(I) = 90.0 
IF (ITHETA(I) .EQ.3) THETA(I) = 75.0 
WRITE (5,25) N , THETA ( I ) 

THETAF ( I ) = THETA ( I ) a PI/ 180.0 
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20 CONTINUE 

A 

WRITE (13, A) ' OUTPUT : ' ,FNOUT 

A 

A TRANSFORMATION OF STRESSES TO PRINCIPAL AXES STRESS ; 

A 

WRITE (13,11) 'NEL' , ' SX' , # SY' , ' SZ ' , ' SXY ' , ' SXZ ' , ' SYZ ' , ' ANGLE ' 

DO 55 1=1, NEL 

READ (15, A) I EL, NOD, SX (I) , SY ( I ) ,SZ(I) , SXY ( I ) , SXZ ( I ) , SYZ ( I ) 
THETA ( I ) = THETAF(I) A180. 0 / PI 

WRITE (13,61) IEL , SX ( I ) , SY ( I ) , SZ ( I ) , SXY ( I ) , SXZ ( I ) , SYZ ( I ) , 

& THETA ( I ) 

A 

L=COS (THETAF ( I ) ) 

M=SIN ( THETAF ( I ) ) 

A ( I , 1 ) = SX(I)A(LAA2) + SY(I)A(MAA2) + 2ASXY(I)alam 
A( I , 2 ) = SX(I)A(MAA2) + SY(I)A(LAA2) - 2ASXY(I)ALAM 
A ( 1 , 3 ) = SZ(I) 

A ( 1 , 6 ) = -SX(I)*L*M + SY(I)AMAL + SXY(I)A(LAA2 - MAA2) 

A( I , 5 ) = SYZ(I)AM + SXZ(I)AL 
A( 1 , 4 ) = SYZ(I)AL - SXZ(I)AM 
55 CONTINUE 

WRITE (13,10) ' NEL ' , ' S 1 ' , ' S 2 ' , ' S 3 ' , ' S 6 ' , ' S 5 ' , ' S 4 ' 

DO 80 1=1, NEL 

WRITE (13,60) I , A( I , 1 ) ,A ( I , 2 ) , A( I , 3 ) , A( I , 6 ) , A( I , 5 ) , A( I , 4 ) 

80 CONTINUE 

A 

a SUBROUTINES FOR CALCULATING THE RELIABILITY ; 

A 

CALL BETA( A,B,NEL) 

WRITE ( 5 , a ) 'INPUT ALFA' 

READ ( 5 , a ) ALF 

WRITE (13, A) 'ALFA : ' , ALF 

CALL ALFA ( A, B, NEL, ALF) 

CALL RELI ( A , B ,NEL ) 

CALL ELRE( B ,NEL) 

A 



10 


FORMAT 


( // ,1X,A10,6(A7,4X) 


/ 


IX , 80 ( ' 


) ) 


11 


FORMAT 


( // , IX , A10 , 7 ( A7 , 4X) 


/ 


IX , 80 ( 1 


) ) 


15 


FORMAT 


(IX, ' THETA ( ' ,12, ' ) = 


' 


) 




25 


FORMAT 


(IX, ' THETA ( ' ,12, ' ) = 


' 


,F6. 1) 




60 


FORMAT 


( / ,4X,I6,6E11. 3) 








61 


FORMAT 


( / ,4X,I6,6E11.3,3X, 


F6 


. 1) 





* A 

STOP 

END 



AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE BETA( A , B ,NEL ) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( ROWS =100 , COLS =6 ) 

REAL*8 A(ROWS, COLS) ,B(ROWS, COLS) 

REAL X(COLS) ,X1C,X2C,X3C 

A 

A PLY STRESS - STRENGTH, RATIO CALCULATION : 

A 

WRITE (13,10) ' NEL ' , 'SRI' , 'SR2' , 'SR3' , 'SR6' , 'SR5' , 'SR4' 
*A. FOR THE PARENT STRUCTURE, 

A 
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DO 90 1=1, HELP 

: STRENGTH PROPERTIES FOR GRAPHITE EPOXY, 



X(l) = 15.0E4 
X(2) = 0.4E4 

X ( 3 ) = X ( 2 ) 

X ( 6 ) = 0.68E4 

X ( 4 ) = X(6) 

X ( 5 ) = X(6) 

X1C = -15.0E4 
X2C = -2.46E4 
X3C = X2C 

IF ( A ( I , 1 ) . LT . 0 ) X(l) = X1C 
IF ( A ( I , 2 ) . LT . 0 ) X( 2 ) = X2C 
IF ( A ( I , 3 ) . LT . 0 ) X( 3 ) = X3C 
DO 80 J=l,COLS 

B(I,J) = ABS ( A ( I , J ) / X( J) ) 

80 CONTINUE 

WRITE (13,60) I , B ( I , 1 ) , B( I , 2 ) , B ( I , 3 ) , B( I , 6 ) , B ( I , 5 ) , B( I , 4 ) 
90 CONTINUE 

*B. FOR THE PATCH , 

A 

N=NELP+1 
DO 100 I=N,NEL 

C STRENGTH PROPERTIES FOR E-GLASS EPOXY, 

C 

X(l) = 10.34E4 
X( 2 ) = 0.276E4 

X ( 3 ) = X ( 2 ) 

X ( 6 ) = 0.414E4 

X ( 4 ) = X( 6 ) 

X ( 5 ) = X ( 6 ) 

X1C - -10.3E4 
X2C = -1.38E4 
X3C = X2C 

IF ( A ( I , 1 ) . LT . 0 ) X(l) = X1C 
IF ( A ( I , 2 ) . LT . 0 ) X ( 2 ) = X2C 
IF (Ad, 3) . LT . 0 ) X ( 3 ) = X3C 
DO 95 J= 1 , COLS 

B( I , J) = ABS ( A ( I , J ) / X ( J ) ) 

95 CONTINUE 

WRITE (13,60) I,B(I,1) ,B(I,2) , B ( I , 3 ) ,B(I,6) , B( I , 5 ) ,B(I,4) 
100 CONTINUE 

10 FORMAT ( // ,1X,A10,6(A7,4X) / 1X,B0('-')) 

60 FORMAT ( / , 4X , 16 , 6E1 1 . 3 ) 

RETURN 

END 

AAAAAAAAAAAAA'AAAAAAAAAAAAAAAAAAAA 

SUBROUTINE ALFA ( A , B , NEL , ALF ) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( ROWS= 100 , COLS=6 ) 

REALA8 A (ROWS, COLS) , B ( ROWS , COLS ) 

REAL AREA ( 100) , AMIN 

A 

* PLY STRESS - STRENGTH, RATIO RISED BY ALFA CALCULATION : 

A — — — — — — — — — — — — 

WRITE (13,10) ' NEL ' , ' VI ' , ' V2 ' , ' V3 ' , ' V6 ' , ' V5 ' , ' V4 ' 

DO 90 1=1, NEL 

READ (17,*) IIEL, AREA( I ) ,AMIN 
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AM I N = 4 0 0 . 

WRITE (5,95) IIEL, AREA( I ) , AMIN 
WRITE (3,95) I I EL , AREA( I ) , AMIN 
DO 80 J = 1 , COLS 

a A(I,J) = B(I,J)aaalf 

A ( I , J ) = ( B( I , J ) aaaLF) a AREA ( I ) /AMIN 
80 CONTINUE 

WRITE (13,60)I,A(I,1) ,A(I,2) ,A(I,3) ,A(I,6) ,A(I,5) , A ( 1 , 4 ) 

90 CONTINUE 

10 FORMAT ( // ,1X,A10,6(A7,4X) / 1X,80('-')) 

60 FORMAT ( / , 4X , 16 , 6E11 . 3 ) 

95 FORMAT ( / / IX, ' AREA( ',12 , ' ) = ',2F10.2) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE RELI ( A , B , NEL ) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( ROWS= 100 , C0LS=6 ) 

REAL*8 A(ROWS, COLS) ,B(ROWS, COLS) 

A 

* PLY STERSS TENSOR RELIABILITY CALCULATION (WEI BULL) : 

A : 

WRITE (13,10) ' NEL ' , ' R1 ' , ' R2 ' , ' R3 ' , ' R6 ' , ' R5 ' , ' R4 ' 

DO 90 1=1, NEL 
DO 80 J = 1 , COLS 

B ( I , J ) = EXP ( -A ( I , J ) ) 

80 CONTINUE 

WRITE (13,60) I , B ( I , 1 ) ,B(I,2) ,B(I,3) ,B(I,6) , B ( I , 5 ) ,B(I,4) 

90 CONTINUE 

10 FORMAT ( // ,lX r A10,6(A7,4X) / 1X,80('-')) 

60 FORMAT ( / , 4X , 16 , 6E1 1 . 3 ) 

RETURN 

END 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

SUBROUTINE ELRE ( B , NEL ) 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

PARAMETER ( ROWS = 100, COLS = 6 , COLS 1 = 7) 

REALA e B( ROWS, COLS) , C( ROWS , COLS1 ) , RT 

A 

AELEMENTS RELIABILITY CALCUL. , AS WELL AS "FIBER" (Rl) AND "MATRIX" 
ARELI ABILITY CALCUL. & TOTAL STRUCTURE RELIABILITY CALCUL. 

A 

WRITE (13,10) 'NEL' , ' R1 ' , ' R2 ' , ' R3 ' , ' R6 ' , ' R5 ' , 'R4 ' , ' RE ' 

RT =1.0 

C(NEL+1 , 1 ) =1 . 0 
DO 90 1=1, NEL 
C ( I , COLS1 ) = 1.0 
DO 80 J=1 ,COLS 

C( I , J ) =B( I , J ) 

C ( I , COLS1 ) = C ( I , COLS1 ) a C(I,J) 

C ( NEL+1 , J ) = C ( NEL+1 , J ) a C(1,J) 

80 CONTINUE 

RT = RT a C ( I , COLS1 ) 

WRITE (5,60) I,C( 1 ,1) ,C( 1 ,2) ,C( 1,3) ,C( 1 ,4 ) ,C( 1,5) ,C( 1 ,6) ,C( 1 ,7 ) 
WRITE (13,60) I,C(I,1) , C ( I , 2 ) ,C(I,3) ,C(I,4) ,C(I,5) ,C(I,6) , C ( I , 7 ) 
90 CONTINUE 

WRITE (13,10) ' TOTAL REL . : ' , ' RT1 ' , ' RT2 ' , ' RT3 ' , ' RT6 ' , ' RT5 ' , 

& ' RT4 ' , ' RT ' 

WRITE (13,65) C ( NEL+1 , 1 ) , C ( NEL+1 , 2 ) , C ( NEL+1 , 3 ) , C( NEL+1 , 4 ) , 
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10 

60 

65 

70 



& C(NEL+1,5) ,C(NEL+1,6) ,RT 

WRITE (5,70) RT 
WRITE (13,70) RT 

FORMAT ( // ,1X,A10,7(A7,4X) / 1X,95('-')) 
FORMAT ( / ,4X,I6,7E11.3) 

FORMAT ( / ,11X,7E11.3) 

FORMAT (//, IX, 'TOTAL STRUCTURE RELIABILITY 

RETURN 

END 



, El 1 . 



/ IX, 4 0 ( ' 
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INPUT STRESS FN = THM3D4 . STR 
INPUT EL EM. ANGLE FN = THM3D4.EL 
INPUT AREAS OF ELEMENTS FN =-THM3D4 . ARE 
OUTPUT FN =THM3D4 . REL 
OUTPUT :THM3D4 . REL 



NEL SX 


SY 


SZ 


SXY 


sxz 


SYZ 


ANGLE 


1 


- 0 . 277E+03 


0 .391E+05 


0 . 187E+02 


-0 . 235E+04 


0 . 909E+00 


0 . 622E+01 


90.0 


2 


0 . 279E+03 


0. 209E+05 


C . 160E+02 


0. 236E+04 


0. 849E+00 


'0 . 604E+01 


75.0 


3 


0 . 279E+03 


0.209E+05 


0. 160E+02 


0 . 236E+04 


-0. 849E+00 


-0 . 604E+01 


75.0 


4 


-0 . 277E+03 


0 . 391E+05 


0 . 187E+02 


-0 . 235E+04 


-0/909E+00 


-0.622E+01 


90.0 


5 


-0. 278E+03 


0.391E+05 


0 . 248E+01 


-0 . 23 5E+04 


0. 399E+00 


0.675E+01 


90.0 


6 


0 . 279E+03 


0. 209E+05 


0 . 217E+01 


0 . 235E+04 


0.361E+00 


0 . 683E+01 


75.0 


7 


0 . 279E-I-03 


0 . 209E+05 


0.217E+01 


0 . 2 35E+04 


-0 . 361E+00 


-0.683E+01 


75.0 


8 


-0.278E+03 


0 . 39 1E+05 


0 . 248E+01 


-0.235E+04 


-0 . 399E+00 


-0.675E+01 


90.0 


9 


- 0 . 278E+03 


0 . 391E+05 


0 . 653E-01 


-0.235E+04 


0 . 552E+01 


-0.2 16E+02 


90.0 


10 


0 . 2 78E+03 


0. 209E+05 


-0. 746E+00 


0 . 235E+04 


0 . 553E+01 


-0 . 216E+02 


75.0 


11 


0 . 278E+03 


0 . 209E+05 


-0. 746E+00 


0 . 235E+04 


-0 . 553E+01 


0 . 216E+02 


75.0 


12 


-0. 278E+03 


0 . 391E+05 


0 . 653E-01 


-0 . 235E+04 


-0. 552E+01 


0. 216E+02 


90.0 


13 


-0 . 278E+03 


0.391E+05 


0.11 3E+00 


-0 . 235E+04 


0.559E+01 


-0. 210E+02 


90.0 


14 


0 . 278E+03 


0 . 209E+05 


-0 . 692E+00 


0 . 235E+04 


0 . 559E+01 


-0.210E+02 


75.0 


15 


0 . 278E+C3 


0 . 209E+05 


-0.6S2E+00 


0. 235E+04 


- 0 . 559E+0 1 


0.210E+02 


75.0 


16 


-0 . 278E+03 


0 . 391E+05 


0.113E+00 


-0. 235E+04 


-0 . 559E+01 


0.210E+02 


90.0 


17 


-0.278E+03 


0. 391E+05 


0 . 250E+01 


-0 . 2 35E+04 


-0. 992E+00 


-0.434E+01 


90.0 


18 


0 . 279E+03 


0. 209E+05 


0.218E+01 


0 . 235E+04 


-0 . 967E+00 


-0. 438E+01 


75.0 


19 


0.279E+03 


0 . 209E+05 


0.21 8E+01 


0 . 235E+04 


0. 967E+00 


0.438E+01 


75.0 


20 


-0 . 278E+03 


0 . 391E+05 


0 . 250E+01 


-0.23 5E+04 


0 . 992E+0 0 


0 . 434E+01 


90.0 


21 


-0 . 278E+03 


0 . 391E+05 


0. 190E+02 


“0 . 2 36E+04 


-0 . 608E+00 


- 0 . 566E+01 


90.0 


22 


0 . 277E+03 


0. 209E+05 


0.163E+02 


0 . 2 35E+04 


-0. 550E+00 


-0.573E+01 


75.0 


23 


0.277E+03 


0 . 2 09E+05 


0 . 163E+02 


0 . 235E+04 


0 . 550E+00 


0 . 573E+01 


75.0 


24 


-0 . 278E+03 


0.39 1E+05 


0. 190E+02 


-0 . 236E+04 


0 . 608E+00 


0 . 566E+01 


90.0 


25 


-0 . 31 8E+03 


0 . 3 90E+05 


“0.11 5E+03 


-0 . 235E+04 


-0. 598E+01 


-0. 109E+01 


90.0 
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26 


0.243E+03 


0. 209E+05 


- 0 . 982E+02 


27 


0. 243E+03 


0.209E+05 


-0. 982E+02 


28 


-0.318E+03 


0.390E+05 


-0. 115E+03 


29 


-0. 320E+03 


0 . 390E+05 


-0 . 1 15E+ 03 


30 


0. 241E+03 


0.20 9E+05 


-0. 980E+02 


31 


0. 241E+03 


0. 209E+05 


-0 . 980E+02 


32 


-0. 320E+03 


0. 390E+05 


-0. 115E+03 


1 


-0. 281E+03 


0.391E+05 


-0. 594E+00 


2 


0. 275E+03 


0. 209E+05 


-0. 107E+01 


3 


0.27 5E+03 


0.209E+05 


-0. 107E+01 


4 


-0. 281E+03 


0. 391E+05 


-0 . 594E+0Q 


5 


-0 . 281E+03 


0 . 391E+05 


-0.646E+00 


6 


0. 275E+03 


0 . 209E+05 


-0. 117E+01 


7 


0 . 27 5E+03 


0 . 2 09E+05 


-0.117E+01 


8 


-0 . 281E+03 


0 . 391E+05 


-0 . 646E+0 0 


9 


-0.281E+03 


0. 391E+05 


0. 320E+00 


10 


0 . 275E+03 


0. 209E+05 


0. 274E+00 


11 


0 . 2 75E+03 


0 . 209E+05 


0. 274E+00 


12 


-0. 281E+03 


0 . 391E+05 


0. 320E+00 


13 


-0 . 282E+03 


0 . 391E+05 


0.641E+01 


14 


C. 274E+03 


0. 209E+C5 


0 . 552E+01 


15 


0. 274E+03 


0 . 209E+05 


0. 552E- L 01 


16 


-0.282E+03 


0. 391E+05 


0 . 64 1E+01 


17 


-0. 294E+03 


0. 391E+05 


-0. 422E+02 


18 


0. 263E+03 


0 . 209E+05 


-0. 361E+02 


19 


0 . 263E+03 


0 . 209E+05 


-0. 361E+02 


20 


- 0 . 2 94E+03 


0 . 391E+05 


-0 . 422E+02 


21 


-0. 295E+03 


0. 390E+05 


-0. 416E+02 


22 


0. 264E+G3 


0 . 209E+05 


-0. 356E+02 


23 


0. 264E+03 


0. 209E+05 


- 0 . 356E+02 



. 235E+04 


-0 . 607E+01 


-0. 414E+00 


75.0 


. 235E+04 


0. 607E+01 


0. 414E+00 


75.0 


. 235E+04 


0 . 598E+01 


0. 109E+01 


90.0 


. 2 3 5E+04 


-0. 634E+01 


-0. 996E+00 


90.0 


. 235E+04 


-0.621E+01 


-0.305E+00 


75.0 


. 235E+04 


0 . 621E+01 


0 . 305E+00 


75.0 


. 235E+04 


0.634 E+ 0 1 


0. 996E+00 


90.0 


. 235E+04 


0 . 2 0 5E+01 


-0. 802E+01 


90.0 


. 2 3 5E+04 


0.204E+01 


-0. 800E+01 


75.0 


. 235E+04 


-0. 204E+01 


0.800E+01 


75.0 


. 235E+04 


- 0 . 2 05E+01 


0 . 802E+01 


90.0 


. 235E+04 


0 . 206E+01 


-0. 778E+01 


90.0 


. 235E+04 


0 . 206E+01 


-0. 781E+01 


75.0 


. 235E+04 


-0 . 206E+01 


0 . 781E+01 


75.0 


. 2 35E+04 


-0.206E+01 


0. 778E+01 


90 . 0 


. 235E+04 


-0. 370E+00 


-0. 167E+01 


90.0 


. 235E+04 


-0. 359E+00 


-0. 170E+01 


75.0 


. 235E+04 


0.359E4-00 


0. 170E+01 


75.0 


. 2 3 5E+04 


0. 370E+00 


0. 167E+01 


90.0 


. 2 3 5E-+-04 


-0. 167E+00 


-0 . 2 11E+01 


90 . 0 


. 2 3 6E+04 


-0. 170E+00 


-0. 218E+01 


75.0 


. 236E+04 


0. 170E+00 


0. 218E+01 


75.0 


. 235E+04 


0. 167E+00 


0. 211E+01 


90.0 


. 2 36E+04 


-0. 238E+01 


-0. 128E+C1 


90.0 


. 235E+04 


-0 . 229E+01 


-0.4 81E+00 


75.0 


. 235E+04 


0 . 229E+01 


0 . 4 81E+00 


75.0 


. 236E+04 


0.238E+01 


0.128E+01 


90 . 0 


. 235E+04 


- 0 . 202E+01 


-0.119E+01 


90.0 


. 236E+04 


-0 . 220E+01 


-0.44 6E+00 


75.0 


. 236E+04 


0. 220E+01 


0 . 4 46E+00 


75.0 



0 

0 

-0 

-0 

0 

0 

-0 

-0 

0 

0 

-0 

-0 

0 

0 

-0 

-0 

0 

0 

-0 

-0 

0 

0 

-0 

-0 

0 

0 

-0 

-0 

0 

0 
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24 


-0.295E+03 


0 . 390E+05 


- 0 . 4 16E+02 


25 


-0 . 281E+03 


0. 391E+C5 


0 . 629E+01 


26 


0 . 275E+03 


0 . 209E+C5 


0 . 535E+01 


27 


0 . 275E+03 


0.209E+C5 


0 . 535E+01 


28 


-0 . 281E+03 


0 . 391E+05 


0 . 629E+01 


25 


-0.28CE+03 


0 . 39 1E+05 


0 . 293E+00 


30 


0 . 276E+03 


0 . 209E+05 


0.287E+00 


31 


0 . 276E+03 


0.209E+05 


0. 287E+00 


32 


-0.280E+C3 


0.391E^05 


0. 293E+00 


NEL 


SI 


S2 


S3 


i 


0. 2S1E+05 


-0 . 277E+03 


C.1E7E+02 


2 


0 . 207E+05 


0 . 4 85E X 03 


0 . 160E 4 *02 


3 


C . 207E+05 


0 . 485E+03 


0. 160E+02 


4 


0 . 391E+05 


-0.277E+03 


0.157E+02 


5 


0.391E^05 


-0 . 278E+03 


0.248E+01 


6 


C. 207E+05 


0. 466E+03 


0. 217E+01 


7 


C . 20 7E J -C5 


0. 486E+03 


0 . 217E+01 


6 


0 . 3 9 IE* 05 


-0 . 278E+C3 


C.248E-01 


5 


C.391E+05 


-0. 278E+03 


0 . 653E-01 


10 


0. 207E*C5 


0.485E^03 


-0. 746E+00 


11 


0 . 2 07E+05 


0 . 4 85E+C3 


-C . 746E+00 


12 


C . 3 51E+G5 


-0. 276E+03 


0. 653E-01 


13 


0. 391E-05 


-0. 278E-03 


0 . 113E+00 


14 


C. 207E-"05 


G. 4S5E+03 


-0 . 692E+00 


15 


0.207E-05 


0 . 4 85E*03 


-0 . 692E+00 


16 


0. 391E+05 


-0. 276E+03 


0. 113E+00 


17 


0. 391E+05 


-0 . 276E+03 


0. 250E+01 


ie 


0 . 207E+05 


0.485E+03 


0 . 218E+01 


15 


0 . 207E+C5 


C.485E+03 


0 . 218E+01 



. 235E+04 


0 . 202E+01 


0.119E+01 


90.0 


. 2 35E+04 


0 . 323E+00 


0.17 9E+01 


90.0 


. 235E+04 


0 . 29 3E+00 


0.211E+01 


75.0 


. 235E+04 


-0 . 293E+00 


-0.211E+01 


75.0 


. 235E+04 


-0.323E+00 


“0 . 179E+01 


90.0 


. 235E+04 


0 . 13 OE+OO 


0 . 284E+01 


90.0 


. 235E+04 


0. 140E+00 


0.243E+01 


75.0 


.235E+04 


*0. 140E+00 


-0.24 3E+01 


75.0 


. 235E+04 


-0. 130E+00 


-0 . 284E+01 


90.0 


S6 


S5 


S4 




. 235E+04 


0 . 622E+01 


-0.909E+00 




.313E-04 


0 . 605E+01 


0.743E+00 




. 313E+04 


-0 . 605E+01 


-0.74 3E+00 




. 235E+04 


-0 . 622E+01 


0. 909E+00 




. 235E+04 


0 . 675E+G1 


-0. 399E+00 




. 313E+04 


0. 669E+01 


0.14 2E+01 




. 31 3E+04 


-0 . 66SE+01 


-0 . 14 2E+01 




. 235E+04 


-0. 675E+01 


0 . 399E+00 




. 2 35E+04 


- 0 . 216E+02 


-0 . 552E+01 




. 3 1 3E+04 


-0. 194E+02 


“0. 105E+02 




. 3 1 3E+04 


0. 194E+02 


0. 105E+02 




. 235E+04 


0 . 216E+02 


0. 552E+01 




. 235E+04 


“0 . 21 0E+02 


-G . 559E+C1 




. 313E+04 


-0.1 89E+02 


“0 . 108E+C2 




. 313E+04 


0 . 1 89E+02 


0 . 1 08E+C2 




. 235E+04 


0. 210E+02 


0 . 559E+01 




. 235E+04 


“ 0 . 4 34E+01 


0. 992E+00 




. 313E+04 


-0.44 8E+01 


-0 . 198E+00 




. 313E+04 


0. 448E+01 


0. 198E+00 





0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 



0 . 391E+05 
0. 391E+05 
0. 207E+05 
0 . 207E+05 
0 . 391E+05 
0 . 390E+05 
0 . 207E+0 5 
0 . 207E+05 
0. 390E+05 
0.39CE+05 
0. 207E+05 
0. 207E+05 
0 . 39 0E+05 
0. 391E+05 
0. 207E+05 
0 . 207E+05 
0 . 391E+05 
0. 381E+05 
0 . 207E+05 
0. 2C7E+05 
0. 391E-C5 
0. 391E+05 
0. 207E+05 
0 . 207E+05 
0. 391E+05 
0 . 391E X 05 
0 . 207E+05 
0. 207E+05 
0. 391E+05 
0. 391E+05 



-0.278E+03 
-0 . 278E+03 
0 . 486E+03 
0 . 486E+03 
-0.278E+03 
-0.318E+03 
0.450E+03 
0.450E+03 
-0 . 318E+03 
-0. 320E+03 
0 . 449E+03 
0. 44 9E+03 
-0.320E+03 
-0. 2eiE+03 
0. 482E+03 
0 . 4 82E+03 
-0 . 281E+03 
- 0 . 261E+03 
C.462E+03 
0 . 4 62E+03 
- 0 . 281E+C3 
- 0 . 281E+C3 
0 . 4 82E+03 
0. 482E+03 
-0 . 281E+03 
-0 . 282E+03 
0.481E+03 
0. 481E+03 
- 0 . 282E+03 
-0 . 294E+03 



0. 250E+01 
0. 190E+02 
0 . 163E+02 
0 . 163E+02 
0 . 19 0E+02 
-0.115E+03 
-0 . 982E+02 
-0. 982E+02 
-0.1 15E+03 
-0. 115E+03 
-0. 980E+02 
-0. 980E-1 02 
-0. 115E+03 
-0 . 594E+00 
-0. 107E+01 
-0. 107E+01 
-0. 594E+00 
-0 . 646E+00 
-0. 117E+01 
-0. 117E+01 
-0 . 646E+00 
0. 320E+00 
0. 274E+00 
0 . 274E+0 0 
0. 320E+00 
0 . 641E+01 
0. 552E+01 
0 . 552E+0 1 
0 . 641E+01 
- 0 . 422E+02 



0 . 235E+04 
0 . 2 36E+04 
0.31 3E+04 
0 . 313E+04 
0.236E+04 
0 . 235E+04 
0.313E+04 
0. 313E+04 
0 . 235E+04 
0.235E+04 
0 . 313E+04 
0 . 3 13E+04 
0 . 235E+04 
0 . 235E+04 
0 . 313E+04 
0. 313E+04 
0.23 5E+04 
0 . 235E+04 
0.313E+04 
0. 313E+04 
0 . 235E+04 
0. 235E+04 
0. 313E+04 
0. 313E+04 
0.235E+04 
0. 235E+04 
0 . 3 13E+04 
0. 313E+04 
0 . 235E+04 
0 . 235E+04 



0. 434E+01 
-0 . 566E+0 1 
-0. 568E+01 
0. 568E+01 
0. 566E+01 
-0. 109E+01 
-0 . 197E+01 
0 . 197E+01 
0. 109E+01 
-0. 996E+00 
-0 . 1 90E+01 
0. 190E+01 
0 . 996E+00 
-0. 802E+01 
-0. 720E+01 
0 . 720E+01 
0 . 802E+0 1 
-0.778E+01 
-0.7 0 1E+01 
0 . 7 0 1E+0 1 
0. 778E+01 
-0. 167E+01 
-0. 173E+01 
0. 173E+01 
0. 167E+01 
-0. 211E+01 
-0.215E+01 
0 . 215E+01 
0.21 1E+01 
-0. 128E+01 



-0.992E+00 
0 . 608E+0 0 
-0 . 952E+00 
0. 952E+00 
-0 . 608E+0 0 
0. 598E+01 
0 . 575E+01 
-0.575E+01 
-0. 598E+01 
0.634E+01 
0 . 592E+0 1 
-0. 592E+01 
-0. 634E+01 
- 0 . 2 05E+01 
- 0 . 4 04E+0 1 
0. 404E+01 
0.205E+01 
-0.206E+01 
-0.401E+01 
0 . 4 0 1E+01 
0 . 206E+C1 
0.370E+00 
-0. 927E-01 
0.S27E-01 
-0. 370E+00 
0. 167E+00 
-0 . 4 01E+00 
0. 401E+00 
-0. 167E+00 
0 . 238E+0 1 
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50 0 . 207E+05 

51 0 . 207E+05 

52 0.391E+05 

53 0 . 39 0E+05 

54 0 . 207E+05 

55 0 . 207E+05 

56 0 . 390E+05 

57 0 . 391E+05 

58 0 . 207E+05 

59 0 . 207E+05 

60 0 . 39 1E+05 

61 0 . 391E+05 

62 0.207E+05 

63 0 . 207E+05 

64 0.391E+05 

ALFA : 5.000000 

NEL R1 

1 0 . 997E+00 

2 0.100E+01 

3 C.100E+01 

4 0 . 997E+00 

5 0 . 998E+0 0 

6 0.100E+C1 

7 0 . 1 00E+01 

• 6 0.998E+00 

9 0 . 997E+0 0 

10 0 . 100E+01 

11 0 . 100E+01 

12 0 . 997E+00 



0 . 4 7 1E+0 3 


-0 . 361E+02 


0. 471E+03 


-0 . 361E+02 


0 . 294E+0 3 


-0 . 422E+02 


0. 295E+03 


-0 . 416E+02 


0. 470E+03 


-0. 356E+02 


0.470E+03 


-0.356E+02 


0 . 295E+03 


-0 . 416E+02 


0 . 281E+03 


0.629E+01 


0 . 482E+03 


0. 535E+01 


0 . 4 82E+03 


0 . 535E+0 1 


0. 281E+03 


0 . 629E+0 1 


0 . 280E+03 


0.293E+00 


0 . 4-83E+03 


0. 287E+00 


0. 463E+03 


0 . 287E+00 


0. 280E+03 


0.293E+00 



P2 


R3 


0. 100E+01 


0. 100E+01 


0. lOOE-t-Ol 


0.100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0.100E+01 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 



0 . 313E+04 


-0. 106E+01 


0 . 313E+04 


0. 106E+01 


0 . 235E+04 


0 . 12BE+01 


0 . 235E+04 


-0 . 11 9E+01 


0 . 313E+04 


-0.10 0E+01 


0 . 313E+04 


0.100E+01 


0 . 235E+04 


0.119E+01 


0. 235E+04 


0. 179E+01 


0. 313E+04 


0. 212E+01 


0. 313E+04 


-0 . 212E+01 


0 . 235E+04 


-0. 179E+01 


0 . 23 5E+04 


0 . 284E+01 


0.313E+04 


0 . 238E+01 


0. 313E+04 


-0 . 238E-+-01 


0 . 235E+04 


-0. 284E+01 



R6 


R5 


0.100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0.100E+01 


0. 100E+01 


0 . 100E+C1 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 1 00E+01 


0.100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


0 . 1 00E+01 



0. 208E+01 
'0 . 208E+01 
-0 . 238E+01 
0 . 202E+0 1 
0 . 201E+01 
-0 . 201E+01 
-0. 202E+01 
-0. 323E+00 
0 . 264E+00 
-0 . 264E+00 
0. 323E+00 
-0. 130E+00 
0.493E+00 
-0 . 493E+00 
0.130E+00 



R4 


RE 


0.974E+00 


0. 971E+00 


0. 898E+00 


0.898E+00 


0 . 898E+00 


0. 898E+00 


0. 974E+00 


0. 971E+00 


0. 983E+00 


O’. 981E+00 


0. 931E+00 


0. 931E+00 


0 . 931E+00 


0. 931E+00 


0 . 983E+00 


0. 981E+00 


0. 974E+00 


0 . 971E+00 


0. 898E+00 


0. 898E+00 


0 . 898E+00 


0. 898E+00 


0 . 9 74E+00 


0. 971E+00 
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13 


0 . 998E+00 


0. 100E+01 


0.100E+01 


0. 100E+01 


0. 100E+01 


0. 983E+00 


0 . 981E+00 


14 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


0. 931E+00 


0 . 931E+00 


15 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 931E+00 


0. 931E+00 


16 


0. 998E+00 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 9B3E+00 


0 . 981E+00 


17 


0. 997E+00 


0. 100E+01 


0 . 100E+01 


0 . 100E+01 


0. 100E+01 


0 . 974E+00 


0 . 971E+00 


18 


0 . 100E+01 


0. 100E+01 


0 . 100E+01 


0.100E+01 


0 . 1 OOE+Ol 


0. 898E+00 


0 . 898E+00 


19 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


0. 100E+01 


0. 898E+00 


0 . 89BE+00 


20 


0. 997E+00 


0. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


0. 100E+01 


0. 974E+00 


0 . 971E+00 


21 


0 . 998E+00 


0 . 1 00E+01 


0 . 100E+0 1 


0.100E+01 


0 . 1 OOE+Ol 


0 . 9 83E+00 


0. 981E+00 


22 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


0 . 931E+00 


0 . 931E+00 


23 


0. 100E+01 


0 . 1 00E+01 


0. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


0. 931E+00 


0 . 931E+00 


24 


0. 998E+00 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 983E+00 


0 . 981E+00 


25 


0. 997E+00 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


0. 100E+01 


0 . 9 74E+00 


0. 971E+00 


26 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 898E+00 


0 . 898E+00 


27 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 1 00E+01 


0. 100E+01 


0 . 898E+00 


0.898E+00 


28 


0. 997E+00 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 1 OOE+Ol 


0. 974E+00 


0. 971E+00 


29 


0. 998E+00 


0. 100E+01 


0.100E+01 


0. 100E+01 


0. 100E+01 


0. 983E+00 


0 . 981E+00 


30 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 1 OOE+O 1 


0. 100E+01 


0 . 931E+00 


0. 931E+00 


31 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 931E+00 


0. 931E+00 


32 


0. 998E+00 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 983E+C0 


0. 981E+00 


33 


0. 999E+00 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0.99 5E+00 


0. 994E+00 


34 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 1O0E+O1 


0. 979E+00 


0. 978E+00 


35 


0. 100E+01 


0. 100E+01 


0 . 1 00E+01 


0. 100E+01 


0. 100E+01 


0. 979E+00 


0. 978E+00 


36 


0. 999E+00 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 995E+00 


0. 994E+00 


37 


0. 999E+00 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 1 OOE+Ol 


0 . 995E+00 


0. 994E+00 


38 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 979E+00 


0. 978E+00 


39 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 979E+00 


0. 978E+00 


40 


0 . 999E+00 


0. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


0 . 100E+01 


0 . 995E+00 


0. 994E+00 


41 


0 . 999E+00 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 995E+00 


0. 994E+00 


42 


0. 100E+01 


0. 100E+01 


0. 100E+01 


0 . 1 OOE+Ol 


0. 100E+01 


0 . 978E+00 


0. 978E+00 
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43 


0. 100E+01 


0. 100E+01 


0. 100E+01 


44 


0. 999E+00 


0. 100E+01 


0. 100E+01 


45 


0. 999E+00 


0. 100E+01 


0 . 100E+01 


46 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


47 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


48 


0. 999E+00 


0. 100E+01 


0 . 100E+01 


49 


0. 999E+00 


0. 100E+01 


0. 100E+01 


50 


0. 100E+01 


0. 100E+01 


0 . 100E+01 


51 


0. 100E+01 


0. 100E+01 


0. 100E+01 


52 


0 . 999E+00 


0. 100E+01 


0. 100E+01 


53 


0 . 999E+00 


0 . 100E+01 


0. 100E+01 


54 


0. 100E+01 


0. 100E+01 


0. 100E+01 


55 


0. 100E+01 


0. 100E+01 


0. 100E+01 


56 


0.99SE+00 


0. 100E+01 


0. 100E+01 


57 


0. 999E+00 


0. 100E+01 


0. 100E+01 


58 


0. 100E+C1 


0. 100E+01 


0 . 100E+01 


59 


0 . lOOE+Oi 


0. 100E+01 


0. 100E+01 


60 


0. 999E+00 


0. 100E+01 


0. 100E+01 


61 


0 . 999E+00 


0. 100E+01 


0. 100E+01 


62 


0. 100E+01 


0. 100E+01 


0. 100E+01 


63 


0 . 1 00E+01 


0. 100E+01 


0 . 1 00E+01 


64 


0. 999E+00 


0. 100E+01 


0. 100E+01 


TOTAL REL. 


RT1 


RT2 


RT3 




0. 944E+00 


0. OOCE+OO 


0. OOOE+OO 



TOTAL STRUCTURE RELIABILITY = 0.103E+00 



. 100E+01 


0. 100E+01 


0.97 BE+OO 


0 . 978E+00 


. 100E+01 


0. 100E+01 


0 . 995E+00 


0. 994E+00 


. 100E+01 


0. 100E+01 


0 . 995E+00 


0. 994E+00 


. 100E+01 


0 . 100E+01 


0.97 8E+00 


0. 978E+00 


. 100E+01 


0 . 1 00E+01 


0 . 978E+00 


0.97 8E+00 


. 100E+01 


0. 100E+01 


0 . 995E+00 


0.994E+00 


. 1 00E+01 


0. 100E+01 


0 . 995E+00 


0 . 994E+00 


. 100E+01 


0. 100E+01 


0.97 8E+00 


0. 978E+00 


. 100E+01 


0. 100E+01 


0. 978E+00 


0. 978E+00 


. 100E+01 


0. 100E+01 


0 . 995E+00 


0. 994E+00 


. 100E+01 


0. 100E+01 


0 . 995E+00 


0. 994E+00 


. 100E+01 


0. 100E+01 


0 . 979E+00 


0 . 979E+00 


. 100E+01 


0. 100E+01 


0 . 9 79E+00 


0.979E+00 


. 100E+01 


0 . 100E+01 


0 . 995E+00 


0. 994E+00 


. 100E+01 


0. 100E+01 


0. 995E+00 


0. 994E+00 


. 100E+01 


0. 100E+01 


0. 978E+00 


0. 978E+00 


. 100E+01 


0. 100E+01 


0.978E+00 


0. 978E+00 


. 100E+01 


0 . 100E+01 


0. 995E+00 


0. 994E+00 


. 100E+01 


0. 100E+01 


0 . 995E+00 


0. 994E+00 


. 100E+01 


0. 100E+01 


0.979E+00 


0 . 978E+00 


. 100E+01 


0. 100E+01 


0. 979E+00 


0. 978E+00 


. 100E+01 


0. 100E+01 


0 . 995E+00 


0. 994E+00 


RT6 


PT5 


RT4 


RT 


0. OOOE+OO 


0. OOOE+OO 


0. OOOE+OO 


0. 103E+00 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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ADINA-IN 



INPUT 



FILE 



A 
A 

A THA3D4 . IN , EIGHT SYMMETRIC LAYERS , REPAIR W/ DIFFER. LAYUP ANGLE. 

A 

FILEUNITS LIST=8 L0G = 7 ECHO =7 
CONTROL PLOTUNIT=CM 
COLORS BCODE=BLUE 

A 

DATABASE CREATE 
WORKSTATION DEVICE=0 

A 

HEADING ' THA3D4 . IN; ***EIGHT LAYERS , DIFFERENT PATCH ANGLE, £90,9011 ***' 
A 

MASTER IDOF=OOOH1 



PRINTOUT VOLUME= MAXI MUM IPD= 


1 IPRIC= 


o 

11 


STRAINS 


PORTHOLE FORMATTED=YES FILE= 


60 






A 








COORDINATES / ENTRIES NODE 


X 


Y 


Z 


1 


-20. 


-20. 


0.1 


2 


0. 


-20. 


0.1 


3 


20. 


-20. 


0.1 


4 


-20. 


0. 


0.1 


5 


0. 


0. 


0.1 


6 


20. 


0. 


0.1 


7 


-20 . 


20. 


0.1 


8 


0. 


20. 


0 . 1 


9 


20. 


20. 


0.1 


21 


-20. 


-20. 


0. 


22 


0. 


-20. 


0. 


23 


20. 


-20. 


0. 


24 


-20. 


0. 


0 . 


25 


0 . 


0. 


0. 


26 


20 . 


0. 


0. 


27 


-2 0. 


20, 


. 0. 


28 


0. 


20, 


. 0. 


29 


20. 


20 . 


. 0 . 


101 


5. 


0, 


. 0. 


102 


0. 


10, 


. 0. 


103 


5. 


10 


. 0. 


104 


15. 


2 


.68 0. 


105 


7.68 


10 


0. 


106 


15. 


10 


0. 


107 


15. 


-10 


0. 



A 



*N0TE THE MAXES MEAN : 1=90DEG 2=15 3=75 4=45 5 = -45 

A= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =r= = = =r = = =r = = =r=r = = = = = = = 

MATERIAL 1 ORTHOTROPIC 181. E05 10.3E05 10.3E05 0.0159 0.0159 0.28, 

GAB=7 . 17E05 7.17E05 7.17E05 

A 

SYSTEM 1 CYLINDRIC X=0 Y=0 Z=0 PHI=0 THETA=-90 . XSI=90. 



COORDINATES / ENTRIES NODES 


R THETA 


XL 


10 


10 


. 0 . 


0.1 


11 


10 


. 45. 


0.1 


12 


10 


. 90. 


0.1 


13 


10 


. 135. 


0.1 


14 


10 


. 180 . 


0.1 


15 


10 


. 225. 


0.1 


16 


10 


. 270. 


0.1 


17 


10 


. 315. 


0.1 
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18 


10 *. 


360. 


0.1 


30 


10. 


0 . 


0 . 


31 


10. 


45. 


0 . 


32 


10. 


90. 


0 . 


33 


10. 


135. 


0 . 


34 


10. 


180. 


0 . 


35 


10. 


225. 


0 . 


36 


10. 


270. 


0 . 


37 


10. 


315. 


0 . 


38 


10. 


360. 


0 . 



A 



LINE 


CYLINDRIC 


11 


10 


EL = 1 


MIDNODES =1 


NC=A 


LINE 


CYLINDRIC 


12 


11 


EL = 1 


MIDNODES= 1 


NC = A 


LINE 


CYLINDRIC 


13 


12 


EL=1 


MIDN0DES=1 


NC=A 


LINE 


CYLINDRIC 


14 


13 


EL = 1 


MIDN0DES=1 


NC =A 


LINE 


CYLINDRIC 


15 


14 


EL=1 


MIDNODES = 1 


NC=A 


LINE 


CYLINDRIC 


16 


15 


EL= 1 


MIDNODES =1 


NC=A 


LINE 


CYLINDRIC 


17 


16 


EL = 1 


MIDN0DES=1 


NC = A 


LINE 


CYLINDRIC 


18 


17 


EL = 1 


MIDN0DES=1 


NC=A 


LINE 


CYLINDRIC 


31 


30 


EL = 1 


MIDNODES =1 


NC=A 


LINE 


CYLINDRIC 


32 


31 


EL = 1 


MIDNODES = 1 


NC=A 


LINE 


CYLINDRIC 


33 


32 


EL = 1 


MIDNODES =1 


NC=A 


LINE 


CYLINDRIC 


34 


33 


EL = 1 


MIDNODES= 1 


NC = A 


LINE 


CYLINDRIC 


35 


34 


EL=1 


MIDNODES =1 


NC=A 


LINE 


CYLINDRIC 


36 


35 


EL = 1 


MIDNODES= 1 


NC = A 


LINE 


CYLINDRIC 


37 


36 


EL= 1 


MIDN0DES=1 


NC = A 


LINE 


CYLINDRIC 


38 


37 


EL = 1 


MIDNODES = 1' 


NC=A 



A 



EGROUP 1 THREEDSOLID MATERIALS TABLES 

GVOLUME 7 4 18 17 27 24 38 37 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 4 1 11 10 24 21 31 30 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 1 2 12 11 21 22 32 31 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 2 3 13 12 22 23 33 32 EL1 = 1 EL2 = 1 EL3 = 4 NODES = 27 

GVOLUME 3 6 14 13 23 26 34 33 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 6 9 15 14 26 29 35 34 EL1=1 EL2=1 EL3 = 4 NODES=27 

GVOLUME 9 8 16 15 29 28 36 35 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 8 7 17 16 28 27 37 36 EL1=1 EL2=1 EL3=4 NODES=27 

GSURFACE 7 8 28 27 EL1=1 EL2=4 NO=9 

GSURFACE 8 9 29 28 EL1=1 EL2=4 NO=9 

AXES-ORTHOTROPIC 1 NODI=101 N0DJ=103 NODK=102 

AXES-ORTHOTROPIC 2 NODI=101 NODJ=104 NODK=102 

AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 N0DK=102 

AXES-ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 

AXES-ORTHOTROPIC 5 NODI=101 NODJ=107 NODK=102 

ASTRESSTABLE 112345678 21 26 27 
STRESSTABLE 1 21 

EDATA / ENTRIES EL TABLE MAXES 









1 


1 


1 


STEP 


4 


TO 


29 


1 


1 








2 


1 


3 


STEP 


4 


TO 


30 


1 


3 








3 


1 


3 


STEP 


4 


TO 


31 


1 


3 








4 


1 


1 


STEP 


4 


TO 


32 


1 


1 



A 
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LOADS ELEMENT 

9 2 -3.0E04 

TO 

16 2 -3.0E04 

A 

FIXB 123 LINES / 8 28 

FIXB 23 SURFACES / 7 8 28 27 / 8 9 29 28 

A 

A 

EGROUP 2 THREEDSOLID MATER I AL = 1 TABLES 

GVOLUME 11 12 5 5 31 32 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 12 13 5 5 32 33 25 25 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 13 14 5 5 33 34 25 25 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 14 15 5 5 34 35 25 25 EL1=1 EL2=1 EL3=4 N0DES=27 

GVOLUME 15 16 5 5 35 36 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 16 17 5 5 36 37 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 17 18 5 5 37 30 25 25 EL1=1 EL2=1 EL3=4 NODES=27 

GVOLUME 10 11 5 5 30 31 25 25 EL1=1 EL2=1 EL3=4 N0DES=27 

AXES -ORTHOTROPIC 1 N0DI=101 NODJ=103 NODK=102 

AXES -ORTHOTROPIC 2 NODI=101 NODJ=104 N0DK=102 

AXES-ORTHOTROPIC 3 NODI=101 NODJ=105 NODK=102 

AXES -ORTHOTROPIC 4 NODI=101 NODJ=106 NODK=102 

AXES-ORTHOTROPIC 5 N0DI=101 NODJ=107 NODK=102 

ASTRESSTABLE 112345678 21 26 27 
STRESSTABLE 1 21 

EDATA / ENTRIES EL TABLE MAXES 









1 


1 


1 


STEP 


4 


TO 


29 


1 


1 








2 


1 


1 


STEP 


4 


TO 


30 


1 


1 








3 


1 


1 


STEP 


4 


TO 


31 


1 


1 








4 


1 


1 


STEP 


4 


TO 


32 


1 


1 



A 

A 

EGZONE NAKE=THREEDSOL ID / 1/2 

FRAME' HEADING=UPPER XSF = 1.0 YSF=0.5 XFMAX=22.5 YFMAX=17.0 
MESH ZONE=THREEDSOLID NODES =0 ELEMENT =1 H I DDEN= DASHED BCODE=ALL 

• * 

ADINA 

A 

END 
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INPUT STRESS FN = THA3D4 . STR 
INPUT ELEM. ANGLE FN = THA3D4 . EL 
INPUT AREAS OF ELEMENTS FN =THA3D4 .ARE 
OUTPUT FN =THA3D4 . REL 
OUTPUT -.THA3D4.REL 



NEL SX 


SY 


SZ 


SXY 


sxz 


SYZ 


ANGLE 


1 


-0 . 3 36E+02 


0 . 371E+05 


-0 . 470E+00 


-0.197E+04 


0 . 142E+01 


0 . 316E+01 


90.0 


2 


0 . 690E+03 


0. 216E+05 


-0. 521E+00 


0. 295E+04 


0. 133E+01 


0. 322E+01 


75.0 


3 


0 . 69 0E+03 


0 . 216E+05 


-0. 521E+00 


0.295E+04 


-0.133E+01 


-0. 322E+01 


75.0 


4 


-0. 336E+02 


0 . 371E+05 


-0. 470E+00 


-0.19 7E+04 


-0. 142E+01 


-0.316E+01 


90.0 


5 


-0.722E+03 


0 . 382E+05 


0. 137E+02 


-0. 203E+04 


-0. 123E+01 


0 . 284E+01 


90.0 


6 


-0.759E+02 


0 . 2 15E+05 


0. 158E+02 


0.27 9E+04 


-0.666E+00 


0. 504E+01 


75.0 


7 


- 0 . 759E+02 


0. 215E+05 


0. 15BE+02 


0. 279E+04 


0 . 666E+00 


-0. 504E+01 


75.0 


8 


-0. 722E+03 


0. 382E+05 


0. 137E+02 


-0. 203E+04 


0. 123E+01 


-0.284E+01 


90.0 


9 


-0 . 530E+03 


0. 393E+05 


0. 120E+02 


-0. 228E+04 


0 . 801E+00 


-0. 160E+02 


90.0 


10 


0.468E+02 


0 . 213E+05 


0. 137E+02 


0 . 248E+04 


0 . 433E+01 


-0 . 165E+02 


75.0 


11 


0 . 468E+02 


0 . 2 13E+05 


0 . 1 37E+02 


0. 248E+04 


-0 . 433E+01 


0. 165E+02 


75.0 


12 


-0. 530E+03 


0 . 393E+05 


0 . 120E+02 


-0 . 228E+04 


-0. 801E+00 


0 . 160E+02 


90.0 


13 


-0 . 189E+03 


0. 397E+05 


-0. 200E+02 


-0. 241E+04 


0. 390E+01 


- 0 . 153E+02 


90.0 


14 


0. 380E+03 


0. 213E+05 


-0. 166E+02 


0 . 238E+04 


0 . 4 07E+01 


- 0 . 1 66E+02 


75.0 


15 


0.380E+03 


0. 213E+05 


-0 . 166E+02 


0. 238E+04 


-0. 407E+01 


0. 166E+02 


75.0 


16 


-0. 189E+03 


0. 397E+05 


-0. 200E+02 


-0.24 1E+04 


-0. 390E+01 


0 . 1 53E+02 


90.0 


17 


0. 166E+02 


0. 378E+05 


-0. 109E+02 


-0. 223E+04 


-0. 117E+01 


-0.141E+01 


90.0 


18 


0 . 627E+03 


0. 208E+05 


-0. 882E+01 


0. 250E+04 


-0. 172E+01 


-0 . 371E+01 


‘ 75.0 


19 


0 . 627E+03 


0. 208E+05 


-0. B82E+01 


0. 250E+04 


0. 172E+01 


0 . 371E+01 


75. 0 


20 


0. 166E+02 


0. 378E+05 


-0. 109E+02 


-0. 223E+04 


0.117E+01 


0. 141E+01 


90.0 


21 


-0. 564E+03 


0. 377E+05 


0 . 34 7E+02 


-0. 209E+04 


0. 981E+00 


-0. 526E+01 


90. 0 


22 


0.41 1E+02 


0. 209E+05 


0. 344E+02 


0 . 260E+04 


0. 167E+00 


- 0 . 618E+01 


75.0 


23 


0. 411E+02 


0. 209E+05 


0.344E+02 


0. 260E+04 


-0. 167E+00 


0. 61BE+01 


75.0 


24 


- 0 . 564E+03 


0 . 377E+05 


0. 347E+02 


-0. 209E+04 


-0 . 981E+00 


0 . 526E+01 


90.0 


25 


-0. 293E+03 


0. 375E+05 


-0.1 01E+03 


-0.243E+04 


-0 . 309E+01 


-0 . 679E+01 


90.0 
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26 


0. 146E+03 


0 . 192E+05 


-0 . 844E+02 


0 . 167E+04 


-0. 457E+01 


- 0 . 428E+01 


75.0 


27 


0. 146E+03 


0. 192E+05 


-0. 844E+02 


0 . 187E+04 


0.457E+01 


0.428E+01 


75.0 


28 


-0 . 293E+03 


0 . 375E+05 


-0. 101E+03 


-0.24 3E+04 


0.309E+01 


0 . 679E+01 


90.0 


29 


0. 120E+02 


0.418E+05 


-0. 140E+03 


- 0 . 252E+04 


-0 . 267E+01 


-0 . 4 87E+01 


90.0 


30 


0 . 664E+03 


0 . 227E+05 


-0. 122E+03 


0 . 265E+04 


-0 . 502E+01 


-0 . 298E+01 


75.0 


31 


0 . 664E+03 


0 . 227E+05 


-0 . 122E+03 


0 . 265E+04 


0 . 502E+01 


0 . 298E+01 


75.0 


32 


0. 120E+02 


0.41 8E+05 


-0. 140E+03 


-0 . 252E+04 


0 . 267E+01 


0 . 4 87E+01 


90.0 


1 


- 0 . 1 23E+03 


0 . 308E+05 


-0 . 208E+02 


-0 . 538E+03 


-0 . 223E+01 


0 . 212E+01 


90.0 


2 


-0. 127E+03 


0.308E+05 


-0 . 350E+02 


-0 . 538E+03 


-0 . 696E+00 


0 . 1 55E+01 


90.0 


3 


- 0 . 1 27E+03 


0 . 308E+05 


-0. 350E+02 


-0 . 538E+03 


0 . 696E+0 0 


-0 . 155E+01 


90.0 


4 


-0 . 123E+03 


0. 308E+05 


-0. 208E+02 


-0. 538E+03 


0 . 223E+0 1 


-0 . 212E+01 


90.0 


5 


- 0 . 874E+02 


0. 301E+05 


0 . 748E+02 


- 0 . 690E+03 


-0. 302E+01 


0 . 302E+01 


90.0 


6 


- 0 . 91 3E+02 


0 . 3 0 1E+05 


0 . 613E+02 


-0 . 690E+03 


-0. 329E+00 


0 . 357E+01 


90.0 


7 


-0.913E+02 


0. 301E+05 


0 . 613E+02 


-0 . 690E+03 


0. 329E+00 


- 0 . 357E+01 


90.0 


8 


-0. 874E+02 


0. 301E+05 


0 . 748E+02 


-0 . 690E+03 


0. 302E+01 


-0. 302E+01 


90.0 


9 


-0 . 698E+02 


0. 300E+05 


0. 499E+02 


-0.775E+03 


-0. 306E+01 


0. 358E+01 


90.0 


10 


-0 . 728E+G2 


0 . 3 O0E+O5 


0 . 3 93E+02 


-0. 775E+03 


-0 . 124E+01 


0 . 4 57E+01 


90.0 


11 


-0. 728E+02 


0.300E+05 


0. 393E+C2 


-0 . 775E+03 


0.124E+01 


-0.457E+01 


90.0 


12 


-0 . 698E+02 


0 . 300E+C5 


0. 499E+02 


- 0 . 775E+03 


0. 306E+01 


-0.358E+01 


90.0 


13 


-0.488E+G2 


0 . 2S5E+05 


-0. 410E+02 


- 0 . 81 0E+03 


0. 150E+01 


- 0 . 539E+01 


90.0 


14 


-0. 517E+02 


0. 2S9E+05 


-0.514E+02 


-0. 810E+03 


0. 943E+00 


-0 . 355E+01 


90.0 


15 


-0. 517E+C2 


0. 299E+05 


-0 . 514E+02 


-0. 810E+03 


-0. 943E+00 


0 . 3 55E+01 


90.0 


16 


-0 . 488E+02 


0 . 299E+05 


-0 . 410E+02 


- 0 . 810E+03 


-0. 150E+01 


0.539E+01 


90.0 


17 


-0. 822E+02 


C . 302E+05 


-0. 545E+02 


- 0 . 768E+03 


0 . 1 78E+01 


-0. 117E+02 


90. 0 


18 


-0. 849E+02 


0. 301E+05 


-0 . 640E+02 


-0. 769E+03 


0. 309E+00 


-0 . 996E+01 


90.0 


19 


-0. 849E+02 


0 . 301E+05 


- 0 . 64 OE+02 


-0. 769E+03 


-0 . 309E+00 


0 . 996E+01 


90.0 


20 


-0 . 822E+02 


0. 302E+05 


-0.54 5E+02 


-0. 768E+03 


-0. 178E+01 


0 . 117E+02 


90.0 


21 


-0. 128E+03 


0. 309E+05 


0 . 374E+02 


- 0 . 668E+03 


0. 724E+00 


- 0 . 125E+02 


90.0 


22 


-0. 130E+03 


0. 309E+05 


0. 298E+02 


- 0 . 668E+03 


0 . 211E+00 


- 0 . 1 16E+02 


90.0 


23 


- 0 . 1 30E+03 


0. 309E+05 


0 . 298E+02 


- 0 . 668E+03 


-0. 211E+00 


0 . 1 16E+02 


90.0 
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24 


-0. 128E+03 


0. 309E+05 


0 . 374E+02 


25 


-0. 162E+03 


0. 315E+05 


0 . 563E+02 


26 


-0. 165E+03 


0. 315E+05 


0 . 456E+02 


27 


-0.165E+03 


0 . 315E+05 


0 . 4 56E+02 


28 


-0. 162E+03 


0 . 315E+05 


0 . 563E+02 


29 


-0. 163E+03 


0 . 314E+05 


-0 . 605E+02 


30 


-0 . 166E+03 


0 . 314E+05 


-0 . 699E+02 


31 


-0. 166E+03 


0. 314E+05 


-0 . 699E+02 


32 


-0 . 163E+03 


0 . 314E+05 


-0 . 60 5E+02 


NEL 


SI 


S2 


S3 


1 


0.371E+05 


-0 . 336E+02 


-0 . 470E+00 


2 


0. 216E+05 


0 . 612E+03 


-0 . 521E+00 


3 


0 . 216E+05 


0 . 612E+03 


-0. 521E+00 


4 


0. 371E+05 


-0.336E+02 


-0. 470E+00 


5 


0. 382E+05 


-0. 722E+03 


0. 137E+02 


6 


0. 214E+05 


-0. 249E+02 


0. 158E+02 


7 


0. 214E+05 


-0. 249E+02 


0. 158E+02 


6 


0 . 382E+05 


-0. 722E+03 


0. 137E+02 


9 


0. 393E+05 


-0. 530E+03 


0. 120E+02 


10 


0 . 211E+05 


0. 229E+03 


0. 137E+02 


11 


0. 211E+05 


0 . 229E+03 


0. 137E+02 


12 


0. 393E^-05 


-0 . 530E+03 


0. 120E+02 


13 


0. 397E+05 


-0. 189E+03 


-0. 200E+02 


14 


0 . 21 1E+05 


0 . 589E+03 


-0. 166E+02 


15 


0 . 211E+05 


0. 589E+03 


-0. 166E+02 


16 


0. 397E+05 


-0. 189E+03 


-0.200E+02 


17 


0.37 8E+05 


0. 166E+02 


-0.109E+02 


18 


0 . 2O7E+05 


0. 730E+03 


“0 . 882E+01 


19 


0. 2O7E+05 


0. 730E+03 


-0. 882E+01 



. 668E+ 03 


-0. 724E+00 


0. 125E+02 


90.0 


. 620E+03 


0. 193E+01 


-0.617E+01 


90.0 


. 620E+03 


0. 117E+01 


-0.485E+01 


90.0 


. 620E+03 


-0.117E+01 


0 . 485E+01 


90.0 


. 620E+03 


-0. 193E+01 


0 . 617E+01 


90.0 


. 589E+03 


-0.297E+01 


0 . 613E+01 


90.0 


. 589E+03 


-0.114E+01 


0 . 337E+01 


90.0 


. 589E+03 


0 . 114E+01 


-0 . 337E+01 


90.0 


. 589E+03 


0.297E+01 


- 0 . 613E+01 


90.0 


S6 


S5 


S4 




. 197E+04 


0 . 316E+01 


-0. 142E+01 




. 266E+04 


0 . 345E+01 


-0.451E+00 




. 266E+04 


-0 . 34 5E+01 


0.451E+00 




. 197E+04 


-0.316S+01 


0 . 142E+01 




. 203E+04 


0. 284E+01 


0 . 1 23E+01 




. 298E+04 


0.47 0E+01 


0.195E+01 




. 298E+04 


-0.47 0E+01 


-0 . 195E+01 




. 203E+04 


-0. 284E+01 


-0. 123E+01 




. 228E+04 


-0. 160E+02 


-0. 801E+00 




. 3 16E+04 


-0.148E+02 


-0 . 845E+01 




. 316E+04 


0. 148E+02 


0 . 845E+01 




. 228E+04 


0. 160E+02 


0.801E+00 




.241E+04 


-0. 153E+02 


-C . 390E+01 




. 3 16E+04 


-0. 150E+02 


-0 . 823E+01 




. 316E+04 


0. 150E+02 


0 . 823E+01 




. 24 1E+04 


0. 153E+02 


0 . 390E+01 




. 223E+04 


-0. 141E+01 


0.117E+01 




. 289E+04 


-0.403E+01 


0 . 7 04E+00 




. 289E+04 


0.403E+01 


-0. 704E+00 





■0 

0 

•0 

0 

•0 

0 

0 

■0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 



0. 378E+05 
0.377E+05 - 
0. 208E+05 
0 . 208E+05 
0.37 7E+05 - 
0 . 375E+05 - 
0. 189E+05 
0 . 189E+05 
0 . 375E+05 - 
0 . 4 1BE+05 
0. 226E+05 
0 . 226E+05 
0 . 4 18E+05 
0 . 3 08E+05 - 
0 . 308E+05 - 
0 . 308E+05 - 
0 . 308E+05 - 
0 . 3 01E+05 - 
0 . 3 01E+05 - 
0 . 3 01E+05 - 
0 . 3 0 1E+05 - 
0.30 CE+05 - 
0.30 0E+05 - 
0. 300E+05 - 
0 . 300E+05 - 
0 . 299E+05 - 
0 . 299E+05 - 
0 . 299E+05 - 
0 . 299E+05 - 
0 . 3 02E+05 - 



0. 166E+02 
0 . 564E+03 
0. 140E+03 
0. 140E+03 
0 . 564E+03 
0 . 293E+03 
0.488E+03 
0 . 48BE+03 
0 . 293E+03 
0. 120E+02 
0.B1BE+03 
0 . B18E+03 
0. 120E+02 
0. 123E+03 
0 . 1 27E+03 
0 . 127E+03 
0. 123E+03 
0 . 87 4E+02 
0. 913E+02 
0 . 913E+02 
0. B74E+02 
0.69BE+02 
0 . 728E+02 
0. 72BE+02 
0 . 69BE+02 
0 . 4 B8E+02 
0. 517E+02 
0. 517E+02 
0 . 4 8BE+02 
0 . 822E+02 



-0.109E+02 
0. 347E+02 
0 . 34 4E+02 
0 . 344E+02 
0. 347E+02 
-0. 101E+03 
-0 . B44E+02 
-0. 844E+02 
-0. 101E+03 
-0 . 140E+03 
-0. 122E+03 
-0 . 1 22E X 0 3 
-0. 140E+03 
-0. 20BE+02 
-0 . 250E+02 
-0 . 350E+02 
-0. 208E+02 
0. 74BE+02 
0.613E+02 
0 . 613E+02 
0. 74BE+02 
0 . 499E+02 
0 . 393E+02 
0 . 393E+02 
0. 499E+02 
-0 .^4 1 0E+0 2 
-0. 514E+02 
-0. 514E+02 
-0. 410E+02 
-0. 545E+02 



0 . 223E+04 
0 . 209E+04 
0. 297E+04 
0 . 297E+04 
0 . 209E+04 
0 . 243E+04 
0.315E+04 
0 . 315E+04 
0.243E+04 
0. 252E+04 
0 . 323E+04 
0. 323E+04 
0. 252E+04 
0. 536E+03 
0. 538E+03 
0 . 53BE+03 
0. 53BE+03 
0. 690E+03 
0 . 690E+03 
0 . 690E+03 
0 . 690E+03 
0. 775E+03 
0. 775E+03 
0. 775E+03 
0. 775E+03 
0.810E+03 
0. 810E+03 
0 . 810E+03 
0. B10E+03 
0 . 76BE+03 



0.141E+01 
-0. 526E+01 
-0 . 592E+01 
0 . 592E+01 
0 . 526E+01 
-0.679E+01 
-0.532E+01 
0 . 532E+01 
0 . 679E+01 
-0 . 487E+01 
-0.41BE+01 
0 . 4 18E+01 
0 . 4B7E+01 
0. 212E+01 
0. 155E+01 
-0. 155E+01 
-0. 212E+01 
0 . 302E+01 
0 . 357E+0 1 
-0 . 357E+0 1 
-0. 302E+01 
0 . 356E+01 
0. 457E+01 
-0. 457E+01 
-0. 35BE+01 
-0. 539E+01 
-0. 355E+01 
0 . 355E+01 
0 . 539E+01 
-0.117E+02 



-0.117E+01 
-0 . 9B1E+00 
-0.176E+01 
0 . 1 76E+01 
0 . 981E+00 
0 . 309E+01 
0 . 331E+01 
-0 . 331E+01 
-0. 309E+01 
0 . 267E+01 
0 . 407E+01 
-0 . 4 07E+0 1 
-0. 267E+01 
0 . 223E+01 
0 . 696E+00 
-0.696E+00 
-0 . 223E+01 
0. 302E+01 
0 . 329E+00 
-0 . 329E+00 
-0.302E+01 
0. 306E+01 
0. 124E+01 
-0.124E+01 
-0. 306E+01 
-0. 150E+01 
-0. 943E+00 
0 . 943E+00 
0.150E+01 
-0 . 178E+01 
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ALFA 



50 


0. 301E+05 


-0. B49E+02 


-0.64 0E+02 


0 . 769E+03 


-0.996E+01 


-0 . 309E+00 




51 


0. 301E+05 


-0 . 849E+02 


-0. 640E+02 


0. 769E+03 


0. 996E+01 


0 . 309E+00 




52 


0 . 302E+05 


-0 . B22E+02 


-0 . 545E+02 


0 . 768E+03 


0 . 117E+02 


0 . 178E+01 




53 


0 . 309E+05 


-0. 128E+03 


0 . 374E+02 


0 . 668E+03 


-0. 125E+02 


-0.724E+00 




54 


0. 309E+05 


-0. 130E+03 


0 . 298E+02 


0 . 668E+03 


-0 . 1 16E+02 


-0. 211E+00 




55 


0 . 309E+05 


-0. 130E+03 


0 . 298E+02 


0 . 668E+03 


0 . 116E+02 


0. 211E+00 




56 


0. 309E+05 


- 0 . 128E+03 


0 . 374E+02 


0 . 668E+03 


0. 125E+02 


0. 724E+00 




57 


0 . 315E+05 
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